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Supervisors' Foreword 



The main physics event of 2012 was the announcement on July the 4th of the 
discovery of a new boson at CERN with a mass close to 126 GeV. 

The on-going analysis of the properties of the particle discovered last year 
points with increasing evidence to the long awaited for particle related to the 
Brout-Englert-Higgs mechanism introduced in the 1960s. The main focus of the 
investigations at the LHC experiments now shifts to the possibility of finding other 
high-energy phenomena such as supersymmetry or extra dimensions, in particular 
in the forthcoming run at 13 TeV energy in 2015. While evidence of new physics 
at the TeV scale is mounting through astrophysical/cosmological observations and 
neutrino experiments, the challenges for a direct observation at the LHC are 
formidable. The discovery of the new boson is a major breakthrough which was 
also based on how precise knowledge and predictions of its properties in the 
Standard Model. Yet it took almost 50 years. Looking for new physics is much 
more difficult as options are wide open and one in many, even unexpected pos- 
sibilities could turn to be the one chosen by nature. That is the challenge at the 
basis of this thesis which, starting from one of such class of models, faces the 
problem of how to search for it via accurate and flexible computer simulations. 

After an introduction to the physics of the Standard Model and beyond, this 
work presents the methods of investigation of the LHC data in the context of 
extensions of the Standard Model featuring gravitons and gravitinos. The pro- 
duction of these particles in association with jets is studied as one of the most 
promising signatures to discover new physics at the LHC. In this work, Priscila 
uses the most advanced techniques, such as algorithms allowing the computation 
of Feynman graphs and helicity amplitudes, to help in the development of new 
simulation tools specialized for graviton/gravitino-like particles. These tools are 
subsequently employed revealing its power. 
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Supervisors' Foreword 



For interpreting the data from the LHC, such tools are currently being used and 
have shown to be crucial for the characterization of the properties of the new 
boson. They will become even more relevant if new physics will be detected at the 
LHC. Any student in theoretical physics as well as experimentalist in high-energy 
physics will find in this work a lively description of the strategies as well as of the 
cutting-edge technologies employed in the search for new physics in collider 
experiments. 
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Abstract 



The Standard Model of fundamental particles and their interactions is one of the 
most successful theories in physics. In particular, up to the weak scale (a few 
hundreds of GeV) it agrees to a great degree with a large set of experimental data. 
However, there are several theoretical reasons, such as the so-called "Hierarchy 
Problem", as well as experimental ones, such as the neutrino masses and the 
evidence for dark matter in the Universe, to expect that something new (particles 
and/or interactions) could lie at TeV scale. Hints and/or answer(s) to these fun- 
damental questions will be provided by the Large Hadron Collider (LHC), 
a proton-proton collider running at high energies. 

The present thesis aims to explore new physics at the LHC through phenome- 
nological studies that employ simulations and computational tools to directly link 
theories with experimental data. In particular, the focus is on Beyond Standard 
Model theories that can incorporate a quantum description of gravity, such as extra 
dimensional theories, Supersymmetry and yet a simplified model where the constant 
of Newton is scale dependent. Phenomenological analyses are performed in which 
graviton and gravitino emission in combination with multiple jets are investigated at 
hadron colliders. Inclusive samples are generated by merging matrix element with 
parton shower descriptions, and validated by a comparison against next-to-leading 
order QCD calculations. Predictions for relevant observables at the LHC and 
Tevatron are obtained. 
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Chapter 1 
Introduction 



The Standard Model (SM) of fundamental particles and their interactions is one of the 
most successful theories in physics. It provides an astonishingly accurate description 
of phenomena in a wide range of scales. In particular, up to the weak scale (a few 
hundreds of GeV) it agrees to a great degree with the experimental data collected, 
also predicting the existence of a new scalar state, the Higgs particle. 1 
There are several theoretical reasons, such as the so-called "Hierarchy Problem" 
and the "fermion mass hierarchy", as well as experimental ones, such as the neutrino 
masses and the evidence for dark matter in the Universe, to expect that "new physics" 
could lie at the TeV scale. Hints and/or answer(s) to these fundamental questions are 
expected to be provided very soon by the Large Hadron Collider (LHC), a proton- 
proton collider running at high-energies. 

In particular, the search for a quantum field theory aiming at unifying gravity with the 
other three known forces of nature (electromagnetic, weak and strong forces) has been 
the main subject of several researches for many years. If the SM is not the complete 
picture, it may be possible to construct effective theories that incorporate gravity as a 
field theory. Higher dimensional theories have been an attempt of such scenarios since 
the beginning of the twentieth century, when it was realized that the 5 -dimensional 
generalization of the theory proposed by Einstein could describe both gravitational 
and electromagnetic interactions. However, because of the weak characteristic of 
gravity, inconsistencies arose when the mathematical theory was formulated. Such 
theories were therefore left aside until the appearance of the theories of Supergravity 
and Superstrings in the early 1970s. 

More recent ideas using extra-dimensions have been proposed in the context of 
minimal extensions of the SM. These theories are created to address SM limitations, 
and have a particular behaviour at high-energies. They however maintain all the low- 
energy features of the SM. The philosophy of such extra-dimensional theories lies 
in the idea that the SM is a 4-dimensional projection of a theory living in higher- 
dimensional space-time. Different scenarios have been built with such principles. We 



The most recent results from the LHC and Tevatron experiments provide, however, the first 
experimental indications of its existence. 

P. de Aquino, Beyond Standard Model Phenomenology at the LHC, 1 
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can cite for instance theories with large extra dimensions [1—3], and theories with 
warped extra dimensions [4, 5]. 

Furthermore, there are different extensions to the theory that are based on more effec- 
tive approaches to include gravity. For example, general descriptions of 4-dimensional 
theories can account for gravitational interactions with SM particles through the exis- 
tence of new particles living in a hidden sector [6—11]. 

In addition, further extensions can be cited, such as Supersymmetry, which is based 
on symmetry principles rather than on extra-dimensions. Supersymmetric theories 
start from the assumption that there must exist an extra symmetry that transforms 
fermions into bosons, and vice-versa. As a consequence, all particles of the SM 
must have super-partners, i.e. particles that differ from SM particles by half a spin. 
Variations of such theories can predict, for example, the existence of the graviton 
particle, a hypothetical spin-2 particle that mediates the force of gravitation, and its 
super-partner, the spin-3/2 particle named gravitino. 

The observation of new fundamental particles at hadron colliders can therefore be 
one indication of the existence of Beyond Standard Model (BSM) theories. Typically, 
each new physics model has its own characteristic predictions. Yet, it is often the case 
that very different theories might lead to similar signatures at colliders. This is the 
case, for example, with models that predict the existence of a dark matter candidate, 
weakly interacting, whose emission would just generate missing transverse energy in 
energetic multi-jet events. Gravitinos and neutralinos, predicted by supersymmetric 
theories, or gravitons, in specific variations of extra-dimensional models, are just a 
few examples of particles and models giving possibly equivalent signatures at the 
LHC. 

With the LHC running, and a high luminosity being achieved, chances are that new 
physics will be emerging soon. However, the fact that so far no BSM particle has been 
detected, implies that more careful, accurate and dedicated studies will be needed. 
If evidence for new physics will emerge, it will be urgent yet certainly not trivial to 
understand its nature and to clearly identify it. The typical final states produced at 
the LHC involve a large number of jets, heavy-flavour quarks, missing energy and 
leptons. Therefore, such an overwhelming number of possible signs of new physics 
need a very careful understanding of the detector as well as an accurate modelling 
of the data themselves. 

For this reason, precise simulation and computational tools play a key role in the 
current analyses. They have become indispensable in the effort of connecting theory 
with experimental data sets. 

The possibility of automatically obtaining predictions at colliders from any BSM 
theory shows substantial progress. With such, it is possible to test and validate theo- 
ries before an experimental analysis is performed. Tests and validations account for 
the specifications of the experiment and are very time consuming. This step already 
considerably reduces the number of theories suitable for correctly describing exper- 
imental data once these are at our disposal. 

Until a few years ago, such possibility was not available. The implementation of 
complex BSM models in existing general purpose event generators was often a 
process susceptible to several calculation errors. That was because the derivation 
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of Feynman rules to describe new interactions and their implementation into event 
generators depended on several aspects, going from convention definitions to com- 
putational difficulty. This situation has dramatically changed with the introduction 
of new simulation tools that can automatically link arbitrary lagrangians to a par- 
ticular Monte-Carlo event generation that can be chosen by the user. Moreover, the 
development of several alternative simulation techniques (such as the merging proce- 
dure employed to match simulations performed by matrix element event generators 
with those obtained from the use of parton shower description) have substantially 
improved the precision of predictions at high energy hadron colliders. The evolution 
of such tools now provides a user-friendly manner to simulate BSM theories, and 
does not require a deep knowledge and/or strong computer skills to be applied. 
The work presented in this thesis aims at performing phenomenological investiga- 
tions at the LHC particularly focused on theories that incorporate gravity by means 
of a set of simulation tools and techniques developed for this purpose. We demon- 
strate how it is possible to perform complete phenomenological investigations that 
cover all the steps to link any theory with experiments. Furthermore, we explicitly 
show the particular efforts we have made for the development of automatic tools. As 
a natural step, we make use of such tools and design search strategies that can be 
applied in different studies, including those at a more experimental level. 
This thesis is divided in five main chapters. In the following chapter, we present an 
introduction to the Standard Model and why extensions thereof are considered. In 
Chap. 3, gravitation as a quantum field theory is introduced as well as BSM theories 
that include the graviton particle. Theories with extra-dimensions and a particular 
theory, in which Newton's constant vary according to the effective scale of the theory, 
are considered. The Large Hadron Collider is introduced in Chap. 4, together with 
simulation techniques that will be employed through the present study. In Chap. 5, 
we present the main phenomenological analyses for graviton production at hadron 
colliders. We consider the emission of one graviton in combination with jets. Finally, 
in Chap. 6, we introduce Supersymmetry and perform phenomenological investiga- 
tions. Despite the large number of new particles predicted within this theory, we 
focus on the production of gravitino in combination with jets. Similar techniques to 
those applied for graviton production analyses are employed for the gravitino. We 
conclude with a complete comparison of all the results obtained, and an outlook. 
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Chapter 2 

The Standard Model 



The Standard Model (SM) is a successful description of particle physics nowadays. 
It is a quantum field theory that explains the dynamics of our universe through matter 
and forces. It incorporates three fundamental forces: the weak, the strong and the 
electromagnetic forces, and includes different types of particles. 
Mainly, the SM states that all the matter is composed by fermionic fields interacting 
among each other via vector fields. In order to understand the field content of the 
SM, it is necessary to begin with the definition of the symmetry group. 



2.1 Symmetries, Gauge Theory and Particle Content 

Symmetries have always played an important role in physics. The creation of the SM 

has not followed a different path: it is a theory based on a local symmetry. 

The first theorem of Noether assures that any differentiable symmetry of the action 

of a physical system leads to a corresponding conservation law [ 1 ] . Any conservation 

law of physics can be interpreted as resulting from the symmetries of a particular 

theory. 

One example is the theory of Quantum Electrodynamics (or, QED). The invariance 
under local transformation implies the existence of gauge fields with specific prop- 
erties. To demonstrate such remarkable consequence, let us begin by considering a 
free theory of fermions. The free lagrangian can be written as 

C = WY ll d* -m)f, (2.1) 

with \jr being the spin-1/2 field, %[r — 4r^y°, and y^ the Dirac matrices. The 
lagrangian is invariant under a global gauge transformation \[/ — >• e~""4r, with a 
a constant phase, which implies the conservation of the Dirac current j^{x): 

9„/=0; j^fy^f. (2.2) 
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However, Eq. (2.1) is not invariant under local gauge transformations, where now a 
depends on the space-time coordinates, and suggests that the derivative should be 
redefined as 

9 M _^ D n = d ,x + iqA n (2.3) 

where A 11 is a vector field. The lagrangian (2.1) with the replacement 3^ — > D >1 is 
now invariant under the following local transformations 

f{x) -» fix) = e- ia{x) fix), (2.4) 
A^ix) -+ A* 4 ' = A^ix) + -d^aix), (2.5) 

q 

and is finally given by 

£=fiiy li d tl -m)ir-q4ry ll A ll f. (2.6) 



The second term has been derived from the imposition of local invariance. It describes 
the interaction of the fermionic matter f already existent on the theory, with the gauge 
field A' 1 . Therefore, a theory that had in principle only matter fields, needs vector 
fields to provide interactions amongst fermions. 

One could thus generalize this principle for all interactions, i.e. generate interaction 
terms for the weak, electromagnetic, strong (and later gravitational) forces through 
specific symmetries imposed on the theory. This is basically the idea for constructing 
the SM initially proposed by Glashow [2] and independently by Salam and Ward [3], 
extended later by Weinberg [4] and Salam [5]. Of course the complexity is larger 
than previously explained, since it works perfectly for abelian theories, however the 
world is not only described by them. Before going ahead, it is worth reminding that 
the SM of particle physics was created as a puzzle, and each piece was put together 
at different moments of History. 

Initially, the electroweak theory was developed, and the gauge group SU (2)l x 
U (l)y was used to relate electric charge with the isospin and the leptonic hypercharge 
of a particle. Besides matter fields, four gauge bosons are included due to the gauge 
symmetry: one triplet associated with the SU (2) group and one singlet associated to 
the Uil). The corresponding lagrangian, restricted to the leptonic fields, is given by 

C = Ly^D^L + e R y^e R - -W""'w£„ - -B^B^ (2.7) 
where the covariant derivatives are defined by 

D u = d u + i -a' Wl +i—B u , 

M 2 m^ 2 m (2.8) 

D^ = d„ + ig'B IM 



2.1 Symmetries, Gauge Theory and Particle Content 



7 



and L is the isospin doublet that contains the left-handed neutrino and electron, e« 
is the right-handed electron, 1 a' are the Pauli matrices, and g and g' are the coupling 
constants. W^ v and B^ v are the field-strength tensors: 

B^v — d/iB v — d v B (1 . 

The matter content of the world is not only formed by leptons but also by hadrons, 
and it has been discovered that hadrons are composed by quarks. The symmetry that 
relates the color charge in quarks is SU (3) f . The SM is consequently defined to 
have SU (3) c x SU (2)l x t/(l)y as the gauged symmetry group. Each gauge boson 
field is associated to the generator of the algebra of each group. Therefore, eight 
coloured spin-1 particles associated to the SU C (3) gauge group exist and are also 
known as gluons. In addition, there are still the four un-coloured particles, W' and 
B^v quoted previously, that will mix to form the massive W ± and Z° gauge bosons 
and the massless photon. 

In a nutshell, the final particle content of the SM is summarized in Fig. 2. 1 [6]. The 
masses and most significant quantum numbers of each particle are stated in the figure. 
One important piece of the SM, the scalar field, has not yet been mentioned. It has 
been stated that conservation laws are a consequence of imposed symmetries. The SM 
theory is therefore based on gauge symmetries, and lagrangians are constructed from 
the assumption of massless fields. The mechanism to generate mass to the particles 
requires the existence of at least one scalar particle, the Brout-Englert-Higgs boson. 2 
This is the so-called BEH mechanism [7-9], which will be explained in what follows. 



Fig. 2.1 An illustration of 
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Standard Model, excluding 
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Right-handed neutrinos are not included in the theory. 
2 From the present moment on we shall abbreviate Brout-Englert-Higgs boson by Higgs boson, 
which is the usual shortened form used in the literature. 
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2.2 The Brout-Englert-Higgs Mechanism 

It is known from Noether's theorem that symmetries imply conservation laws. It is 
a consequence of the invariance of both the lagrangian and the vacuum of a theory. 
However, there could be a situation where the lagrangian is invariant under a symme- 
try in which the vacuum is not. If such circumstance occurs, the symmetry is defined 
as to be broken, or hidden. 

Consider for instance a scalar field (p whose lagrangian is given by [10]: 

£0 = \{^4>W4>) - V(4>), with v(4>) = X -^ 2 <p 2 + l -x<p\ (2.10) 

where (p is real and A > 0. The lagrangian is clearly invariant under the transformation 
(p — >• —(p. When the vacuum expected value (vev) is calculated, two separated 
situations arise: 

(a) If /x 2 > 0, the vacuum is invariant under such symmetry, 

(cP) = (O|0|O> = ( M 2 >0). (2.11) 

(b) However, if fj? < 0: 



<0)o = ±y = ±j= (M 2 < 0). (2.12) 

The vacuum state in this case is degenerated, and it depends on the choice 
between +v and —v. 

Both situations are illustrated in Fig. 2.2. On the left side, the potential is shown as a 
function <j> for fj? > 0. As Eq. 2.1 1 revealed, only one vev is obtained. However, on 
the right, the plot of same potential is shown for the choice of n 2 < 0. Equation 2.12 
exhibits a degenerate vacuum state, displayed by the two local minimums in Fig. 2.2b. 




Fig. 2.2 Potential V(<j>) of Eq. 2.10 as a function of <j> for the two situations presented: a with 
jj? > and b for fj. 2 < 0. For case (a), the minimum is at 0o = and for case (b) (j>o = ±u, where 
the vacuum is degenerated 
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We could therefore choose one state, for instance (0) o = +v, and re-define the field 
in order for the vacuum to be at the origin: 

|(jc) = 4>(x) - Wo = Hx) - v. (2.13) 

Now (|>o = 0, and the lagrangian (2.10) becomes 

A = ^OmDO^) - 2 - ^ 3 - \^ 4 - (2-14) 

This is the lagrangian of a free scalar field £ with mass term 3 m$. — v '—2p 2 . The sym- 
metry was therefore hidden behind the primary definition of the field. By redefining 
the field it was possible to obtain the massive £ . This concept is named spontaneous 
symmetry breaking and it illustrates the main mechanism used to generate the mass 
of fermionic and bosonic fields of the SM. 

In the SM, it is necessary to define a doublet composed of two complex scalar fields. 
Suppose the doublet can be defined in a Hermitian basis as [11-13] 

0= Uv = viU-W (2 - 15) 

The scalar lagrangian to be added to (2.7) will be similar to the one given in (2.10). 
However, in order to maintain the gauge invariance under SU(2)i x U(\)y, the 
covariant derivative (2.8) should be used: 

Cbeh = \D„4>\ 2 - V(<j>). (2.16) 

In the new basis, the potential is 



V(4>) 



^ 2 fe>A+^fe>?) • (2-17) 



We can choose the position of the minimum as (f>\ — 02 = 04 = and 03 = v, 
being v the Higgs vev. A new field h can be defined to be a radial excitation around the 
vev. The symmetry SU ( 2)i x U ( l)y is broken to U (l) em if we choose the specific 
value for the vev to be *J—(i 2 /2X. The field can thus be re-written as 

(218) 

The potential V (0) of the lagrangian becomes 



3 Notice the mass is positive, since we are analysing the case where fi 2 < 0. 
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Vbeh = —kv 1 * + kv^ + kvh 3 + -Xh*. (2.19) 

The second term indicates the mass of the Higgs field, m 2 h — 2v 2 X. The third and the 
forth terms, i.e. the terms related to h 3 and h A , indicate self-couplings of the Higgs 
boson. 

The kinetic term of the lagrangian, | <f> \ 2 includes terms with the W'^ and fields, 

+ *| = ^T^~ [ff 2 (<) 2 + 9 2 (W 2 ) 2 + {-gWl + 5 %) 2 ] . 

(2.20) 

If four new vector fields are defined as: 

Z° = - 1 = (gWf. - g'Bg) , W^r — -7= (W l =p i W 2 ) , 

(2.21) 



fj 



Equation 2.20 can be written in terms of the new fields. The boson masses can be 
identified by the following terms of the lagrangian 

m 2 w W+W-> i + - {m 2 Z„Z^ + m^A") . (2.22) 

With the above expression, it is straightforward to find the masses related to the W jjL , 
Z M and fields [11]: 

m w = m z - + g' 2 and m A = 0. (2.23) 

Notice however that covariant derivatives can now be written in terms of the new 
bosonic fields. Consequently, it is possible to re- write Eq. 2.8 in a more useful format, 
in terms of the electromagnetic coupling, given by 



99' 



J<F+g 



'2 



(2.24) 



or in terms of the weak mixing angle, 



9 9 
g e — gsin6 w , cos9 w — — , and sin^u, = — (2.25) 

Vg 2 + g' 2 Vg 2 + g' 2 

The mass of the W ± and Z bosons are therefore related through the expression 

m w — mz cos 6 W . (2.26) 
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As a result, all effects of W ± and Z exchange processes (at tree level) can be exhibited 
as a function of the g e , 9 W and mw parameters. 

Finally, to find the value for the vev v we can replace the expression found for the 
mass of the W ± bosons in terms of the Fermi constant: 



G F = V2-^ = ^. (2.27) 
8 mi, 2 v 2 



Because the Fermi constant is experimentally known with a very good accuracy, 
G F = 1.16637(1) x 1CT 5 GeV" 2 [14], the value for the Hi ggs vev can be deduced 
as v ~ 246 GeV. This is the value where the SU (2)l x U (l)y is broken. How- 
ever, because the parameter X is not known, the mass of the Higgs boson cannot be 
predicted. 

A similar mechanism is used to obtain the mass of the fermions. The invariant 
lagrangian that should be added to (2.7) can be exhibited as 

£yuk — L(j> e R - \ d Ql <t>dR - X u Ql4>ur + h.c. (2.28) 

where Qi is the isospin doublet that contains the left-handed up and down quarks, 
ur and cIr are the right-handed up and down quarks, (j> — io2<p* , and ai is one of 
the Pauli matrices. After spontaneous symmetry breaking, the mass of the fermions 
(except neutrinos) are generated as m f = k/ v/2, and neutrinos continue to be 
massless. The couplings Xf are called Yukawa couplings [11], and are determined 
depending on the experimental values of fermionic masses. 



2.3 Experimental Facts and Standard Model Constraints 

As it has been stated, the SM is one of the most tested theories in physics. It has 
been experimentally investigated through different methods, and nowadays particle 
colliders are used for such task. 

The final missing piece is, with no doubt, the experimental observation of the scalar 
sector of the SM through the Higgs particle. Very recent developments have been 
made in order to limit and constrain the mass of the Higgs particle using data from the 
Large Hadron Collider (LHC), the revolutionary proton-proton collider experiment 
situated in Switzerland. At this stage, we limit ourselves to display the most advanced 
SM results obtained from LHC data. A detailed description related to such collider 
will be presented in Chap. 4. 

The LHC has started colliding protons for physics analyses in 2009, breaking energy 
records ever since. For more than a year starting in 2010, it has collided protons with 
a center-of-mass energy of 7 TeV, achieving a considerably high sensitivity. For this 
reason, many of the SM parameters have already been very precisely measured by 
the LHC. A summary of the most recent values due to electroweak constraints can be, 
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Fig. 2.3 Higgs branching ratio for contributing decay channels as a function of the Higgs mass. 
These plots are made by the Higgs cross section working group, extracted from [17] 



for example, found in Ref. [15]. As ithas been stated, within the SM the only missing 
aspect is now the Higgs boson mass. However, it has been further constrained at the 
LHC by the end of 2011. 

For a Higgs boson with a mass above 125 GeV, the search is dominated by bosonic 
modes, where Higgs decays to W + W~ and ZZ. For lower masses, the most important 
channels are bb, xx and yy [16]. This is explicitly displayed on Fig. 2.3 extracted 
from [17] that shows the Higgs branching ratio of each decay channel as a function 
of its expected mass. 

The largest decay mode of the SM Higss bosons is h — > W + W~ . If the Higgs boson 
would have a mass between 160 and 170 GeV, an excess would most probably (95 % 
of confidence level) have already been observed by the Tevatron, proton-anti-proton 
hadron collider situated in the United States of America. Within the W + W~ channel, 
purely leptonic (h — >• — ► Ivlv) and semi-leptonic (h — > W + W~ —> Ivqq) 

decay modes are being searched at the LHC. The decay of the Higgs boson to a pair 
of Z bosons has also been very attractive due to the production of a very clean signal. 
If the Z boson would decay as Z — >• 11, h — >• ZZ would produce 41 signature, which 
is very uncommon for the SM, leading to very small background. The production of 
h — > ZZ — > liqq and h -¥ ZZ -¥ llqq are also being analysed. At the low mass 
region, direct analyses of h — »■ yy, h — »■ bb and h — »• xx have been taken care of by 
the collaborations of the LHC experiments. Limits for the Higgs boson mass could 
therefore be set taking into account these analyses. 

In July 20 1 2 two collaborations of the LHC experiment have announced the discovery 
of a new boson state that, in principle, could be the Higgs boson [18, 19]. However, 
it is yet very early to state whether it is the so-expected Higgs particle, or a new 
bosonic particle with different properties. 

We are living therefore in a historical moment for particle physics. If the Higgs is 
not discovered by the LHC very soon, the SM will not have a mechanism to explain 
the generation of mass for fermions and bosons. Fortunately, it has been shown that 
the LHC has very strong indicators for a low mass Higgs boson of approximately 
126 GeV. For the moment, only time can provide us these answers. 
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In the previous sections, it has been shown why the SM of particle physics is one 
of the most successful theories in physics, which explains with great accuracy the 
experimental data we have so far. However, the SM suffers from few shortcomings, 
which led theorists to believe there could be new physics lying at the TeV scale. 
Here, we present some of these issues, which will be addressed in further chapters 
by the introduction of extensions of the SM. 



2.4.1 The Hierarchy Problem 

Many features of the SM arise from the BEH mechanism, the procedure to generate 
masses for the particles observed. However, only calculations at the lowest level of 
perturbation theory, also known as tree-level contributions, have yet been discussed. 
In reality, radiative corrections to the mass of a particle should be computed when 
one-loop corrections to the propagator are considered. 

Take for example a toy-theory with a single fermion t/r coupled to a massive scalar 
4> [20]. The lagrangian is given by 

£<j> = i fY^t + |3„ <t>\ 2 ~ m^\ct>\ 2 - Xfffcp. (2.29) 

As it has been done for the SM, we suppose the symmetry is spontaneously broken 
and the fermion obtains the mass m f = XfV /V2 at the lowest order. 
To compute corrections to the fermion mass within this toy-theory, it is necessary to 
consider one-loop contributions to the fermion propagator due to the scalar particle, 
as illustrated in Fig. 2.4a. The renormalized fermion mass is given by 

m f — M r + 8m f (2.30) 

where m J — X fv/*f2, and the correction to the mass depends on the introduction 
of a cutoff A to the theory: 

3X 2 f m f a 2 
5m/ = — rf-f-log — + ... (2.31) 



Fig. 2.4 a Example of a 
fermionic loop on a scalar 
particle, b illustration of a 
scalar loop on a fermionic 
particle 



(a) 



(b) 
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where ... indicates terms that vanish when A — >• oo or terms independent of the 
cutoff. Notice that the correction to the mass of the fermion depends explicitly on 
m f . In the limit where the mass of the fermions are very small, the lagrangian (2.29) is 
invariant under chiral transformations. Therefore, decreasing the mass of the fermion 
increases the symmetry of the theory. It is thus usually said that fermionic masses 
are protected by chiral symmetry. 

Similarly, we could compute corrections to the scalar propagator due to the fermion 
loop, as shown on Fig. 2.4b. When this is done, the mass of the scalar obtains a Sm s 
correction, which can be written as: 



Differently from corrections for fermionic masses, the mass correction for a scalar 
particle is found to be quadratically divergent. For this reason, nothing can actually 
protect the mass of the scalar particle of having very large corrections. If we assume 
the theory to be valid until the Planck scale, i.e. scale at which quantum corrections 
to gravity become important, corrections to the scalar particles will be of such order 
of magnitude. Unless the mass of the particle and its correction are at most of the 
same order, the theory is considered to have a naturalness problem. 
For the SM Higgs boson the same fact occurs. Obviously, corrections to the propaga- 
tor due to all types of particles on the loop must be considered. However, the largest 
contributions come from fermionic loops, where quadratic divergences appear. 
Since unitarity requires the Higgs boson to have a mass smaller than 1 TeV, quadratic 
divergences would indicate a large problem, considering corrections would have a 
much larger order of magnitude than the mass at lowest order. A counter term to 
cancel these quadratic divergences could solve the problem. However, such large 
cancellations would be highly un-natural. Conversely, one could adjust the cutoff 
down to 1 TeV. The result would be a theory that ceases to being valid at this energy 
scale. This is known as the hierarchy problem, and it is one of the reasons for which 
it is believed that new physics could indeed lie above the TeV scale. 



2.4.2 Fermionic Mass Hierarchy 

Another issue not well explained by the SM is related to the order of magnitude 
of the fermion masses. We have seen that fermionic lagrangian is defined by the 
Yukawa term, or Eq. 2.28, and they should acquire mass from spontaneous symmetry 
breaking. Concentrating on the quark terms, we have: 




A 2 - 6 ml log h 2m 2 f + ... 



(2.32) 




Ki Ql Hr - K Ql<Pu r + h.c, 



(2.33) 



which can be decomposed as 
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&yuk = K ~c'r ?r)Mv i^cpj + (4 S' R b' R )M D +h.C. 

where we have used the prime notation to represent the weak eigenstates. However, 
there is no reason, in principle, for the matrices M. u — -j= x',1 and M d — -j= X'J to 
be diagonal. In fact, the weak eigenstates are linear superposition of mass eigenstates 
and the transformation can be given as 



= V u c , Is' \= V d Is] (2.34) 



where V u and V d are unitary matrices that diagonalize the mass matrices Ai u and 
M d . 

In the SM lagrangian, charged currents (processes intermediated by a charged parti- 
cle, such as W ± ) are proportional to 

( u' L c' l i' L ) y» ( 4 \ = ( u L c L t L ) (Vj V d )y« (£\ (2.35) 

while neutral currents are proportional to 

( u' L c' L t' L ) I s[ \ = ( u L c L i L )(V}v u )y» (s[\. (2.36) 

Notice that Vu V u — 1, and there is no flavour mixing terms (at tree level) for neutral 
currents. However, flavour changing terms appear for charged currents coming from 
the CKM matrix defined to be Vqkm = V'uV d [10]. 

Experimentally, the CKM matrix is constrained and its parameters are defined accord- 
ing to the observed mass values for the fermions. It has been also realized that the 
CKM matrix follows a hierarchy and can be parametrized following the Wolfenstein 
parametrization [21]. 

In principle there is no natural explanation for the hierarchy observed in the CKM 
matrix. This can be related to the hierarchy observed in the experimental values 
obtained for the mass of the fermions. In particular, note that fermionic masses are 
generated from the same mechanism. It is therefore un-natural to obtain masses with 
so different order of magnitude. For instance the ratio wieiectron/wtop is of the order 
of 3 x 10" 6 . This issue is usually known as the fermionic mass hierarchy, and it 
could be one more indication of new physics at the TeV scale. 
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2.4.3 CP Violation 

CP symmetry is the combination of charge (C) conjugation and parity (P) symmetries. 
These symmetries state that the laws of physics must be the same if: 

• charge is conjugated, i.e. particles are replaced by anti-particles; 

• parity is exchanged, i.e. left-handed particles are replaced by right-handed particles. 

Both symmetries hold for electromagnetic interactions. It has been discovered how- 
ever that they can be separately violated for weak interactions [22, 23]. The CP 
symmetry has been proposed to solve the problem of individual C and P violations, 
with the statement that the CP should be the real symmetry between matter and 
anti-matter. 

Violation of CP has been experimentally discovered in 1964 by Cronin and Fitch 
[24, 25], through experiments made with K°. Later on, CP violation in B mesons 
has also been observed [26-28]. 

The main question regarding the subject nowadays is why is the strong force CP- 
invariant? The QCD lagrangian contain one term that could be responsible for a 
possible CP violation. This term is set to zero due to the non-existence of such events. 
The issue is also known as strong CP problem, and presents one more evidence for 
new physics at TeV scale. Besides, CP violation is often used to explain the unbalance 
of matter/anti-matter of the Universe. 

The most famous solution has been proposed by Peccei and Quinn [29, 30]. This 
theory postulates the existence of a new particle named axion, generated from the 
spontaneous symmetry breaking of a new symmetry. 



2.4.4 Dark Matter 

A different indication of new physics at the TeV scale is the evidence of dark matter 
existence. Dark matter is defined to be an undetermined type of matter that does not 
emit or reflect electromagnetic radiation. In fact, the subatomic composition of this 
new type of matter is not yet known, neither are its interactions. 
The first evidence was postulated by Jan Oort in 1932, in a study related to the search 
of orbital velocities of stars in the Milky Way. In 1933, when calculating orbital 
velocities of galaxies and clusters, Fritz Zwicky discovered that part of the mass 
expected by his calculations was missing, and denominated it dark matter [31]. 
Nowadays, several observations have confirmed dark matter evidence [32]. For 
instance, analyses combining galaxy cluster dynamics, supernova data and Cosmic 
Microwave Background (CMB) radiation with Big Bang nucleosyntesis are good 
indicators for its existence. Besides, precise measurements have been performed 
and recently confirmed the existence of dark matter at larger scales than the size 
of galaxies and clusters. It is believed that approximately 5 % of the known matter 
in the Universe is formed of atoms, while 23 % is dark. A relation of the matter 
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Fig. 2.5 WMAP data shows that the content of the present Universe is formed by approximately 
5 % of baryonic matter, i.e. atoms that compose planets and stars; 23 % of dark matter; and 72 % 
by dark energy, subject not discussed within this study. Dark energy should be different from dark 
matter, and responsible for the experimentally observed acceleration of the Universe. This image is 
based on the one produced by the NASA/WMAP Science Team [33] 



content of the Universe as of today, extracted from WMAP results, are displayed on 
Fig. 2.5 [33]. 

In principle, dark matter could be formed of either baryonic, non-baryonic content or 
yet by a combination of both. It could consist of planets, black-holes or even burnt- 
out stars in the galatic halo, for example. These are so-called MACHOS, i.e. massive 
compact halo objects, and are the most probable constituent of the dark matter if it 
is found to be formed of baryonic content [34]. 

Nevertheless, there are reasons to believe that dark matter could not consist of only 
baryonic content. Studies have shown that there is more dark matter in the Universe 
than the maximum number that can be calculated. The theory of Big Bang nucle- 
osyntesis and acoustic peaks in the CMB can be used to accurately determine the 
fraction of baryons in the Universe, and therefore to predict how much dark matter is 
expected to be baryonic. For this reason, it is widely believed that a significant part of 
the dark matter content is likely to be exotic, i.e. composed by non-baryonic content. 
From the SM of particle physics, only the neutrino would fill in all the requirements 
to be a non-baryonic dark matter candidate. Since it does not carry electric charge 
and it interacts very weakly to the particles of the SM, the cross section for neutrino 
production is very small. However, the mass density of neutrinos in the Universe 
has been calculated and is not enough to explain the matter fraction of the cosmic 
average density. Besides, there are studies of galaxies structure affirming that neu- 
trinos with a very small mass cannot be entirely responsible for the measured dark 
matter composition. As a solution, there could exist different types of neutrinos not 
yet observed, such as sterile neutrinos, to account for part of the dark matter content. 
In this case, clearly new physics is expected. 

There are usually three categories to separate possible non-baryonic dark matter 
content: (i) hot dark matter, which encompass theories with sterile neutrinos, (ii) 
warm dark matter and (iii) cold dark matter. These categories are related through the 
velocity of propagation of the particles, rather than their actual temperature. Particles 
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moving ultra relativistically fit in the first category, particles moving relativistically 
fit in the second category and finally, particles moving at classical velocities are 
categorized as cold dark matter. 

Different (and most probable) options are available. If dark matter is non-baryonic 
and not formed by neutrinos, it should be composed by new particle(s). The most com- 
mon suggestions are axions and WIMPs (or weakly interacting massive particles). 
The axion is the hypothetical particle postulated by Peccei and Quinn to solve the 
strong CP problem that appears in QCD [29, 30], as explained previously. They are 
supposed to be scalar particles generated from the spontaneous symmetry breaking 
of a new symmetry. 

WIMPs are massive particles, predicted by some theories Beyond the Standard Model 
(BSM), i.e. theory extensions to the SM. Examples of such theories are supersym- 
metry or theories with extra dimensions, both to be detailed on subsequent chapters. 
Candidates for WIMPs should interact only via the weak or gravitational forces, and 
not through the strong force. These are for example the neutralino or the gravitino 
in different variations of supersymmetric (SUSY) theories. SUSY theories not only 
predict a stable particle to be the dark matter candidate, but also solve several other 
issues presented here (for example the hierarchy problem). It is therefore one of the 
most studied BSM theories giving rise to different branches. More details on SUSY 
theories will be presented in Chap. 6. Finally, WIMPs candidates could also be for 
example new and unpredicted particles coupled to the SM only via gravitational 
interaction. 

All the options briefly discussed here to find the best candidate for dark matter largely 
indicate the existence of new constituents of fundamental particles. Moreover, an 
extended version of the SM would be necessary to incorporate such exotic particles. 
If this is true, chances are that these particles will be at the LHC energy range. 
Therefore, hints of new physics existence could appear very soon. 
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Chapter 3 

Gravitation in Beyond Standard Model Theories 



It has been shown that the SM is a successful theory that unifies three forces of 
Nature. However, general relativity is still left apart, and for many years several 
analyses seeking a general unification have been made. Throughout this chapter we 
shall briefly present the standard description used to quantize gravity, and assemble 
a quantum field theory that unifies all forces. The intention is also to display the 
limitation of the standard quantum gravity theory and introduce how new theories 
could well address SM shortcomings. 



3.1 Gravity as a Field Theory 

Initially, to introduce gravity to the theory, one could start from the renowned 
Einstein-Hilbert action for gravity [1]: 



J gra 







' R ' 


j d A x s 













(3.1) 



where k 2 = 167rGjv, g — det[g lw ], g^ u is the metric tensor and the curvature scalar 
R = g^ v R jw . The Ricci tensor is defined as 



where, 



r od = -7T ( d a9pa + d, 3 g ary - d a g al3 ) 



(3.2) 



(3.3) 



Varying S gra vity will result in the Einstein field equation [2, 3]: 

1 K 2 

Rll V _ _ R= T yv 

2 y 2 



(3.4) 
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where TP" is the matter energy-momentum tensor. 

For a quantum theory of gravity, it is necessary to consider quantum fluctuations for 
the gravitational field. One possibility generally chosen is: 

9 lily = ^liv ~t~ ft h\xv 

(3.5) 

gi Lv _ _ K hl iv + K 2 h » h Xu + __ 

where the flat space is represented by the metric r\ iLV — diag ( 1 , — 1 , — 1 , — 1 ) . Notice 
we work using natural units, i.e. h — c — 1. If one applies Eq. (3.5) in (3.1), the 
action can usually be expanded in orders of k and takes the form: 



where 



'-'grav 



Sqrav 



J d 4 x [C° grav + C\ rav + C 2 grav + . . .] (3.6) 



c 



Irav = \ 3* V d ° h " U ~\d a h d a h + d a h i) : h" : - i)J,„ : d 7,'"' 



+ n (^h 2 - \h^h^ + (hi^Kx - hh^ -RT, (3.7) 

where indices are raised/lowered with t]^ and in which 1Z — rj^TZ^, and 
Tinu = R^u with g fiu -* 

We are looking for a quantum theory to unify gravity with the other three forces of 
Nature, or in other terms, to incorporate gravity in the SM. For this reason, besides 
the graviton lagrangian we should consider all the other types of fields: 

S(cf>, ^ A„) = j d 4 x sf=~gC{4>, A Afj), (3.8) 

where the action S can be the action for either a scalar field cj), or a fermionic field 
ip, or yet a vectorial field A^. 

As it has been performed for the gravity action, it is necessary to consider quantum 
fluctuations in the action above. As a consequence, an interaction term between the 
gravity field and the SM field is obtained: 

S = J d A x(-^h llv T»\ (3.9) 

where 7^„ is the energy-momentum tensor corresponding to the matter content of 
the theory, i.e., 
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(3.10) 



The final action for a quantum theory of gravity, where SM fields can couple to the 
quantum field of gravity, should contain therefore a total lagrangian which is the sum 
of Eqs. (3.1) and (3.9). The inclusion of matter in combination with Eq. (3.4) results 
in the total lagrangian that can be re- written as: 



MG 



d a h d a h + d a h dnh 



ail 



(3.11) 



where the first term is also known as the Fierz-Pauli lagrangian [4] and we take 
11 = 0. 

Consider for example the action of a scalar field: 



Using the expansion of gP" in Eq. (3.12), we obtain: 



with 



7^ = d^A ~ (dpcfidP^ - m 2 4> 2 ) 



(3.12) 



(3.13) 



(3.14) 



From Eqs. (3.13) and (3.14) we can therefore construct Feynman diagrams related 
to the interaction between scalar fields and the graviton tensor represented by in 
these equations. 

We can similarly apply the same technique for fermions and vectors in order to obtain 
the interaction vertex between these particles and the graviton field. We thus consider 
the generic action that include vectors and fermions, 



d 4 x. 



1 



(3.15) 



where e is the vierbein, tp is the fermionic field, and is the electromagnetic tensor 
that contains the vector field An. In addition, we consider quantum fluctuations spec- 
ified by Eq. (3.5). Writing the lagrangian expression to be compatible with Eq. (3.9), 
the energy-momentum tensor for vector and fermionic fields are 
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(3.16) 



As a result, it is possible to obtain the total Feynman rules for the present theory. 
These rules describe the effective theory in which a graviton field couples to the 
usual vector, fermion and scalar fields. They are summarized in Fig. 3.1. 
We can see from Fig. 3.1 that if gravity couples to the SM fields, it is very weakly 
coupled. In fact, Eq. (3.9) reveals that the interaction is proportional to the square 
root of Newton's constant, or in other terms, to the inverse of the Planck mass. The 
interaction of the graviton with the SM fields will therefore be so small that it would 
be very hard to be tested by collider experiments. 

Nevertheless, the problem for quantizing gravity is more extensive. When comput- 
ing quantum corrections to the graviton propagator, we account for loop diagrams. 
To calculate quantum corrections at one-loop, we need to consider the lagrangian 
to quadratic order. When this is done, because gravity coupling is dimension-full, 
i.e. proportional to k, and taking into account the non-linear nature of the theory, 
divergences appear that cannot be re-absorbed by any fundamental parameters that 
we have introduced. Because of the dimension-full coupling and explicit powers of 
momentum in the self coupling of the graviton, divergences gets worse as more loops 
are considered. It is believed therefore that the physics we know ceases to being valid 
above the Planck scale, when quantum corrections become important. 
The dream of unifying gravity with all three forces of nature is not yet reached, 
however we have seen that progress has been made for the understanding of the 
underlying theory. Currently, one promising approach is to treat relativity as an 
effective field theory [1, 6]. Quoting J. Donoghue in Ref. [1]: "Since small quantum 
fluctuations at ordinary energies behave normally in perturbation theory, it is natural 
to try to separate these quantum corrections from the problematic (and most likely 
incorrect) high energy fluctuations. Effective field theory is the tool to accomplish 
this separation." 

Different extensions to the SM aiming to treat this problem using the effective field 
theory approach have been proposed. This is the case for example of the theories 
suggested by G. Dvali et al. [7-9] and by X. Calmet et al. [10-12]. Additionally, 
other BSM theories that intend to address different SM issues (such as the hierarchy 
problem or indirect signs of dark matter) may also achieve the unification of all the 
forces of Nature, including gravity. The following sections will present an overview 
of some particular SM extensions and describe their solution to the issues cited above. 
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+ C^px (k 3 - k 1 ) a + 

+ C^ax (k 2 - k 3 ) p + F^pvx] 




= if j [r c r dh Gp Vpa x5+ = -ig 2 {T",T b } 

fea6 recd/~i i read rebcs~t 

J J ^iivpXaS T J J ^pvpaS\\ 



Fig. 3.1 Feynman rules for a theory containing general spin-0, spin- 1/2, spin-1 and spin-2 fields, 
described by the lagrangian summarized in Eq. (3.9). The notation is coherent with [5], and the 
parameters Cp„ pa , Dp Vpa , and Gp, /piT xs are defined in Table 3.1 
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Table 3.1 Definition of the parameters C livpa , D^ p(T , E^ fXT , F fwpa \ and G llupa \s using the 
notation introduced in Ref. [5] 





— '/// o'/jyiT ~\~ iliifT' lun 'liiU'lnfTt 
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= V^k\ a k 2p - [n^ahvkip - r\ppk\ a k 2v + 'f] P ak\pk 2v + (n -*■ v)]; 


E pupa{k\ i ^2) 


= VpiAkipkia + k 2p k 2rT + k\pk 2a ) 
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= VppV<T\(k2u - h v ) + il lla Vp\(hv - k\ v ) 




+ i)p\ripo{k\ v - k 2v ) + (fi -* v)\ 


G pvpa\& 


= VnvVpaW ~ ViwVpSVaX + [VnpfluSVXa 




+ V f i\Vi'aVpS ~ VppVwVXS ~ VpiWvSVpa + (/•< ~* vj\. 



3.2 The Theory with a Scale-Dependent M p i 

The Planck scale M p i is generally treated as the fundamental scale of a theory. This 
is consequence to the fact that this is the scale where quantum gravitational effects 
become important. However, it could happen that these effects become large at a 
much smaller scale, such as the TeV scale. Consequences would thus be several. The 
trivial outcome would be that standard theories would not be any longer valid above 
this scale and extensions would have be considered. 

The primary extension we shall present is the consequence of a presumed renormal- 
ization of the Newton's constant, or in other terms, the Planck mass. The following 
will be based on Ref. [12]. 

It is known that the Newton' s constant is inversely proportional to the Planck mass: 



1 



K 

or, — = 

2 M p i 



7T 



1 



M 



(3.17) 



Pi 



where M p \ = M p //V87r ~ 2.4 x 10 8 GeV is the reduced Planck mass. Recall 
that the action of a general field coupled to gravity is given by Eq. (3.9). If we 
consider a generic massless scalar field, the resulting action shall be the combination 
of Eqs. (3.1) and (3.12): 



S — Sgrav ~r~ 



d^x ■ 



-9 



R 



1 



(3.18) 



Expanding g pl/ according to Eq. (3.5), the resultant action should be the sum of 
Eqs. (3.6) and (3.13). However, notice that Newton's constant is incorporated by 
definition in the expansion through the parameter k. If this dependence is isolated 
from hn V , the expansion becomes simply 



9fiv — Vpv ~*t~fo\xv- 



(3.19) 
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Fig. 3.2 Illustration of the graviton propagator and its one-loop correction due to the existence of 
scalar fields coupled to gravity in the theory 



The leading order term of the lagrangian displayed in Eq. (3.18) can be displayed as: 

R 1 h h 

-3 ~ -^hUh = (3.20) 

Hi Hi Hi Hi 

Since k was not absorbed in the definition of small fluctuations, loop correction to the 
propagator of the graviton could be re-interpreted as a renormalization to Newton's 
constant. 

In order to be more precise and consistent with the example we have been employing, 
consider our toy theory with a scalar field coupled to gravity. The graviton propagator 
must receive 1-loop corrections from the scalar particle in the loop, as illustrated in 
Fig. 3.2. Neglecting the index structure, the graviton propagator including the 1-loop 
correction, can be displayed as 



Pqrav (<? ) 



q 2 k 2 q 2 k 2 q 2 k 2 



(3.21) 



However, similarly to what has been done to calculate the Higgs boson mass 
corrections, it is found that the integral included in E is quadratically divergent. 
The authors of [12] proposed to absorb this part into the redefinition of Newton's 
constant, obtaining: 



1 1 c 



A 2 , (3.22) 



G r ™ G h N are 16n 2 



where G' ™ is the renormalized constant measured by experiments, c is a new para- 
meter, and A is the ultraviolet cutoff. 

If the Planck mass is scale-dependent, then the true scale at which quantum effects 
become important is /i* : 

Gtv(m*)~m; 2 - (3.23) 

Requesting G r ™ to be of the order of the Planck scale, G^ re (A) ~ A and the 
ultraviolet cutoff of the order of (i*, results on c ~ 0(1O 32 ). 

If we now assume the possibility that spin- 1/2 and spin-1 particles also couple to 
gravity 1 [13, 14], and compute the total number of particles required to generate the 



Consider the simplest lagrangian to describe such theory, where the coupling with fermions is 
treated with the vierbein formalism and the vectorial coupling employs the traditional field strength 
tensor. 
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Planck scale from the fundamental scale jtx*, we obtain: 

, 2 = — ^ (3 24) 

where Ni — + Nj, — 4-Ny, is the number of scalar fields in the theory, N$ is 
the number of fermion fields and Ny vectorial fields. 

Equation (3.24) indicates that roughly 10 32 light particles coupled to gravity (within 
a new hidden sector), are necessary to address the hierarchy problem and generate 
the Planck scale from the fundamental scale, here assumed to be the electroweak 
scale. In principle, this observation may appear to be excessively strong. However, 
more generally accepted theories, such as the theory with Large Extra Dimensions 
[15-17] also predict the existence of a large number of new particles, as we shall 
confirm in Sect. 3.3. 



The Effective Lagrangian 

A main consequence of the theory presented above is that there must exist a large 
number of particles coupled to gravity through a hidden sector. Because the theory 
contains a massless graviton, from the present moment on, we shall refer to it as 
massless graviton model, or MGM for short. To obtain the effective lagrangian of 
the theory, one must proceed by defining the action of each field and couple it to the 
gravity by considering quantum fluctuations for the metric. Equation (3.8) exhibits 
explicitly the procedure. 

Feynman rules have been calculated for a general theory with spin-2 particles cou- 
pled to SM fields. They are displayed in Fig. 3.1. For the MGM theory, each SM 
particle that couples to gravity within the hidden sector shall follow these rules. If 
one integrates over the total number of new particles, the effective scale of the theory 
can reach the electroweak scale through the scale dependence of Newton's constant. 
Therefore, Eq. (3.23) can be translated to the Planck scale, 

M p l(fj,*) ~ M*, with M p i = M p ,/V^. (3.25) 

As a consequence, the MGM theory will have the same Feynman rules as displayed 
in Fig. 3.1, with (3.25) included in the k parameter. Because M p /(/i t ) ~ /i*, the 
gravity coupling to the SM fields becomes important at the electroweak scale once 
[1+ is chosen to be O(TeV). 

We therefore have a complete theory that incorporates the gravity force, and in which 
all types of fields can only couple to gravity through a hidden sector. Now, graviton- 
related processes can have reasonably large cross sections at colliders, and they need 
to be taken into account. Cross-sections and decay-rates based on this particular 
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theory can be calculated in order to impose limits for its parameters taking into 
consideration electroweak constrains. 

It is important however to stress that this is an effective theory, and we have calculated 
vertices up the first order of h^ v . It signifies that there is at most one graviton field at 
each vertex. We can start by computing cross-section at colliders for processes with 
one graviton emission, as for example, proton + proton — > gluon + graviton. 
This calculation shall be performed in Sect. 3.5, followed by full phenomenological 
analyses at the LHC. In addition, we shall compare these cross sections calculated for 
the MGM theory against cross sections computed based on different BSM theories. 
With this aim, an introduction for different BSM theories that contain spin-2 particles 
are presented below. 



3.3 Extra Dimensions 

Theories with extra spatial dimensions to work as an extension to the SM theory, have 
become extremely popular in the last decade. However, these ideas are not new. They 
have started to appear in the beginning of the twentieth century just after Einstein's 
theory of relativity had been formulated. 

The first idea was proposed by Theodor Kaluza, a german physicist, who studied the 
consequences of the theory of relativity of Einstein assuming a 5 -dimensional space- 
time metric, instead the usual four. For his own surprise, Kaluza discovered that it was 
possible to unify Einstein's theory of gravitation with Maxwell's electromagnetism 
by the addition of one dimension [18]. A few years later, Oskar Klein reformulated 
the theory proposed by Kaluza and demonstrated that very small and compact spatial 
extra dimensions could exist. However, these theories predict that gravity should be 
very weakly coupled to SM fields, and for this reason they were, for a very long time, 
left aside. 

In the search for unified theories, in the earlies 80's, these ideas started to re-born 
and gave rise to different types of new theories, from strings to extra-dimensional 
theories, the latter being the subject we will now focus on. 

Theories with extra dimensions start from the idea that we could live in a world 
described by more than 4 dimensions. However, different from string theories, where 
in total 10 or 1 1 dimensions are required, the assumption is that only a small number 
of extra-dimensions are added to the usual four. The idea is to look for an extension 
to the SM by adding very small and compact extra dimensions. 
What would be the 4-dimensional consequences of having a Z)-dimensional space- 
time? This was one of the questions asked by N. Arkani-Hamed, S. Dimopoulos and 
G. Dvali, authors of the first proposition to employ Kaluza and Klein calculations 
for the search of a SM extension. Nowadays, this theory is also known as ADD 
theory [15-17], or theory with large extra dimensions, name originated from the 
prediction related to the size of the radius of the extra dimension. Several branches 
of extra dimensional theories have been developed ever since. We can cite for instance 
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the theory with warped extra dimensions, proposed by L. Randall and R. Sundrum 
[19, 20]. Both SM extensions are discussed in the following sections. 



3.3.1 Large Extra Dimensions 

Arkani-Hamed, Dimopoulos and Dvali have proposed to use the framework of extra 
dimensions in principle to address the hierarchy problem. It is known that there are 
at least two significant energy scales in nature: the Planck scale, where quantum 
corrections related to gravity become relevant; and the electroweak scale, where 
masses of the SM particles lie. One of the key motivations for creating extensions 
to the SM theory is, as explained previously, the large value for the ratio M p i /M ew . 
Although the electroweak scale has been well tested by recent experiments, it is not 
easy to explore the gravitational sector due to the large order of magnitude of the 
M p ] scale. 

For this reason, the fundamental scale of the SM theory has been assumed to be 
M p i . However, similarly to what has occurred to the MGM theory, a reinterpretation 
of some parameters could bring the value for the fundamental scale down to the 
electroweak scale, and the naturalness of the theory would be recovered. 
Consider a D-dimensional theory, where D — 4 + 5 and 8 is the number of extra 
dimensions. Einstein's action for gravity is now, 



where x ft are the usual 4-dimensional coordinates, y are the extra-dimensional coor- 
dinates and g is the higher dimensional metric corresponding to g. If one impose 
the so-called KK reduction, i.e. if we integrate the action over the extra dimensional 
coordinate and impose that the 4-dimensional action should be the same as Eq. (3.1), 
we obtain a relation between n and nrj. This relation is clearer when written in terms 
of the Planck constant: 



where V& is the volume formed by the extra dimensional space, proportional to the 
radius of the extra dimension(s). 

The fundamental scale of the theory is now Mo and can be of the order of the 
electroweak scale, according to Arkani-Hamed, Dimopoulos and Dvali's philosophy. 
The Planck scale is consequently generated by the fundamental scale through the 
volume of the un-observed extra dimensions. 

Values of the radius of the extra dimensions R are related to the Planck mass according 
to Eq. (3.27). If the fundamental scale is chosen to be of the order of 1 TeV, a theory 
with one extra dimension, i.e. 5 — 1, is excluded, since R ~ 10 11 m. For 5 = 2, 




(3.26) 




(3.27) 
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R ~ 0.2 mm, which is within experimental limits (low gravity experiments) [21]. 
For this reason, from now on we consider ADD theories with 6 > 2. 



A Realistic Model 

Now, it is clear that it would be possible to address the hierarchy problem within 
the framework of extra dimensional theories. We can therefore proceed to present a 
realistic theory where gravity couples to SM fields. 

It is usual to name the full dimensional space as the bulk, and the 4-dimensional 
space as the brane. For this reason, theories with extra dimensions are also known 
as brane-world models. We shall also make use of these terms. 
The first concern should be the accessibility of the fields to the extra dimensional 
coordinate. If a general field <j> is free to propagate through the bulk, the consistent 
action should be [21] 

Sbulk — J d 4 xd 5 y^/^£<p(x,y). (3.28) 
Rather, if the field must be constrained on the 4-dimensional brane, 

Sbrane = J d A x d 8 y ^ U to 5 s (y - y ). (3.29) 

To find out the consequences of a general field propagating in the bulk, consider 
once again a toy-theory described by a generic scalar field coupled to gravity. Now 
the space-time is 5-dimensional and the scalar field is free to propagate through the 
bulk. 

The relevant action for such theory is 

-.■¥]. (330, 

being y the coordinate of the extra dimension, and the upper case indices running 

through all the dimensions, i.e. M, N = 0, 1, 2, 3, 5. 

The compactness of the manifold is revealed by the periodicity of the field 

cj)( x f, y ) = (f,( x », y + ItiR), (3.31) 

and a Fourier expansion can be employed: 

<M*'\ y) = i — - yoto + y] \v>n cos (^) + <p„ sin (^)l • (3.32) 

V2tt/? ^ VttR l k R j \RJl 
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(a) 



— . — — x 




Fig. 3.3 3-dimensional spheres represented by the infinite sum of the 2-dimensional circles along 
the y-axis. The left sphere has been constructed using the software Mathematica® 



To apply the KK reduction, we must integrate the 5 -dimensional action and relate 
it to the known 4-dimension action for a scalar field. Introducing Eq. (3.32) in the 
5-dimensional action and integrating it over y, we obtain 



2 



d fl ip„ + mltplj 

.n=0 n=\ 

(3.33) 

where m„ — m 2 + . Equation (3. 12) is the corresponding 4-dimensional action to 
(3.33). Comparing both equations, we can deduce that the higher dimensional scalar 
field is observed in the effective theory as an infinite tower of fields with masses m n . 
These excited modes will be named KK excitations, referring to the calculation first 
performed by Kaluza and Klein. 

The consequence of a field free to propagate in the bulk is thus the appearance of a 
tower of 4-dimensional fields with the same quantum numbers but increasing masses. 
In order to visualize how such consequence would be possible, we shall make an 
analogy with a lower-dimensional theory. Consider the three-dimensional sphere 
illustrated in Fig. 3.3a. How could one describe it in a two-dimensional theory? 
The sphere could, for instance, be characterized as a sum of an infinite number of 
2-dimensional circles, located in the vicinity of each other, as displayed in Fig. 3.3b. 
These circles should have different radii, forming the sphere when placed alongside. 
The KK-excited modes of the scalar field could be interpreted following the same con- 
cept. In this sense, the zero-mode field could be identified with the SM field. The sub- 
sequent modes should therefore exist with higher masses. Because the D -dimensional 
space-time metric considered is flat, the expansion of the field has been performed 
in terms of trigonometric functions. Therefore, the mass separation between two 
successive KK modes is 1 /R for the ADD theory. 
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The concept of having a D dimensional object projected into the D — 1 dimension 
is also known as compactification. It can be applied for any type of field, (i.e. for 
spin-0, spin- 1/2, spin-1 and spin-2 fields) provided that it propagates through the 
bulk. The consequence will consistently be the same: the reduced theory will have a 
similar spectrum with larger degeneracy. 



The Effective Theory 

In the ADD theory, only the graviton field is allowed to propagate through the bulk 
while the SM fields are kept in the 4-dimensional brane. The theory has been exten- 
sively discussed by several authors [5, 22, 23]. Here we summarize how to obtain 
the 4-dimensional effective theory starting from the D-dimensional action. 
Equation (3.28) describes the action of a particle propagating in the bulk. Therefore, 
in the ADD theory, this is the action of the graviton. Quantum fluctuations for the 
metric must be adapted to the higher dimensional space-time: 

9MN — VMN + HD h[4N- (3.34) 

The lagrangian of a particle in the brane is given by Eq. (3.29), where now <f> 
should refer to the SM fields. The energy momentum tensor has also to be defined in 
£)-dimensions, and when the gravity field is weak it becomes: 

Tmn = T M1/ jf M rf N 8 s (y). (3.35) 

In 4-dimensions, Eq. (3.9) determined the effective coupling of the graviton with a 
SM field. In £>-dimensions, the corresponding lagrangian is given by: 

Hnt = ~Y h MN T (3.36) 

In order to start the dimensional reduction, assume the extra dimension is compact- 
ified in a torus. By doing so, the higher dimensional volume can be approximated 
as Vg ~ (2nR) s , and the graviton field, which is periodic h^ tu (y) — h lw (y + 2nR), 
allows for the expansion 

u CvM „\ MN \ > Any R n\ 
hMN(x,y)= > 7=^ e ■ (3.37) 



Applying the expansion (3.37) in the total action, and performing the KK reduction, 
the part of the lagrangian proportional to the graviton field should be 

= J2>$ r ' 8 '' (3-38) 
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in addition to an extra term related to the dilaton field, which we will not focus on 
within the present study. 

The zero mode of the graviton field should be massless and coupled to the SM 
particles proportionally to the inverse of the Planck mass, giving rise to the weak 
nature of gravity as we know it. The following KK modes should be massive, and 
the mass splitting between consecutive modes 2 should be Am — l/R. 
Notice the similarity between the lagrangian obtained for the 4-dimensional massless 
graviton (MGM) model presented in Sect. 3.2 and the effective lagrangian of a higher 
dimensional graviton described by Eq. (3.38). For a generic coupling constant, both 
lagrangians are equivalent. The Feynman rules for the ADD theory can thus be 
identified with the rules specified in Fig. 3.1. The difference lies in two factors: i) 
for the MGM theory only one massless graviton is predicted, while for the ADD 
theory also a tower of massive gravitons is expected; ii) for the MGM theory, the 
coupling constant is scale dependent and can be of the order of the electroweak scale 
because of the existence of a large number of invisible particles that can only feel the 
gravitational force. Conversely, for the ADD theory, the effective scale can reach the 
TeV scale due to the existence of extra dimensions, which results in the prediction 
of the KK tower of the graviton. 

As it has already been stated, the mass value for each KK excited mode is proportional 
to \n\/R. Therefore, the mass splitting between two successive modes is so small that 
it cannot be measured. It would require a large number of states in order for the mass 
splitting to be at experimental range. In other words, a large number of KK modes 
of the graviton would be necessary to compensate the factor 1 / M p i coming from 
each interaction vertex. Therefore, for cross section and decay rates calculations, it 
is more relevant to obtain the contribution of the sum over all the KK states. One 
can thus integrate the number of KK-modes between | n | and | n \ + dn [23] . The total 
number will be given by 



dN = 



2tt^ 2 

r(<5/2) 



\n\ s ~ l dn = 



F(<5/2) 



M 



pi 



M 2 + 5 8^ 



-dm 



(3.39) 



using that m — \n\/R and Eq. (3.27) to relate the higher dimensional radius R with 
M p i and Mq. As consequence, the cross section calculation for graviton production 
will have an extra integration over the KK mass states. 



d L a 
dt dm 



2tt^ 2 ' 
F05/2) 



da„ 



Mjf s 8tt dt 



(3.40) 



where da m /dt is the differential cross section for producing one KK graviton of 
mass m. 

In the next chapter we shall compute cross sections for producing KK graviton(s) at 
the LHC, and compare them with results obtained by taking into account different 



2 Here we assume a flat higher dimensional space-time metric. 
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BSM models such as the MGM or the RS model, which will be presented in the 
succeeding section. We can, nevertheless, already state that the dependence of the 
cross section on M p i can be replaced by Mp when the emission of the full tower of 
graviton is being considered. Similarly to the MGM theory, the ADD theory predicts 
therefore graviton emission processes at the LHC energy reach. 



3.3.2 Warped Extra Dimensions 

Several branches of extra-dimensional theories have been suggested since ADD was 
first proposed. Various scenarios can be created by considering different particles 
in the bulk, or different higher dimensional space-time metric. Within this line of 
thought, L. Randall and R. Sundrum have proposed an alternative to the ADD scenario 
to address the hierarchy between the electroweak and Planck scales [19, 20]. They 
suggested a theory with 5 -dimensional space-time that could be described by the 
total metric 

ds 2 = g M N dx M dx N ', (3.41) 

where again upper case indices represent the full set of coordinates or M, N — 
0, 1, 2, 3, 5, and 



9MN 



Equation (3.41) can also be written as 

ds 2 = e~ 2k | v| n^ v dx 1 ' dx v - dy 2 , (3.43) 

where x'' are the usual 4-dimensional coordinates, k is of the order of the Planck 
scale and y is the fifth-dimensional coordinate. 

We shall start as customary, by writing the 5-dimensional action for gravity in the 
bulk: 

f 4 r^T R 

Sgravity — / d xdy^-g 

J K 5 



(3.44) 



where K5 — V I671-/M5 is the higher dimensional coupling and g = det[gMN]- 
Employing Eq. (3.42) in g, we obtain the relation ^/—g — e _4 *' v ' . 
For the authors, although the 4-dimensional space-time is flat, the exponential factor 
works as a warping factor for the total metric. The symmetry S l /Z2 is applied, i.e. the 
periodicity of the fifth dimension is used to identify both y y + 2ttR and y — v. 
Both extremity points of the orbifold, y — and y — tt R, do not transform under 
the Z2 symmetry [24]. 
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In principle, only gravity will be allowed to propagate through the bulk while other 
fields must be localized at one of the branes. We can thus choose y — as the initial 
point and define the action at each 4-dimension brane, located at the orbifold fixed 
points v = and y — tt R: 



Svi.s — / d X^J Qvis £vis 



Shid = J d 4 Xy/-g hi d Chid g h Jv = Qiivix^, 0) 



(3.45) 



-2kirR , 



From these equations, we obtain g^f — rj^ v , whereas g^'J — e~ i}^,. 
In order to understand the consequences of this exponential factor, let us suppose the 
Higgs field is locked in the visible brane, i.e. localized at y = irR. The Higgs action 
is: 



rihiggs 



= / d\S- 



9vi 



qZd^d v h 



\\\H\ 



(3.46) 



Applying gJ'J — e zk7TK g^, and re-defining the Higgs field as H — >■ e k,KK H, we 
obtain 



S'^ s = J d 4 x \rTD^D v H - A (\H\ 2 - e~^ R v^f 



(3.47) 



Notice the exponential factor redefines the symmetry-breaking scale, which is now 

(3.48) 



v = e~ M v . 



We can therefore conclude that the theory proposed by Randall and Sundrum predicts 
that any mass parameter at the visible brane is exponentially reduced from its absolute 
value at the hidden brane. If we assume the natural scale of the theory to be M p i at 
the hidden brane, we can require the visible brane to have masses of the order of the 
electroweak scale. This requirement constrains the product kR ~ 12, and address to 
the hierarchy problem in a very natural manner. 

In order to obtain the relation between the Planck scale M p i and the effective scale 
Ms, one should perform the dimensional reduction in the action (3.44). Identifying 
the resulting action with the 4-dimension action, we obtain: 



M 



Mi 



1 



-2knR 



(3.49) 



From (3.49) it is clear that M p i depends strongly on the limit kR. It results that 
the huge gap between the Planck scale and the electroweak scale exists due to the 
warping factor of the 5-dimensional metric. However, different from the ADD theory, 
the fundamental scale of the RS theory is still the Planck scale, and is reduced to the 
TeV scale at the brane localized at y = irR. 
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A Realistic Model 

In order to know how the exponential factor of the metric will influence the dimen- 
sional reduction, consider a 5-dimensional scalar field <f> allowed to propagate in the 
bulk [25-27]. The relevant action is again 

Sf = J d\dy^g l -[g MN dM<j>d N <P-m 2 <p 2 ~\ , (3.50) 

which can be written, considering Eq. (3.43), as: 

Sf = J d\ dy e~ 4kM l - [e 2kM d^d^cp + 85(f) d 5 - m 2 2 ] . (3.51) 

Using the periodicity of the field with the symmetry imposed, we can expand the 
scalar field as 

y) = -7^T/(^) My). (3.52) 



If fn(y) are chosen to satisfy both: 

pirR 

2nR 



l — P dye- 2kR My)f m (y) = 6 nm , and, (3.53) 

?R J-wR 

[-e 2k ^d 5 (e- 4k ^d 5 ) + m 2 e~ 2k ^] f n (y) = m 2 n f n (y), (3.54) 

Equation (3.51) becomes 

S l d = \lLj d * x [v^d^ndutpn ~ mW] ■ (3.55) 

Notice that the zero-mode, or n — 0, represents the action of a scalar field in the 
usual 4-dimensional theory. Equation (3.55) is therefore the dimensionally reduced 
action. 

The solution of Eq. (3.54) is 



My) 



e 2kM r v \ /w B c* ,yl M 



where J a and I a are Bessel functions of order a, a = ^J4 + 'j^, m n are the masses 
of each KK state and the coefficients N n and (3{m n ) are determined by normalization 
and the boundary conditions in the branes: 

4>ix\ y) |y=o,^= 85^, y) \ y=0 ,*R= 0. (3.57) 
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Therefore, 

2J a (f)+ J' a (f) 

0( - m " ) = - 2 Y a (^) + fY' a cfy (3 - 58) 

Due to the warping factor of the metric, the KK mass states m n do not increase 
linearly as happens for the ADD theory. They can be calculated and given by [27] 

m„ ~ (n + | - ^) Trke-* kR , (3.59) 

for m n k, kR ^ 1 and n = 1, 2, 3, .... 

Once more, the consequence of a 5-dimensional field propagating in the bulk is the 
appearance of the KK tower. However, for RS theories, two successive KK mass 
states are not separated linearly, but their separation depends on the exponential 
factor of the metric. 



The Effective Theory 

In order to obtain the effective theory, let us assume that only gravity can propagate 
in the bulk. The SM fields must be constrained in the 4-dimensional brane localized 
at y — 7r R. The action of gravity and SM fields are given by Eqs. (3.44) and (3.45). 
The graviton field must be expanded as it has been performed for cj) in Eq. (3.52), 
and fluctuations of the metric must be taken into account: 

= e~ 2kM rip, + V- ( 3 - 6 °) 

Imposing (3.60) in (3.44) and (3.45) we obtain the interacting lagrangian of the 
graviton and the SM fields: 

C=\hf,T^-e- k ^\h^\ (3.61) 

The massless graviton, localized at the hidden brane, couples to SM fields propor- 
tionally to the Planck mass, whereas the first KK mode can couple to the SM fields 
proportionally to the TeV scale, accounting for their localization at the visible brane. 
The zero mode of the graviton can therefore be identified as the massless graviton 
predicted by usual quantum gravity, whereas the first mode can have a mass of the 
order of the TeV scale. The mass of the n-th mode of the graviton is given by [28]: 

m n =kx n e- krR (3.62) 

where x n are the roots of the Bessel function of order 1 (for example x\ — 3.8317). 
Once more, due to the similarity of Eqs. (3.61) and (3.9), the Feynman rules for KK 
modes of the graviton coupled to the SM fields are exhibited in Fig .3.1. However, the 
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generic coupling displayed in Fig. 3.1 is now given by e 



-IcttR I M 



/M p i for the coupling 



of the first mode of the graviton with the SM fields. Notice that because of the 
exponential factor present in this coupling, RS theories are considered to be much 
more strongly interacting than ADD theories with large extra dimensions. 



3.4 The Graviton Propagator 

We have presented so far a complete overview of BSM theories featuring spin-2 
particles, from 4-dimensional theories to extra dimensional theories. Further ahead, 
we shall display the phenomenology of these theories at hadron colliders. We are 
lacking therefore only an analytical discussion on the graviton propagator, which is 
one of the factors differentiating massless graviton theories to the ones containing 
massive spin-2 particles. We shall follow the formalism presented in [29] to derive the 
spin-2 propagators. Nevertheless, the reader could check [30, 31] for the derivation 
of the propagator in full details. 

To obtain the graviton propagator, we start from the Fierz-Pauli lagrangian [4] for 
massive spin-2 particles: 



The last term is related to the mass of spin-2 particles and the relative — 1 coefficient 
between h^ v and h 2 is an assumption named Fierz-Pauli tuning [29]. There is 
no known symmetry that imposes such choice. However, the Fierz-Pauli tuning is 
usually employed to prevent the theory from containing scalar ghosts. 
If (3.63) is integrated by parts, the action/lagrangian can be re- written as 



S, 



'grav — 




(3.63) 




(3.64) 



where O^' 



is given by 




and satisfies the symmetries 




(3.65) 
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To extract the graviton propagator, we should apply the quantum prescription and 
move to the momentum space. The propagator will be the operator V a p,aX which 
satisfies the relation 



(3.66) 



For Eq. (3.65) that is 



p 2 — m 2 



2 ( V* " 


PaP/3\ ( 


PvPQ\ 
m 2 / 




m 2 )V va 


PvPa 
m 2 


"3 (Vp> - 


PnPv\ ( 
2 ) W 


PaP/3 
m 2 



(3.67) 



To derive the massless propagator we follow a similar prescription. We start from 
Eq. (3.63) in the limit where m — >• 0. We can again integrate it by parts obtaining: 



S - f /r - h P» v > a P h r 

^grav — / 2 ^ rl otp 



(3.68) 



m-i-0' 



where S^ap — O^^p 
However, now it is necessary to fix the gauge in order to obtain the propagator. The 
most common choice is the de Donder gauge: 



" \ d « h = o, 



which is achieved by adding the gauge fixing term for the lagrangian: 
C gf = (d a h% - \d,h} gT (d a h° v - \d v h\ . 



(3.69) 



(3.70) 



Notice we use h a a = h. 

The total lagrangian is given by the sum of the gravitation and the gauge fixed 
lagrangians: 

C = C grav + C g f = \ V h a p. (3.71) 



with 



(3.72) 



To obtain the propagator, we once more go to the momentum space and apply the 
quantum prescription. We obtain simply: 
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(3.73) 



When the massless (3.73) and massive (3.67) propagators are compared, one can 
realize a difference in the last term, which appears even in the limit where m — > 0: 



m-»0 q 



rf a r) vti + rfPr) va 



(3.74) 



This is the famous van Dam-Veltman discontinuity [32, 33]. 

To understand how such discontinuity arises, let us consider the study made by Dicus 
et al. in Ref. [34]. They have performed a full analysis on the graviton propagator 
both for the massive and the massless cases. 

First, consider an off-shell graviton particle. If the graviton is massless, it has the 
usual five polarization tensors (J — 2, J z = 0, ±1, ±2) in addition to the (J — 0) 
component: 



piv,,ap „ ^ e/ ll/ (2, A) e a!3 *(2, A) — e^(0, 0) e a3 *(0, 0). (3.75) 
A=±2,±l,0 

However, an on-shell massless graviton has only transverse components: 

pwaf}„ ^ e ^(2,A)e ai3 *(2, A). (3.76) 

A=±2 

The authors of Ref. [34] affirm that besides the cancellation of (/ — 2, J z — 
±1) components, the (/ — 2, J z — 0) is cancelled by (J — 0), leaving only 
(/ — 2, J z — ±2) components. 

However, a massive graviton does not have the (/ = 0) component: 

V 1 ' 1 '-" : - ^ e^(2, A)e Q3 *(2, A). (3.77) 

A=±2,±l,0 

Thus, when the zero-mass limit of Eq. (3.77) is considered, the (/ — 2, J z — ±1) 
and ( J — 2, J z — 0) components are cancelled, but the (/ = 0) component remains. 
Therefore, the discontinuity arises from the fact that a virtual massless graviton 
has a component with zero angular momentum (/ = 0), that serves to cancel the 
(/ — 2, J z — 0) component when the graviton is on-shell, whereas the massive 
graviton does not have the (J — 0) component at all. 

It is clear now that theories containing massless gravitons should be treated separately 
from theories comprising massive such particles. In this context, the similarity of the 
Feynman rules between the MGM theory and extra dimensional theories disappear 
at the propagator level. One could therefore investigate if the discontinuity of the 
graviton propagator can also been seen at the cross section level. 
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With this question in mind, we have analytically computed cross sections of graviton 
emission processes for both cases: emission of a massless graviton and emission of 
a massive graviton, taking the limit where m — s- 0. This will be the subject of the 
next section. 



3.5 Graviton Emission Cross Sections 

On Sect. 3.1, we have shown the Feynman rules for a general theory in which SM 
fields are coupled to the graviton field. These rules are displayed in Fig. 3.1, con- 
sidering a generic coupling presented through the parameter k. We shall now use 
these rules to analytically calculate cross sections for graviton emission processes in 
addition to quarks and gluons, and we shall follow Ref. [35]. 

Quarks and gluons can be produced in combination with a graviton G from partonic 
sub-processes q + q — s- G + g,q+g q + G,q+g q + G,andg + g — > g + G. 
Primarily, we have calculated the leading order contributions at the parton level for 
a theory featuring a massless graviton. Regarding the notation, we use G for the 
massless, and non-tilded G for the massive spin-2 particles. For simplicity, at this 
stage, we treat all quarks as being massless. The sum over the polarizations of the 
graviton is given by Eq. (3.76). We have found that the polarization and color averaged 
cross section for q + q — >• g + G is given by 

da(q q -*■ gG) n 2 a s / s 2 + 2t 2 + 2 1 s 
dco&9 = 72 V ^ 

where a s — g 2 /4w, g s is the strong coupling constant, k is given on Eq. (3.17), s and 
t are the Mandelstam variables: t — — l/2s(l — cos 9), and 9 is the angle between 
the beam and the visible final particle. The cross sections for q + g — >• q + G and 
q + g — s- q + G are given by 

da(qg^qG)_ k 2 a s /2s 2 + 2s t + t 2 
d cos 9 ~ 192~ V s~t 

One can use crossing symmetry from Eq. (3.78) to obtain the corresponding matrix 
element of Eq. (3.79). These cross sections can be deduced from [36]. Finally we 
also obtain the cross section for g + g —> g+G which is given by 

da(gg^gG) = 3 k 2 a s / s 2 + st + f 2 \ 2 
dcos9 128 \sjt(s + t)) 

Notice that these calculations can be substantially simplified considering that the 
graviton is on-shell and therefore any contraction of the style T^ v ka vanishes, where 
T^ v is the energy momentum tensor and kn is the momentum of the graviton. 



(3.78) 



(3.79) 
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Subsequently, we have computed the differential cross section contributions for a 
theory featuring a massive graviton. Now, the sum over the polarizations of the 
graviton is given by Eq. (3.77), the massive graviton propagator. We have obtained 
the following respective results. The cross section for q + q — > g + G is given by: 



dajqq ->;G) 
dcosQ 



9 

K a s 



144 



— 

1 -m 2 /s\\ 



f (s-m 2 ) 2 

2—- -5+4 

v Aut 



Aut \ / ,m 2 ,A 
(s - m 2 ) 2 ) V + ^T^ ) 
4„, + 



(s — m 1 ) 



2\2 



( u — t \ / m% Y 

\s — m 2 ) s 



(3.81) 



where s, t, u are the Mandelstam variables, and t, u — — l/2.s(l — m 2 /s){ \ =pcos 8). 
The cross section for q + g — > q + G can be obtained from this expression by 
crossing s t : 



da (q g — >• q G) k a s (— t/s)(l — m js) 



d cos 6 



384 



(1 -m 2 /t) 2 



(Aus \ ( ,m 2 ^A\ 
j -<fs^)( 1+( t)) 



(it- m 2 ) 2 

+ 2-- -5 + 4 

\ 4ms 

(\ 2 2 
2 — 
t - m 2 ) v t ' 



(t - m 2 ) 2 
Aus 



(t - m 2 ) 2 



)t(-(t> 2 ) 



(3.82) 



As in the massless case, the cross section for q + g — > q + G is also the same as that 
ofq+g—^q + G. For the process g + g — >• g + G, we find 



dcrjgg ->■ gG) 
dcosO 



2 



256 (1 - m 2 /s)(\ - cos 2 6>) 

"(3 + cos 2 0) 2 (l + (^) 4 ) 

„,2 



■4(7 + cos 4 6)- 



(> + <T>') 

■6(9-2cos 2 6» + cos 4 6»)(^-) 2 j . 



(3.83) 



If we take the limit where m — >■ of Eqs. (3.81), (3.82) and (3.83), we obtain the 
same results as Eqs. (3.78), (3.79) and (3.80) respectively. The discontinuity existent 
for the propagator of a spin-2 particle vanishes at cross section level when processes 
with an on-shell graviton are considered. 

The smoothness of this limit has been pointed out in several references, including 
[22]. Particularly, Mirabelli et al. have pointed out that due to helicity conservation, 
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only the polarization modes that correspond to the massless graviton are produced 
in the reaction being analysed. However, in order to understand theoretically the 
smoothness of this limit, we have carried out a small calculation that verifies it for 
the qq — > gG process. The validation follows. 
From Eqs. (3.73) and (3.67) it is possible to obtain the relation 



lim V M u,a/3 — '/ , ., ; ,., ! j 



1 



(3.84) 



On the calculation of the cross section of the process specified, the difference between 
both calculations (using massless or a massive graviton in the m -¥ limit) lies on 
the amplitude \A4\ 2 . Therefore, the corresponding relation to Eq. (3.84) between 
amplitudes computed for each case can be written as 



lim X \M llv , a a\ 2 = lim M*V>^M a j 



m-s-0 



pol 



m-s-0 



= M* u M llv + ^M*/ M c a 



Z Mp, v , a 
pol 



J vl p , a 



(3.85) 



where M wa| g is the amplitude calculated using a massive graviton and M.uv,<x[i 
using a massless graviton. 

To obtain a smooth limit at the cross section level, the second term of Eq. (3.85) 
should be zero: 

'M£2 2 = 0. (3.86) 



To verify the smoothness of the limit we check that Eq. (3.86) is indeed true for the 
particular processes we are interested in. Take for example qq —> gG. Four Feynman 
diagrams contribute for such process, which are illustrated in Fig. 3.4. 
Calculating the amplitude for the first diagram, we obtain: 



(a) (b) 4 (c) (d) 




Fig. 3.4 Illustration of the four Feynman diagrams contributing for the process qq — > gG. The 
straight lines represent a quark or an anti-quark, the curly line represents the gluon and the double- 
curly line the graviton 
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= *°-v la T a ue«. (3.87) 

From the amplitude of the second diagram: 

M(b) = igsK 1 [ 7/(( _2£ 2 + + 7 „(-2t 2 + pi)„ 

8 (*2 - Pi) 
-2^(-2 fe+ ^)] (fe- p^a^ u e«e^ 

= ^^v la T a ue a a . (3.88) 
Similarly, for the third diagram: 

M (.c) = 3 !^H i}laT a ue a (389) 

And finally the forth: 

+ D l _ ll/tCl g{-p2, p\ + p2) + C~ 1 £'^,a0(-P2, Pi + P2)] « 

= ^ v 7 , d r a u i [| (-27 7q3 ) + (-*)T? a/! ] = 0. (3.90) 

It is now trivial to verify that the sum of all contributing diagrams is zero, exactly as 
expected from Eq. (3.86). If one follows the same strategy and computes \M\ 2 for 
the additional processes that can produce one graviton, the smooth limit is always 
obtained. Even for massive quarks in the initial/final state a similar calculation has 
been carried out and the smooth limit prevailed. 

To conclude, it would be significant to notice that the expressions for the cross 
sections displayed within the present study can be applied for the MGM and extra 
dimensional theories if the correct effective coupling are considered. For instance, 
for the MGM theory, Eqs. (3.78), (3.79) and (3.80) can describe correctly the theory 
if the general k parameter is replaced by the effective parameter set in Eq. (3.25). 
Similarly, the strategy can be applied to relate Eqs. (3.81), (3.82) and (3.83) with the 
ADD theory. 

Our subsequent goal is to investigate how to distinguish between MGM and ADD 
theories through phenomenological analyses at colliders, and moreover, if we can 
still identify signs of such theories considering the expected SM background. This 
will be the starting point for Chap. 5 following an introduction to hadron colliders 
and phenomenological simulations which will be presented in the next chapter. 
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Chapter 4 

The LHC and Collision Simulations 



The Large Hadron Collider (or LHC for short) was built to smash protons at a very 
high center-of-mass energy, in order to test predictions of different theories in particle 
physics and to confirm the existence of the Higgs boson. 

For the comparison between experimental data against theory predictions, simula- 
tions play a crucial role. In order to correctly make predictions of the consequences 
of particular BSM theories, it is necessary to use simulation tools. These packages 
can reproduce events expected at the hadron collider, which can be used to compare 
with real data. In this chapter, we thus present an overview of both the LHC itself 
and methods used to simulate events at hadron colliders. There are numerous refer- 
ences on this subject. We, however, shall base our discussion and follow the structure 
presented in Refs. [1-3]. 

We begin with an introduction of the LHC and its characteristics. Further ahead 
we present Monte-Carlo simulations of collisions and conclude by defining matrix 
element/parton shower merging schemes, which will be substantially employed in 
the following chapters for phenomenological analyses of new physics emissions. 



4.1 Introduction to the LHC 

The LHC was constructed at the "Organisation europeenne pour la recherche 
nucleaire" (CERN) near Geneva, in the tunnel which previously comprised the LEP 
experiment, the Large Electron-Positron collider that ran up to the year 2000. 
The LHC tunnel has 27 km of circumference, is buried approximately 100 meters 
under ground and is designed to collide protons with a maximum center-of-mass 
energy of 14 TeV. Collisions have started in the end of 2009, rapidly achieving the 
center-of-mass energy record of 2.36 TeV, exceeding even limits established by the 
Tevatron. 1 



The Tevatron is located at the Fermi National Accelerator Laboratory, and was built to collide 
protons and anti-protons with a center-of-mass energy limit of approximately 2 TeV. It has worked 
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Fig. 4.1 The LHC accelerator complex at © CERN [4] 



On March, 2010, the LHC started colliding protons with *Js — 7 TeV, beating its 
own previous record, and achieving the final center-of-mass energy established for 
the first run. 

To reach such colliding energy, the proton beam has to go through several acceleration 
steps, which are shown pictorially in Fig. 4.1 [4]. First, protons are produced in the 
L1NAC, and accelerated to 1 .4 GeV in the Booster. They are subsequently accelerated 
at the PS (proton synchrotron), and pass to the SPS to achieve the energy of 450 GeV. 
They finally move to the main ring to gain their final energy (at the current stage 
3.5 TeV each proton) and collide within one of the four main detectors: the Compact 
Muon Solenoid (CMS), the Toroidal LHC Apparatus (ATLAS), the LHCb or yet the 
Large Ion Collider Experiment (ALICE). The first two detectors are called "general- 
purpose" experiments. The LHCb is oriented towards analyses on the b-quark sector, 
CP violation and rare decays. ALICE however has been designed to focus on heavy 
ion physics and the study of quark-gluon plasma. 

In a nutshell, the goal for colliding protons is to test theories of fundamental interac- 
tions and their building blocks. Based on theoretical models we can predict chances 
of obtaining specific final states from the collision of two protons. With the LHC 
data, we are able to verify if theoretical expectations are actually fulfilled. With cross 



(Footnote 1 continued) 

from 1983 until September of 201 1, having many important discoveries on its list of achievements 
(including the discovery of the top quark in 1995). The large number of breakthroughs made by the 
Tevatron enabled the proposition of a even stronger particle collider, the LHC. 
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sections provided by theoretical calculations based on specific models, it is possible 
to compute the quantity of such events are expected at the detector using the specifi- 
cations for the collider. In reality, the complexity for linking experimental data with 
theoretical expectations is higher. However, we try to summarize the main aspects 
which will be employed within this study in what follows. 

The rate of collisions R and the luminosity C are related through the value for the 
cross-section a that can be calculated for a particular process: 

dN 

R = = £ct. (4.1) 

dt 

Here, ^M- is the number of events produced per second. 

For different experimental reasons, detectors cannot observe all events produced. We 
should thus account for the efficiency for observing a particular process. We obtain 



Nevents — o X L x £ , for L — J C dt , (4.2) 

where N even ts is the total number of events produced, L is the integrated luminosity 
provided by the experiment and £ is the efficiency of observing such process. 
The choice of building a proton collider was made aiming the achievement of high 
luminosities. The luminosity is controlled by the parameters of the experiment, the 
process however should be parameter independent. Anyhow, because the proton is 
composed by partons, one should consider parton distribution functions (also known 
as PDFs) to the computation of inclusive cross sections, considering the combination 
of different initial partons. 

The cross-section for a hard scattering process with two hadrons in the initial state 
[5] is 

o(P\, Pi) — X! / dXl dXl fit* 1 ' ^ fj( X2 ' °ij (/>!> P2< 2 2 /Mf> Q 2 /^ 2 r) 

ij 

(4.3) 

where Pi and P2 are the four- momenta of the initial hadrons and pi and P2 are the 
four-momenta of the partons present in the hard interaction. They are related by 

pi = xi Pi and p 2 — x 2 Pi- (4.4) 

The quark or gluon distribution functions at a factorization scale p f are exhibited 
by fi t j . The factorization scale can be qualitatively described as the resolution of 
the hadron which is being tested [6]. Q 2 is the characteristic scale of the scattering, 
and jiR is the renormalization scale. The parton-level cross section for the process 
envisaged to be analysed is and depends on Q as well as on the choice of the 
factorization and renormalization scales. 

From Eq. (4.4) it is trivial to observe that partons only carry a fraction of the four- 
momentum of the initial hadron, resulting on a parton-level center-of-mass energy 



52 



4 The LHC and Collision Simulations 



smaller than the total proton-proton «Js energy. They are related by r = x\ x 2 — s/s, 
where S — (pi + p 2 ) 2 and s — (P\ + P 2 ) 2 . We can define a new quantity, the 
differential parton luminosity, as: 

dLt , 1 1 Z" 1 r , 

—7^ = - 7: : „ — r / dxi dx 2 \xifi(x 1 )x 2 fj(x2) + (1 -*■ 2)1 S(x - x x x 2 ). 

(4.5) 

The cross-section will thus be given by: 

, , ^ f 1 dr \\dLi, 
a is) = > / — - — 
^J TQ r Is dr 

i,j u 

where the sum is performed over all pairs of partons. 

If one chooses the set of parton distribution function to be used, it is possible to 
obtain the differential luminosity - -4^ as a function of the partonic center-of-mass 
energy. Examples of such plots are shown for example in Ref. [7]. 
At the LHC with *Js — 14 TeV, we find that, in general, processes produced from 
a gluon-gluon initial state contribute more significantly than quark-quark or quark- 
gluon initial states. Therefore, processes having significant contribution at the Teva- 
tron collider may not be equally relevant at the LHC. 



(4.6) 



4.2 Predictive Methods at Colliders 



In order to correctly predict the consequences of a particular theory at colliders, it is 
necessary to define which observable one should focus on. In this section, we thus 
introduce the relevant observables for the present study. 



4.2.1 Kinematics and Jet Definitions 

Typically, in a collider experiment, the two incoming beams are assumed to be in 
opposite directions on the z-axis. The polar angle is 9 and the azimuthal angle is (p. At 
hadron colliders, the incoming beams have a spectrum of longitudinal momenta deter- 
mined by the PDFs. The center-of-mass of the incoming partons is usually boosted in 
comparison to the center-of-mass of the incoming hadron. Therefore, it can be useful 
to define final states in terms of variables proportional to longitudinal observables. 
Let us start by defining the transverse momentum pj of a particle. The four- 
momentum of a particle can be written as 



P 11 = (E, Px , Py , Pz ). 



(4.7) 
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The transverse momentum is defined to be 



PT = V Px + Pj = P sin6 - ( 4 - 8 ) 
Similarly we can define the transverse mass and transverse energy: 



mT — ^jpj+m 2 , Ej — E sin0. (4.9) 

Two further observables that will be largely employed within our phenomenological 
investigations are the missing transverse energy and the missing transverse momen- 
tum of an event: 

$T— — ^ Et, and $t= — Pt (4-10) 

visible visible 

which can be constructed from the sum of the transverse energy/momentum of visible 
particles within the event. This quantity is one of the smoking-gun observables at 
the LHC, since many new physics signatures are expected to have large E^j. 
Other significant observables are the rapidity y, 

y= H§^) =taniri (f)' (4 - n) 

= -lntan^j. (4.12) 



and the pseudo-rapidity 



In the coordinate system presented above, 8 is not Lorentz invariant, however, the 
pseudo-rapidity is. In the limit where m — > 0, the rapidity and the pseudo-rapidity 
of a particle are the same. It is also interesting to point out that we can calculate the 
pseudo-rapidity of an outgoing particle even not knowing its mass. 
Notice also that the four-momentum of a particle can be written in terms of the 
transverse mass, the rapidity, and the azimuthal angle: 

p 11 — (mj cosh y, /?7-sin</>, pjcoscp, mrsinhy). (4.13) 

Up to this stage we have only defined observables related to incoming or outgoing 
particles. However, sometimes we cannot completely experimentally differentiate 
outgoing particles. That is the case for light quarks and gluons, which due to color 
confinement, hadronize before being detected. These hadronic final states are usually 
numerous and appear as an agglomerate of particles, called jet ofhadrons. An illus- 
tration of the relation between parton-level outgoing particles and jets, i.e. hadronic 
final states, is represented in Fig. 4.2, extracted from [8]. 
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LO parton 



NLO parton 



parton shower 



hadrons 




Fig. 4.2 Illustration of the relation between partonic final states and hadronic final states, which 
gives identical jets in all the cases. Figure based on [8] 

To associate the hadronic final state with the most probable parton-level particle, it 
is necessary to define the jet. However, the reconstruction of a jet is not unique and it 
can vary depending on the parameters chosen. One should specify at least which are 
the particles belonging to the jet, and how to combine their momenta. These choices 
of parameters compose what is called as a jet definition. There are algorithms created 
to this aim, also referred as clustering algorithm. 

These are separated in two main classes: cone algorithms, and sequential recombina- 
tion algorithms. The first is based on the idea that the energy flow of the vertex is not 
modified by the branching 2 of the parton. It is a "top-down" approach in the sense 
that the algorithm defines the size of the jet based on one main particle. Sequential 
recombination algorithms, use a "bottom-up" approach, in which two particles are 
combined in one, sequentially, until only one particle (containing all the information 
of the others) is left. More details of each algorithm is given below. For more details, 
we refer to Refs. [8, 9]. 

Cone algorithms were the first type of jet algorithms proposed intended for the 
treatment of outgoing particles from e + e~ collisions [10]. It starts by classifying all 
the (final) particles according to their pj. The particle with the highest pj is defined 
as the central particle. The algorithm draws a cone of radius R around the central 
particle, in terms of the pseudo-rapidity r\ and the azimuthal angle <f>: 



where A?? = rji — r\ c and A0 = </>,• — <p c , in which i is the index referred to the 
s-th particle, and c to the central particle. Subsequently, the algorithm uses all the 
particles in the cone, and find the direction of the sum of the momenta of these 
particles. If the sum is the same as the momenta of the central particle, a stable cone 
is obtained. Otherwise, the sum of the momenta (of the particles) of the cone is made 
the new central particle and the iteration re-starts until the cone is found to be stable. 
Once the cone is stable, all the particles inside the cone can be removed from the list 
of particles of the event. The iteration stops once all the particles are identified as 
belonging to a specific cone. 



2 Once a parton is emitted, it can multiply and emit other partons. These are usually soft, or collinear 
partons, which can be approximately reproduced in perturbation theory. A more precise explanation 
is given further ahead within this chapter. 




(4.14) 
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The major disadvantage of using such kind of algorithm is the use of the pj of the 
particles to decide what will be the first central particle. The problem is that the 
transverse momentum is not a collinear safe observable, and it could lead to diver- 
gences. One could then use cone algorithms that do not rely on one particle to be 
the central, but perform several tries in which all particles are once central. Several 
stables cones would be obtained in this case, and an optimized merge would be done 
defining the jet. These types of cone algorithms avoid collinear divergences, but are 
not infra-red safe. A full solution to both problems was found and nominated Seedless 
Infrared Safe Cone algorithm [11], being the main cone algorithm used nowadays. 
Sequential recombination algorithms have a different strategy. They are based on a 
definition of "distance" between jets rather than a cone. The first algorithm using 
sequential combination was also developed for e + e~ collisions and called the kj 
algorithm [12]. In this case, a distance between two particles (i and j) is defined as: 

2min (e 2 , E 2 ) (1 - cos£) 

ytj = — q( > ( 4 - 15 ) 

where E is the energy of the particles, Q is the scale of the hard-scattering and 9 
is the angle separation between both particles. The two particles with the smallest 
distance, measured in terms of yn , are recognized and combined together giving rise 
to a "pseudo-particle". This strategy is repeated until the pseudo-particle and all the 
other particles in the event are separated by a larger distance than y cut . When this 
happens, the pseudo-particle is determined as the jet. 

The algorithm was adapted later on for hadron collisions [13, 14]. The distance 
between particles i and j can be fixed in terms of the AR: 




(4.16) 



where AR is defined similarly to the cone algorithm. It is also necessary to define 
the beam distance: 

diB = Pn- (4-17) 

Now, the algorithm searches for the smallest distance, either d/j or dj g ■ If the smallest 
distance is the one separating both particles i and j, the algorithm works as explained 
previously combining both to a pseudo-particle. Otherwise, the particle i is removed 
from the list of particles and the jet is defined. This algorithm is very often employed 
to define a jet. We will make use of it in our phenomenological analyses. 
A branch from the kj algorithm also very commonly used as of today is the so-called 
anti-kj algorithm [15]. Both kj and anti-A:7- algorithms have similar strategies, the 
difference lies on the distance definition. In this case, du and dig are defined as: 
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dij = 



diB — 



max 
1 



(p Z Ti ,p 2 T 



i) 



ARfj 
R 2 



(4.18) 



Once we have selected one method for defining a jet, it can be treated as an unique 
particle. We can thus define observables that are related to jets, such as the Hj . From 
the same class of the tfj observable, the Hj is defined as the sum of the pj of all 
the jets of an event: 



H T = 2 



\Pt 



(4.19) 



i=jets 



The Hj is an observable widely used in experimental searches due to its power on 
measuring the scale of the visible pj of an event. 

As the parton center of mass energy increases, hadron collision events with large 
jet multiplicities become more probable. As a result, accurate simulations require to 
correctly account for the presence of QCD radiation, which might modify predictions 
for these relevant observables that we have just presented. For specific processes, it is 
important therefore to account for extra orders in as when calculating cross section 
contributions. Subsequently, we demonstrate the necessity of such contributions. 



4.2.2 Next-to-Leading Order Cross Section Contributions 

The expression (4.3) is obtained through perturbative QCD. One can obtain cross 
section contributions of specific processes up to the lowest order in as- To achieve 
more accurate results, one should compute higher orders contributions. These can 
be calculated considering the addition of one extra particle (quarks, or gluons) to the 
process being evaluated. The first two additional orders in as are usually named next- 
to-leading order (NLO) and next- to-next- to leading order (NNLO) contributions. 
Take for instance the production of a Z boson [8] . The leading order term is pp —> Z. 
Equation (4.3) for this leading process is given by 

a(Pi, Pi) = X/ dx i dx2 /y(*2. A* 2 ) a ulz- (4- 2 °) 

To add the NLO term, besides taking loops into account, it is necessary to consider 
extra quarks and gluons produced in addition to the Z boson. The additional term 
is given by pp — »■ Z + X, where X is the emitted particle. A pictorial view of 
the process is presented in Fig. 4.3. Accounting for the NLO contribution, the cross 
section is given by: 
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Fig. 4.3 LO and NLO diagrams for pp — > Z production 

c(-Pi, Pi) = ^ j dx\ dx2 fi(xi,fj?) fj(x2, p 2 ) [trijt-z + o-^-z+x] • 

hi 

(4.21) 

In the previous sections, we have verified that the expected rate of observation of 
a process is related to its cross section, which is linked to observables such as the 
transverse momentum or the rapidity of a particle. The inclusion of higher order 
corrections on cross section calculations will thus increase the accuracy of results 
obtained for these relevant observables. This can be seen for example in the calcula- 
tion of the Drell- Yan cross section for the Z boson decaying into a pair of leptons. 
Consider the qq — > process, mediated by a Z boson and a photon. The tree- 

level cross section can be calculated [5]: 

a(qq -> = (q\ - 2Q q VtV q h{s) + (Aj + vj)(A 2 q + V^ 2 (s)) , 

(4.22) 

where s — (pi + pi) 2 , Pi is the four-momentum of the quark, p2 is the four- 
momentum of the anti-quark, Q f is the charge of the fermion, Vf — Tj — 
2Q f sin6 w , Af — Tj, — +1/2 for the up-quark/neutrino families and 
Tl — —1/2 for down-quark/electron families, fi and §2 are defined to be: 




(4.23) 



where Gf is the Fermi constant, a is the electromagnetic coupling, Mz and Tz are 
the mass and width of the Z boson. 

The differential lepton pair cross section da/dM^ t can be calculated, and from 
Eqs.(4.3) and (4.4) we know that s — x\X2S. Accounting for the different parton 
probabilities, the differential cross section is given by 
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= / dx 1 dx 2 y j [f q (x 1 )fg(x2) + Mx 1 )f q (x 2 )]--- r Y . (4.24) 
Jo „ clMfp 



da 

dM U q Ulrl u 

To find the consequences of NLO corrections to a relevant observable, it is possible 
to find the double differential cross section as a function of the rapidity distribution. 
Equation (4.1 1) can be written in terms of the parton momentum fraction: 

y = - log ( - ) , and consequently xi — -Jx e y , x 2 — Vr e~-\ (4.25) 
2 \x 2 J 

where r = Mf./s. The tree-level double differential cross section is thus 

= T-^fq(.Xi)fq{X 2 ) + fq(.Xi)f q (X2). (4.26) 



dM 2 u dy 3s 



'/ 



Once the next-to-leading order is considered, the double-differential cross section 
will have an extra order in as and divergences should be eliminated for the calculation 
to be finite. This is possible through different regularization schemes chosen [5]. For 
the process being considered, the (double differential) rapidity distribution is obtained 
for the LO, NLO and NNLO for instance by Anastasiou et al. in [16]. 
From such NLO and NNLO plots, we can notice one feature of these corrections at 
LHC: the large ctnlo/^lo factor, also known as the K-factor. Besides, 

&NNLO ONLO 

< . (4.27) 

onlo glo 

A similar statement can be made concerning the uncertainty bands, (which decreases 
from LO to NNLO expectations) suggesting that even higher order corrections would 
have minimal relevance at the LHC energy scale. 

Several tools exist to automatically compute cross sections at leading order of 2 —> n 
processes. These tools consider only tree-level diagrams, and the phase-space integra- 
tion is performed in a region where results are finite. Tools to determine automatically 
NLO corrections have a complexity much larger. Because loop computations involve 
divergences, it is not so trivial to obtain NLO correction from numerical calculations. 
A lot of progress has already been made for the automation of such tools [17-19], 
however, we shall not employ this variation of simulation tools within the present 
studies as they are not yet completed or publicly available. 

Nevertheless, leading order calculations have proven to give astonishing results. 
Distinct leading-order tools focused at different regions of phase-space have shown 
to be complementary providing overall very precise results when concerns the shape 
of relevant observables. Simulate events employing the correct tools is therefore the 
basis for obtaining accurate predictions at hadron colliders. Next, we present some of 
these tree-level tools, to subsequently introduce phenomenological analyses focused 
on graviton emission processes. 
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4.3 Monte- Carlo Collider Simulations 
4.3.1 Simulation Tools 

Usually, hadron collisions can be separated mainly in two categories: hard and soft 
collisions. The difference lies on how separated the resulting partons are. Well sepa- 
rated partons belong to the first category, whereas collinear partons are included in the 
second. Likewise, hadron collisions can be simulated using mainly two approaches, 
which are suitable for each category presented above. The first is the matrix element 
technique, in which amplitudes are associated to a particular process, allowing for a 
cross-section calculation and event generation after the necessary integration over the 
phase-space. The second is the parton-shower description that simulates successive 
parton splitting, taking into account the emission of QCD radiation and, possibly, 
hard scatterings. 

There are several Monte-Carlo software that can be used for both descriptions. For 
instance, Sherpa [20], CalcHEP [21], CompHEP [22], MadGraph/MadEvent 
[23-26], Alpgen [27], HELAC [28] or Whizard [29] are based on the matrix- 
element technique. 

To exemplify how simulation within the matrix element technique works, let us use 
one of the above matrix element generator: MadGraph/MadEvent. Initially, the 
process envisaged to be calculated should be specified, i.e. both initial and final state 
particles should be stated. The theory which is likely to be used for the compu- 
tation of the processes has to be declared. As a result, MadGraph generates the 
Feynman diagrams and computes matrix element amplitudes at a randomly chosen 
point of phase-space using helicity amplitude routines. The implementation of helic- 
ity amplitudes and wave-functions have first been performed by the HELAS package 
[30-33] and have recently been automatized [34], as describes Appendix C. These 
matrix element amplitudes are subsequently used by MadEvent, which computes 
cross-sections (or decay widths) for the particularly chosen process, and generates 
simulated events for hadron colliders. 

This approach is very flexible due to the fact that, in principle, any theory or model 
of particle physics can be used within the code, once it has been implemented in 
the event generator. However, the description is only valid at determined regions of 
phase-space, where final-state partons are very well separated from each other. 
When protons collide at high energies, there is a large probability that partons will 
continue basically in the same direction with small transverse momentum. There is 
however some probability that particles with large transverse momentum are gener- 
ated from the collision. This is so-called hard scattering, and can be diagrammatically 
seen in Fig. 4.4. When the scattering is hard, calculations must be done using per- 
turbative QCD [5]. The lowest order is the one used in matrix-element simulators 
described so far. 

Nevertheless, soft processes are more likely to be generated. A emission is soft when 
the particle emitted has very a small fraction of the energy of the particle that emitted 
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it. A particle defined to be collinear when the angle between two particles is very 
small, such as they seem to be almost collinear in the literal sense of the word. 
In the limit where final-state partons are soft, or collinear, one has to take into 
account parton splitting and QCD radiation. These are also refereed as "underlying 
events" and encompass events with initial and final state radiations and interactions 
of remaining partons of the initial proton. In fact, when a collision happens, partons 
can interact with other produced partons (or partons reminiscent from the collision) 
by radiating quarks and gluons. These radiated particles can themselves interact with 
other partons to produce their own QCD radiation. This forms a chain, or a shower 
of radiating quarks and gluons, also known as parton shower. 
The parton shower approach is a treatment of QCD parton splitting which is a per- 
turbative approximation and valid up to the order of some cut-off value. Above this 
value, which is usually fixed in terms of the pj of the final particles, the probability 
of the splitting is low and the hard scattering approach is more suitable to describe 
such process. To briefly summarize how the parton shower approach works, we shall 
strictly follow Refs. [3, 5]. Consider the parton branching illustrated in Fig. 4.5. The 
parton a splits into b and c. We assume that 

t = Pi » Pi and Pi (4-28) 

where t is also known as the virtuality of the incoming parton. The separation angle 
between both outgoing partons is 9, and the energy fraction is z: 
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Fig. 4.5 Illustration of a 
parton branching based on 
Ref. [5] 



In the soft and collinear limit, the cross section to describe a TV -parton process (e>7v) 
is modified by the splitting of one parton as: 



as dt /> 

d(JN+i = da N - P ba (z)dz 

2tt t 



(4.30) 



where P is the Altarelli-Parisi splitting function [5]. As a consequence, sequential 
parton splitting can be simulated. 

Suppose multiple parton branching, such as this represented in Fig. 4.6. It can be 
shown that the probability of a parton evolving from t' to t without resolvable radi- 
ation, i.e. without emitting other partons above a cutoff e, is given by A(f)/A(r'), 
where A(f) is: 



A,-(0 = exp 



z 



' dt' 



dz— — Pji(z) 



(4.31) 



The A;(r) is the probability summed over all possible processes i — > j. It is called 
Sudakov form factor and is the basis of Monte-Carlo simulation for the parton-shower 
approach. These algorithms employ the Sudakov form factor to control sequential 
splittings. They consider the probability of evolving from t\ to without resolvable 
radiation: 

xkH (432) 

where 1Z is a randomly generated number between [0, 1]. If the value of ti is smaller 
than Q 2 , the next branching should have a fraction z — x%/xi of the initial momen- 
tum and the next splitting should be generated. If this is the case, the probability 
distribution for generating a fraction z of the initial momentum will be proportional 
to (a s /2jz)P(z), where P(z) is the splitting function. Otherwise, no further split- 



Fig. 4.6 Illustration of 
sequential parton splitting 
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ting occurs. One can obtain the fraction of the initial momentum after a splitting by 
solving [5]: 



where 1Z' is randomly generated number between [0, 1], and e is the cutoff for 
resolvable radiation. 

This results on a cascade emission of partons up to the step where no more branching 
occurs. At this stage, partons are converted into hadrons following a hadronization 
description. This is used to simulate the creation of stable particles which will be the 
true observed final state at colliders. The hadronized particles will interact directly 
with the detector, and the resulting description is described as being "at detector 



Different Monte-Carlo simulation tools are also available to simulate events fol- 
lowing the parton shower and hadronization descriptions. These are for instance, 
Pythia [36], Herwig [37], and Ariadne [38]. 

The choice of which description to use is of course made depending on the process 
to be analysed. Collisions with a significantly low center-of-mass energy, processes 
that include few jets in the final state are very well described by using matrix-element 
simulation tools. However, as the parton center-of-mass energy increases, events with 
numerous jets in the final state become more likely. It is thus necessary to take into 
account the presence of QCD radiation using parton shower simulation. Because 
parton-shower approach underestimates the radiation of hard jets, the best solution 
is to use both descriptors at the same time. Nevertheless, some care is necessary in 
the limit where both tools are employed. At this region, one could consider the same 
event twice. For avoiding this problem one can use a matching, or merge scheme 
[39-41], which will consistently divide the space-space in two regions and specify 
the description to be used. A more detailed explanation on matching schemes is 
introduced in the following section. 



4.3.2 Matrix-Element/Parton-Shower Matching Scheme 

The idea is to correctly simulate events using Monte-Carlo techniques and account 
for jets produced by initial state radiation. Generally, there are two ways of simulating 
(N + l)-jet events: i) requesting a final state with N partons, where hard emissions 
could lead to the extra jet, or ii) from (N + 1) partons final state, in which soft or 
collinear radiation evolved to generate the extra jets. 

As we have previously stated, the necessity of employing both prescriptions to simu- 
late events with multiple jets in the final state comes from the fact they act at different 
regions of the phase-space. Parton shower descriptions are accurate for collinear/soft 
jet emission, but they can easily break down above a given scale where matrix ele- 
ment approach should be used. However, if both approaches are employed without 




(4.33) 



level". 
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Double Count 




shower shower shower Z + 1-jet 

Z + 1-jets Z + 2-jets + hard gluon 

Fig. 4.7 Pictorial example of a possible double counting of events with different multiplicities if the 
parton shower approach is considered at the same time as the matrix element technique without the 
employment of matching scheme. In the first diagram, it is shown one example for the production 
of a Z boson in association with one parton. The diagram in the middle exemplifies the production 
of a Z boson in addition to 2-partons final state. Finally the third diagram shows the production of 
a Z + 1 -parton process in which the parton radiates a hard-gluon, resulting on the same signature 
as the previous diagram 



some care, we could double count events with different multiplicities. One example 
of an event that could be double counted is illustrated in Fig. 4.7. The first diagram 
shows the production of a Z boson with one additional parton, which showers. The 
second diagram represents the Z boson produced with two partons in the final state. 
The third diagram however illustrates the production of a Z boson in addition to 
one parton which radiates and creates a hard gluon. Considering both approaches 
at the same time would double count events represented by the second and third 
diagrams. We have therefore to apply matching techniques in order to define, on an 
event-by-event basis, which prescription should be employed. 
Different approaches have been proposed to merge the matrix element with parton 
shower descriptions. One can cite for instance the CKKW scheme [39, 40] and the 
MLM scheme [42, 43]. In general these procedures differ in three main aspects: 
the jet definitions used in the matrix element scheme, the employment of accep- 
tance/rejection of the jet configuration and in the way that jets are vetoed inside the 
parton shower approach. Nevertheless, these approaches are based on a similar tech- 
nique, which is basically to divide the phase-space into two regions. The separation 
between these regions, let us name it Q ma tch , should be expressed in terms of the 
AR or pr variables specified on the previous section, or of the k j distance defined 
byEqs.(4.16)and (4.17). 

Here, we will focus on the MLM approach, using the software 
MadGraph/MadEvent for matrix element simulations, and Pythia for the parton 
shower description. To be more precise, let us describe how the MLM description 
works. We use the k j algorithm to cluster jets, and we shall strictly follow the guide- 
line presented in Refs. [41, 43]. A description in details can be found within these 
references. 

The initial stage is to generate processes with N partons in the final state, i.e. to 
produce pp — >• X + A' -jets samples using the matrix element descriptor. They are 
defined following the kinematic cuts: 
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PT > pf", M < rjmax, ARij > Rmin, (4.34) 

where pj and \ r]\ are the transverse momentum and pseudo-rapidity of the final state 
partons and ARji is their minimal separation in the (rj, 0) plane. The parameters 
with a subscript "min", or "max" must be the same for all N partons. 
The tree structure is set according to the color-flow of the event, and the chosen 
scales of as, using the k j value calculated at each vertex. The factorization scale for 
the parton densities is given by the hard scale of the process, i.e. Q 2 — m 2 x + p 2 ~ x . 
Subsequently, the shower is generated by the parton-shower descriptor and evolved 
up to the hard scale of the process, which was defined to be Q 2 . A jet algorithm is then 
applied to the partons produced, specifying the final number of jets. These jets are 
determined by the following parameters: the cone size R c i us , the minimum transverse 
energy Ej US and the pseudo-rapidity r) ch,s . To be consistent with the division of the 
phase space, and to cover all the available range, it is necessary that R c i us > R m in< 

clus < a pclus > min 

'Imax — 'max <" lu c y _ Pj ■ 

The matching scheme starts now. Beginning from the hardest parton, the closest jet 
in the (r/, <p) plane is chosen. If the distance between the parton and the center of 
the jet is smaller than R m i n , the jet is matched. This matching test is done in order 
to match each jet to a parton in the event. If that happens, the event is considered to 
be fully matched. Events which are not fully matched are rejected. Also, for events 
that pass the matching test, it is additionally required no extra jets to be present. 
That will avoid the double counting of the events that are already predicted, and 
more accurately described, by the N + I sample. However, for the highest parton 
multiplicity, these events are not rejected (considering no such N +1 sample will be 
available), maintaining the full inclusiveness character of the sample. 
It is valuable to state that the distribution of observables should not depend on the 
parameters defined on the scheme (such as R, rj or Q ma tch)- Moreover, the transition 
from one description to the other should be smooth. 

The determination of the parameters used for the matching should depend on the 
process to be calculated. To demonstrate the utility of matching schemes, we shall 
present results for the process pp — >• Z + 0, 1, 2-jets. We use this specific process, 
since it is one of the main background channels for the phenomenological analyses 
that we present for graviton emission. The Z boson could decay into a pair of neutri- 
nos, and the signal expected would be missing transverse energy in addition to jets. 
The first parameter we need to define is the matching scale. This is not a physical 
observable, therefore the results should not depend on it. However, the choice of 
such scale should be done smartly, in order to divide coherently the phase-space. 
If the cutoff is chosen too low, divergences in the matrix element simulation could 
appear. However, if the scale is chosen too high, there could be some regions where 
the parton shower would not be able to simulate events. 

The intention is to have a smooth transition from one description to the other. The 
differential jet rate is the best variable to measure how smooth this transition may be. 
This variable is defined as the kj scale at which an event passes from the (N + l)-jet 
to a N jet configuration while the algorithm cluster the showered partons. In reality, 
all the matching procedure do is to group showers according to the k j scale chosen: 
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if the jet is within the kj distance, it belongs to a specific parton and grouped together 
with all the jets belonging to the same parton, otherwise it does not. In consequence, 
in a transition N + \ — > N, all the events with final number of jets smaller than N + 1 
should be below the cutoff, and events with multiplicity equal or larger than N + 1 
should lie above the scale chosen. In Fig. 4.8, we present differential jet rate plots for 
pp — >■ Z + 0, 1, 2 -jets, where we have used Q ma tch — 30GeV, Pythia for the 
parton shower simulation and MadGraph/MadEvent for the generation of events 
based on matrix element calculation. For the employment of the matching scheme 
as well as for the creation of the plots, we have made use of MatchChecker [3], 
code designed for the interpretation of matching analyses. 

The first graph of Fig. 4.8 shows the 1^0 differential jet rate. We can see that the 
curves of Z + 1-jet and Z + 2-jets samples are placed above the cutoff, while the 
Z + 0-jet sample is placed below. For 2 — > 1, we have Z + 2-jets above the Q matc h 
scale, while Z + 0-jet and Z + 1-jets situate below. One can see therefore in the 
three curves presented, the exact division where the straight pink line lies. In the left 
side of the line, mainly events are determined by parton shower approach, and in 
the right side using the description based on matrix element amplitudes. There is no 
need to focus at the 3^-2 differential jet rate, because the transition is given fully 
by the parton shower. 

As it has been stated, the choice of the Q ma tch scale should not modify overall 
results, i.e. the shape of the observed distributions. To prove this statement we show 
in Fig. 4.9 the differential jet rates for several choices of the cutoff scale. As expected, 
the curves lie above each other, predicting the same final result. 
Finally, we can show distributions for some relevant observables. For instance 
Fig. 4. 10 exhibits the differential cross section as a function of both the pj of the Z 
boson (left) and the pj of the highest energetic jet (right). One can easily see how 
parton shower and matrix element descriptions are complementary. At low values for 
the pt, the main contribution comes from lower multiplicities, basically computed 
through the parton shower description. As the energy increases, contributions from 
multiple jets in the final state become more significant and therefore matrix element 
description is more suitable to correctly account for them. 




LogPrfferentialjet Rate 1—0) LogtOirTenential J etfiate 2- It . og(Wferenia1J n us.v i . 



Fig. 4.8 Differential jet rates for the process pp — > Z + 0, 1 , 2-jets. The cutoff can be identified 
by the straight pink line. The curves in the left of the cutoff are produced using parton shower 
simulation. The curves in the right by matrix element simulation. Here, we can easily see the 
smoothness of the transitions 
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Fig. 4.9 The choice of the Qmatch scale should not modify the shape of any physical distribution. 
This statement is demonstrated in this plot of the differential jet rate of a 2 — > 1 transition for 
different values of the Qmatch scale. The curves shown are the final result, i.e. the sum of the 
contribution given by each multiplicity. The process being considered is Z + 0, 1, 2 jets 



One substantial use of the matching technique is that it predicts very accurate shapes 
for distributions. The reason is that besides leading order processes being considered 
with the employment of the matrix element description, QCD radiation is taken 
into account within the parton shower approach. Therefore, even being a leading 
order simulation, the matching technique takes into account real emission processes. 
It therefore provides very precise shapes of the distributions comparable to NLO 
predictions. However, NLO calculation are still necessary for the prediction of the 
K-factors. 

In order to verify how precise matching predictions are compared to NLO calculation 
for graviton emission processes, we have performed complete phenomenological 
analyses. Next chapter will therefore be focused on such phenomenological studies. 




Fig. 4.10 Differential cross section as a function of the pj of the Z boson, left, and of the pj of 
the highest energetic jet, in the right. The process being considered is Z + 0, 1 , 2 jets, at the LHC 
with v^s = 14 TeV, Q match = 20 GeV and CTEQ6L1 PDF 



References 



67 



References 



1. J.M. Campbell, J.W. Huston, W.J. Stirling, Hard interactions of quarks and gluons: a primer 
forLHC physics. Rep. Prog. Phys. 70(1), 89-193 (2007, hep-ph/0611148) 

2. E. Halkiadakis, in Proceedings for TASI 2009 summer school on physics of the large and the 
small: introduction to the LHC experiments, p. 32 (2010, 1004.5564) 

3. S. de Visscher, Observability of an unconventional two-Higgs-doublet model at the LHC. 
PhD thesis, Universite Catholique de Louvain, 2009 

4. C. Lefevre, Illustration of the CERN accelerator complex. (2008) http://cdsweb.cern.ch/ 
record/ 1260465 

5. R.K. Ellis, W.J. Stirling, B.R. Webber, QCD and collider physics. (Cambridge University 
Press, Cambridge, 2003) 

6. F. Maltoni, T. McElmurry, R. Putman, S. Willenbrock, Choosing the factorization scale in 
perturbative QCD, p. 25 (2007, hep-ph/0703156) 

7. C. Quigg, LHC Physics Potential vs. Energy: Considerations for the 201 1 Run. FERMILAB- 
531FN-0913-T, p. 33 (2011, 1101.3201) 

8. G.P. Salam, Elements of QCD for hadron colliders, Proceedings of 2009 European School of 
High-Energy Physics— ESHEP 2009 (2010, 1011.5131) 

9. G.P. Salam, Towards jetography, Eur. Phys. J. C67, 637-686 (2010, 0906.1833) 

10. G. Sterman, Jets from quantum chromodynamics. Phys. Rev. Lett. 39(23), 1436-1439 (1977) 

11. G.P. Salam, G. Soyez, A practical seedless infrared-safe cone jet algorithm. J. High Energy 
Phys. 2007(05), 086 (2007, 0704.0292) 

12. S. Catani, Y. Dokshitzer, M. Olsson, G. Turnock, B.R. Webber, New clustering algorithm for 
multijet cross sections in e + e~ annihilation. Phys. Lett. B 269(3^1), 432^t38 (1991) 

13. S. Catani, Y. Dokshitzer, M.H. Seymour, B.R. Webber, Longitudinally-invariant kj -clustering 
algorithms for hadron-hadron collisions. Nucl. Phys. B 406(1-2), 187-224 (1993) 

14. S.D. Ellis, D.E. Soper, Successive combination jet algorithm for hadron collisions. Phys. Rev. 
D 48(7), 18 (1993, hep-ph/9305266) 

15. M. Cacciari, G.P. Salam, G. Soyez, The anti-k_t jet clustering algorithm. J. High Energy Phys. 
2008(04), 12 (2008, 0802. 1 189) 

16. C. Anastasiou, L. Dixon, K. Melnikov, F. Petriello, High-precision QCD at hadron colliders: 
Electroweak gauge boson rapidity distributions at next-to-next-to leading order. Phys. Rev. D 
69(9) (2004, hep-ex/03 12266) 

17. G. Ossola, C.G. Papadopoulos, R. Pittau, CutTools: a program implementing the OPP reduc- 
tion method to compute one-loop amplitudes. J. High Energy Phys. 2008(03), 042 (2008, 
0711.3596) 

18. R. Frederix, S. Frixione, F. Maltoni, T. Stelzer, Automation of next-to-leading order compu- 
tations in QCD: the FKS subtraction. J. High Energy Phys. 2009(10), 003 (2009, 0908.4272) 

19. V. Hirschi, R. Frederix, S. Frixione, M. Vittoria Garzelli, F. Maltoni, R. Pittau, Automation 
of one-loop QCD computations. J. High Energy Phys. 2011(5), 64 (201 1, 1 103.0621) 

20. T. Gleisberg, S. Hoeche, F. Krauss, A. Schaelicke, S. Schumann, J. Winter, SHERPA l.alpha, 
a proof-of-concept version. J. High Energy Phys. 2004(02), 28 (2003, hep-ph/031 1263) 

21. A. Pukhov, CalcHEP 2.3: MSSM, structure functions, event generation, batchs, and generation 
of matrix elements for other packages, pp. 1-31 (2004, hep-ph/0412191) 

22. E. Boos, V. Bunichev, M. Dubinin, L. Dudko, V. Edneral, V. Ilyin, A. Kryukov, V. Savrin, A. 
Semenov, A. Sherstnev, CompHEP 4.4-automatic computations from Lagrangians to events. 
Nucl. Instrum. Methods Phys. Res., Sect. A 534(1-2), 250-259 (2004, hep-ph/0403113) 

23. J. Alwall, P. Demin, S. de Visscher, R. Frederix, M. Herquet, F. Maltoni, T. Plehn, D.L. 
Rainwater, T. Stelzer, MadGraph/MadEvent v4: the new web generation. J. High Energy 
Phys. 2007(09), 028 (2007, 0706.2334) 

24. J. Alwall, M. Herquet, F. Maltoni, O. Mattelaer, T. Stelzer, MadGraph 5: going beyond. J. 
High Energy Phys. 2011(6), 37 (201 1, 1 106.0522) 

25. T. Stelzer, W. Long, Automatic generation of tree level helicity amplitudes. Comput. Phys. 
Commun. 81(3), 357-371 (1994, hep-ph/9401258) 



68 



4 The LHC and Collision Simulations 



26. F. Maltoni, T. Stelzer, MadEvent: Automatic event generation with MadGraph. J. High Energy 
Phys. 2003(02), 027 (2003, hep-ph/0208256) 

27. M.L. Mangano, F. Piccinini, A.D. Polosa, M. Moretti, R. Pittau, ALPGEN, a generator for 
hard multiparton processes in hadronic collisions. J. High Energy Phys. 2003(07), 001 (2003, 
hep-ph/0206293) 

28. A. Cafarella, C. G. Papadopoulos, M. Worek, Helac-Phegas: A generator for all parton level 
processes. Comput. Phys. Commun. 180(10), 1941-1955 (2009, 0710.2427) 

29. W. Kilian, T. Ohl, J. Reuter, WHIZARD: Simulating multi-particle processes at LHC and 
ILC. Eur. Phys. J. C 71(9) (2007, 0708.4233) 

30. K. Mawatari, Y. Takaesu, HELAS and MadGraph with goldstinos. Eur. Phys. J. C 71(6), 1-10 
(2011) 

31. K. Hagiwara, J. Kanzaki, Q. Li, K. Mawatari, HELAS and MadGraph/MadEvent with spin-2 
particles. Eur. Phys. J. C 56(3), 435^147 (2008, 0805.2554) 

32. K. Hagiwara, K. Mawatari, Y. Takaesu, HELAS and MadGraph with spin-3/2 particles. Eur. 
Phys. J. C 71(1), 1-14 (2011) 

33. H. Murayama, I. Watanabe, K. Hagiwara, HELAS: HELicity amplitude subroutines for 
Feynman diagram evaluations, http://www-lib.kek.jp/cgi-bin/img_index7199124011 (1992, 
KEK-91-11) 

34. P. de Aquino, W. Link, F. Maltoni, O. Mattelaer, T. Stelzer, ALOHA: Automatic Libraries 
of Helicity Amplitudes for Feynman diagram computations, Comput. Phys. Commun. 183, 
2254-2263 (2012,1108.2041) 

35. R.C. Group, Measurement of the inclusive jet cross section using the midpoint algorithm in 
Run II at the Collider Detector at Fermilab (CDF). PhD thesis, University of Florida, 2006 

36. T. Sjostrand, S. Mrenna, P. Skands, A brief introduction to PYTHIA 8.1. Comput. Phys. 
Commun. 178(11), 852-867 (2008, 0710.3820) 

37. G. Corcella, I.G. Knowles, G. Marchesini, S. Moretti, K. Odagiri, P. Richardson, M.H. Sey- 
mour, B.R. Webber, HERWIG 6: an event generator for hadron emission reactions with inter- 
fering gluons (including supersymmetric processes). J. High Energy Phys. 2001(01) (2001, 
hep-ph/0011363) 

38. L. Lonnblad, Ariadne version 4 — A program for simulation of QDC cascades implementing 
the colour dipole model. Comput. Phys. Commun. 71(1-2), 15-31 (1992) 

39. F. Krauss, Matrix elements and parton showers in Hadronic interactions. J. High Energy Phys. 
2002(08), 015 (2002, hep-ph/0205283) 

40. S. Catani, F. Krauss, B. R. Webber, R. Kuhn, QCD matrix elements + parton showers. J. High 
Energy Phys. 2001(11), 063 (2001, hep-ph/0 109231) 

41 . J. Alwall, S. Hoche, F. Krauss, N. Lavesson, L. Lonnblad, F. Maltoni, M. Mangano, M. Moretti, 
C. Papadopoulos, F. Piccinini, S. Schumann, M. Treccani, J. Winter, M. Worek, Comparative 
study of various algorithms for the merging of parton showers and matrix elements in hadronic 
collisions. Eur. Phys. J. C 53(3), 473-500 (2007, 0706.2569) 

42. F. Caravaglios, A new approach to multi-jet calculations in hadron collisions. Nucl. Phys. B 
539(1-2), 215-232 (1999, hep-ph/9807570) 

43. S. Hoeche, F. Krauss, N. Lavesson, L. Lonnblad, M. Mangano, A. Schaelicke, S. Schu- 
mann, Matching parton showers and matrix elements. Proceedings of HERA and the LHC: A 
workshop on the implications of HERA for LHC Physics CERN— DESY 2004/2005 (2006, 
hep-ph/0602031) 



Chapter 5 

Phenomenology of Graviton Production 



In exact sciences, the word "phenomenology" is employed to relate empirical obser- 
vations of phenomena with theories that attempt to describe it. This definition applies 
to particle physics as well. Comparison between theory and experiment is achieved 
through computer simulations. High-energy events are generated following the prob- 
ability distributions expected from a theory, assuming certain initial conditions. Once 
the expectations of a particular theory are known, it is easy to test it and to constrain 
parameters. 

The goal of this chapter is to perform phenomenological analyses on graviton emis- 
sion processes at hadron colliders, taking into account theories which have been 
previously presented as a possible extension to the SM. 

We start with an introductory analysis of graviton emission within processes with 
one jet in the final state [1, 2]. > The events are characterized at colliders through 
its missing transverse energy signal. Two different theories are considered for such 
analysis: the theory with one massless graviton (MGM), and the theory with large 
extra dimensions (ADD), both presented in Chap. 3. The theory with warped extra 
dimensions (RS) is not expected to generate missing energy signature, since the first 
mode of the graviton should quickly decay into SM particles. Therefore, we do not 
consider it for this specific analysis. 

We subsequently present a full phenomenological treatment of graviton emission 
within processes with multiple jets in the final state [3]. 3 For the analysis, most types 
of theories predicting spin-2 particles are included and separated in three represen- 
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tative classes. Predictions of most relevant observables using different simulation 
methods are displayed, compared and validated against next-to-leading order QCD 
calculations. 



5.1 Graviton Emission in Mono- Jet Searches 

The specific four-dimensional model discussed in Chap. 3 is rather different at the 
theoretical level from the ADD extra dimensional model. Besides the fact that the 
graviton propagator has a different expression for each theory (due to the dependency 
of the propagator with the mass of the graviton), gravity could become strong in four 
dimensions if there would be a large hidden sector interacting gravitationally with 
the SM. 

This large hidden sector could contribute to the renormalization group running of 
Newton's constant. In Chap. 3 we have shown that for approximately 10 32 new par- 
ticles, it would be possible to have M p i(l TeV) ~ 1 TeV. For this reason, at the 
effective level of the theory, we could replace the coupling constant of Eq. (3.9) by 
the running Planck mass: 

K i 

~==——, and M p i (fit) ~ /i*, (5.1) 
2 Mp/(/i t ) 

where ~ 1 TeV. The masses of the new particles determine the energy scale 
at which the running of Newton's constant starts. We could choose for example 
ra, ~ 100 GeV. If that is the case, most of the running of the Planck mass occurs at 
the order of 1 TeV. This suggests that the running of the reduced Planck mass could 
be effectively approximated by a step function: 

If fi < 1 TeV, M p i(n) = 2 - 4 x !0 18 GeV, (5.2) 
while for \i > 1 TeV, M pi (ji) ~ 1 TeV. (5.3) 

As a result, the phenomenology expected for the MGM theory is very peculiar, since 
strong gravity breaks at the TeV scale for the MGM, while it does not for the ADD 
theory. 

On the other hand, an infinite tower of KK excited modes of the graviton (coupled 
to the SM fields) are expected within the framework of the ADD theory. Subsequent 
modes of the graviton should have a mass gap of 1 /R, with R being the radius of the 
extra(s) dimension(s). For plausible values of R, excited modes of the graviton have 
a very compact discrete mass spectrum, i.e. subsequent modes should lie very tight to 
each other in absolute values of mass. In order to obtain an effective description of the 
full KK graviton tower, it is necessary to sum over all the states, or approximately, to 
integrate over the entire range of the mass spectrum. The mass integration is defined 
in Eq. (3.40), which results on the (expectation of) observation of one KK graviton 
with mass m. 
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In both theories, the graviton particle is expected to be identified at hadron colliders 
as transverse missing energy or momentum. Therefore, our first goal is to make a 
phenomenological analysis of processes that include invisible particles [ 1 , 2] . A com- 
parison of the signal against the SM expectation for missing energy, which is mainly 
translated into the production of pair of neutrinos, will give us a useful benchmark. 
One extra motivation to the current phenomenological investigation has been pre- 
sented in Sect. 3.4. We would like to investigate if there are phenomenological conse- 
quences of the propagator discontinuity, even if no discontinuity has been observed 
with analytical calculations at cross section level. For this reason, we have chosen to 
scrutinize graviton production in addition to one-jet in the final state, including all 
the processes discussed in Eqs. (3.78)-(3.83). 

Mono-jet (i.e. 1-jet + missing energy) signature is one of the most familiar signal 
of several BSM theories, including both scenarios presented here. It will be crucial 
therefore to establish whether LHC could differentiate between theories. Moreover, 
in case of non-observation of mono-jet signals, one could use data from the LHC to 
constrain parameters of both theories. Particularly, it could set a limit on the value 
of the scale at which quantum gravitational effects become important. 
The production of jets with large Ej in addition to a graviton G can emerge from 
the parton sub-processes q + q ^ g + G, q + g ^ q + G, q + g ^ q + G, and 
g + g — »• g + G. To simulate and generate events for both theories, we use the matrix 
element description through the MadGraph/MadEvent package. The ADD theory 
was implemented in MadGraph/MadEvent version 4 by K. Hagiwara et al. [4] and 
is supported by the package since 2008. 

We follow the same strategy to implement the MGM theory according to our neces- 
sity. First, internally to the MADGRAPHcode (or more precisely to the HELAS code), 
modifications in the graviton wave function routines are made in order to properly 
deal with massless spin-2 particles. In addition, helicity routines are modified and 
the massless spin-2 particle propagator is added to the code. Finally, a model file 
for the MGM is created (including particles, interactions and vertex definitions). In 
Appendix A the implementation is detailed, and checks and validation tests that have 
been accomplished are displayed. 

The running of the Planck mass is implemented with a Heaviside step function in the 
cross section. The Planck mass for collisions at the parton level with Vs > /i* is given 
by M(/i*) = and for less energetic parton level collisions we can approximate 
M p i — > oo. Because it is quadratic in energy, most of the running happens near to 
the scale in which gravity becomes strong, /i*. Therefore, the approximation of the 
running of the Planck mass by a Heaviside step function can be very accurate. 
The discussion begins with the results from predictions obtained by the MGM theory. 



5.1.1 Mono- Jet for the MGM Theory 

As explained, the graviton itself is not detectable and it appears as missing energy. 
One could think that if gravity becomes strong at the TeV scale, the Tevatron would 
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Fig. 5.1 Number of events per bin considering a luminosity of 100 fb 1 . The background consid- 
ered is Z(— > vv) + jet. Figure extracted from [1] 
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Fig. 5.2 Number of events per bin at Tevatron considering a luminosity of 1.1 fb~ l . The back- 
ground is Z(— > vv) + jet. Notice that the y-axis is not displayed using a logarithmic scale in order 
to facilitate the analysis. Figure extracted from [ 1 ] 



have already felt its effect through Drell-Yan processes. However, in this scenario, 
the Drell-Yan process is extreme suppressed since virtual gravitons are decaying 
into the large hidden sector, and thus, the main signature would be that of the present 
analysis, therefore initially compared our predictions to the irreducible background 
process pp{p) — > v v + jet [5-7] both for the LHC and Tevatron. 
We display results for the number of events expected at the LHC in Fig. 5.1 and at 
the Tevatron in Fig. 5.2. The analysis is made at the LHC for *fs — 14 Tev with L = 
100 fb -1 of integrated luminosity and at the Tevatron for «J~s — 1.96TeV with 
L = 1.1 fb -1 . The CTEQ6L1 PDF is employed. In order to have the background 
reduced, we choose a high-fy cut and require the jet to be central: 
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LHC : pi > 500 GeV, W et \ < 3 

It (5A ^ 
Tevatron : p J T et > 100 GeV, \r]J e '\ < 3. 

From Fig. 5. 1 it is possible to deduce that there is a possibility for the LHC to identify 
signals of a massless graviton if p* < 4 TeV. 

The distinctive indication of such theory would be the peak observed at E™' ss — 
p*/2. The reason for this specific location of the peak is understood. Let us fix 
VI = Mo, which is the invariant mass of the final state. The pr of the jet, and in this 
case also the Ej , is ^ sin 9, where 9 is the angle of the jet with respect to the beam. As 
9 varies, for any fixed Mq, the Ej changes and can be in the range < Ej < Mo/2. 
Now, let us concentrate on the cross section for the graviton emission processes. 
As Mq increases, the cross section should diminish fast due to the reduced phase 
space and the rapid decrease in the parton densities. As a consequence, the cross 
section will have a maximum value when Mo is minimum. However, in our case 
Mo > M p i (/a) ~ p+ below which our cross section is zero by definition. This leads 
the Ej distribution to have a maximum at /i*/2, as observed in both Figs. 5.1 and 5.2. 
Therefore, the threshold is clearly a consequence from the imposition of the running 
of the Planck mass through the Heaviside step function. 

From Fig. 5.2 we can infer that the Tevatron reach is quite limited. Combining results 
from the figure and these presented by the CDF Collaboration [7], we find that if 
/i* ~ 500 Gev, experimental hints should already have emerged at Tevatron. This 
limit is also obtained when cosmic rays data is used [8]. For p t > 600 Gev, we 
observe that it would be very hard to differentiate the MGM theory from the SM 
background. 

We now present a comparative analysis of the results obtained with those predicted 
for the ADD theory at the LHC. 



5.1.2 Mono- Jet for the ADD Theory 

We employ two different techniques to simulate ADD at the LHC. Primarily, a 
modified version of the code developed by T. G. Rizzo for a different investi- 
gation [9] is adapted to the present analysis. Results are compared to a Mad- 
Graph/MadEvent fully simulated sample. Same cuts given in Eq. (5.4) are used, 
and the PDF CTEQ6.6 [10] is applied. 

The results of our direct comparison of the ADD predictions with those of the MGM 
theory assuming that Mp/(/i*) = 1 TeV can be found in Fig. 5.3. 4 Here, we observe 
that both new physics models predict mono-jet Ej distributions which are reasonably 
different in shape. The decreasing curve from MGM mono-jet spectrum above the 



4 Reprinted from Physics Letters B , v. 682/4-5 , Xavier Calmet, Priscila de Aquino, Thomas G. Rizzo, 
Massless versus Kaluza-Klein gravitons at the LHC, 446-449, Licence number 3125370259028 
(2010), with permission from Elsevier. 
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Fig. 5.3 The number of mono-jet events expected as a function of the Ej' n distribution at the 
14TeV LHC, with L = 100 fb" 1 of integrated luminosity. A Ef n cut of 500 GeV has been 
applied as well as the \i] je '\ < 3 cut to require the jet to be central. The yellow histogram is the 
expected SM background (Z — > vv), and the red and higher histograms are for the ADD model 
with the number of extra dimensions being 2, 3, 4, etc. The lower black curve is for the MGM model 
with = 1 TeV. The parameters of the ADD curves are in each case adjusted by varying their 
associated Planck scale to produce the same result as does the MGM model at Ej' n = 500 Gev, 
in order to compare the relative shapes for these distributions. Figure extracted from [2] 



\/s — Mi, threshold is lighter than the corresponding ADD distributions for any 
number of extra dimensions. 

Comparing Figs. 5.1 and 5.3, we see the peak in the distributions associated with 
the MGM model, arise from the cross section threshold at VS = /i*, which are 
non-existent for the distributions associated with the ADD theory. In the latter case, 
the Et distribution decreases monotonically, while the threshold in the MGM model 
cross section naturally leads to peak in these corresponding distributions. Therefore, 
the threshold would be the main point in distinguishing predictions of these two 
classes of models at the LHC. Moreover, because we have an indication of a cutoff 
in the cross section at a fixed value of Vs, this feature could allow us to extract the 
cutoff value directly from the LHC experimental data. 

Notice from Fig. 5.3 that the parameters of the ADD theory are adjusted so that its 
curve starts at the same point as that of the MGM theory. This is done in order to 
fairly compare shapes of distributions obtained from both theories. A small differ- 
ence between the black (MGM) and coloured curves (ADD) is observed. It is not a 
powerful difference that would allow to experimentally differentiate between both 
BSM theories. However, both curves lie above the background, and if observed, they 
would annunciate possible signals of new physics. 
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5.2 Graviton Emission in Multi-Jet Searches 

We have analysed so far phenomenological consequences of graviton emission in 
combination with one parton level light jet. However, events with only one jet have a 
relatively low probability to appear at the LHC, considering its large center-of-mass 
energy. Besides, we have exclusively investigated graviton emission processes in 
which the signature is composed of invisible particles. There are however theories, 
such as the RS theory, that predict the decay of the spin-2 particle (into a pair of 
leptons for example). 

Our subsequent goal is therefore to perform a complete analysis on graviton emis- 
sion processes, including processes with more than one jet in the final state, and 
incorporating a full set of theories that contain different classes of spin-2 particles. 
Additionally, we aim on comparing our results against SM background expectations 
at the LHC and at Tevatron. 

The search for this kind of new physics at the TeV scale is today among the major tasks 
of the LHC. Experimental analyses of such signatures require accurate simulations 
not only to design efficient signal selection strategies but also to extract information 
on new physics parameters, such as masses, coupling strengths, structures and spin. 
The main idea is thus to include a wide class of models that involve spin-2 particles, 
such as: 

• weakly coupled (Planck scale suppressed) massive graviton models, with infinite 
number of nearly continuously distributed KK modes of the graviton and the 
overall integrated effects are significant at the TeV scale, represented here by the 
ADD model; 

• TeV scale massless graviton models, in which the graviton effects are negligible 
below the fundamental scale of quantum gravity, i.e. the TeV scale, due to the 
running of the gravitational constant, e.g., in theories with many particles that 
interact with the SM particles only gravitationally. Represented here by the MGM 
theory; 

• strongly coupled (TeV scale) massive graviton models, in which the graviton KK 
modes are widely separated. The lowest excitation gets mass at the TeV scale 
and can decay into the SM particles within detectors. This class of theories is 
represented here by the RS model. 

To show and validate the potential of MadGraph/MadEvent implementation 
choose to study the most important signatures in representative models of the three 
classes above. For the ADD and MGM models the gravitons couple weakly to the 
SM particles and appear as invisible particles at colliders, as previously discussed. 
The main searching channel is pp(p) — > G + jet(s), in which at least one accom- 
panying hard jet is needed. For the RS model, resonance production pp(p) — >■ G is 
expected, with G decaying into e + e~ or for example. For completeness we 

additionally consider graviton in mono-jet channel to compare the RS model with 
the ADD and MGM results. 
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Several tools exist nowadays that are able to perform such simulations [11]. Here, 
we focus on MadGraph/MadEvent, where a spin-2 particle implementation at the 
matrix element level exists [4]. In the appendix A we present the modifications made 
on such implementation to incorporate massless spin-2 particles. For the analysis 
with multiple jets in the final state, modifications in the massive spin-2 routines 
are still necessary in order to obtain a fully automatic event generator. First, Mad- 
Graph version 4 can only generate Feynman diagrams with up to 4-point vertices, 
and for this reason, diagrams containing 5-point vertex (e.g., 4 gluons-graviton) have 
to be added by hand. Although it was fine for previous studies on graviton produc- 
tion with up to 2 jets [12, 13], it is inconvenient for graviton production with higher 
multiplicity jets. Secondly, gravitons are densely distributed in the ADD model and 
hence, it requires a special treatment for the propagator. These modifications have 
been successfully performed and are detailed in Appendix B. 
As the parton center-of-mass energy increases, hadron collision events with large 
jet multiplicities become more probable. As a result, it is necessary that accurate 
simulations correctly account for presence of QCD radiation, which could modify 
leading order predictions for relevant observables. 

Recall that in Chap. 4 it has been explained that there are mainly two approaches to 
simulate hadronic collisions: the matrix element technique and the parton shower 
approach. In parton shower MC programs like PYTHIA [14] and HERWIG [15], addi- 
tional jets are usually obtained in the collinear and soft approximation. Hard and 
widely separated jets are poorly described in this approach. On the other hand, tree- 
level fixed order calculations can provide reliable predictions in the hard region, 
while failing in the collinear and soft limits. For this reason, to correctly take into 
account QCD radiation, matrix element simulators must be employed together with 
parton shower approach making use of matching description detailed in Sect. 4.3.2. 
In this study, we use the kr-MLM matching scheme implemented in MadGraph/ 
MadE vent and interfaced to PYTHIA for parton shower and hadronization. In the kj 
matching scheme, matrix element multi-parton events are produced with a minimum 
separation k j cutoff, as described in Eq. (4.16) for Q ma tch — 2mhr ^ or ever Y event, 
the final-state partons are clustered according to the k j algorithm. Only clusterings 
corresponding to the Feynman diagrams provided by MadGraph are allowed. In 
order to closely mimic the behaviour of the parton shower, the kj value for each 
clustering vertex corresponding to a QCD emission is used as renormalization scale 
for a s in that vertex. For the central hard 2 — s- 1 or 2 — ► 2 process, the transverse 
mass rrij — p\ + m 2 of the particle(s) produced in the central process is used as a 
factorization scale, as well as a renormalization scale. 

Subsequently, this event is passed to the parton shower MC simulator. Before 
hadronizing or decaying, the final partons are clustered into jets using the kj algo- 
rithm with a jet cutoff of Q ma tch > The jets are then compared to partons. 
If dij (parton, jet) is smaller than the cutoff Q ma tch> they match. The event is only 
accepted is all the partons have matched with jets. For events with parton multiplicity 
smaller than the highest multiplicity, the number of jets must be equal to the number 
of partons. 
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Fig. 5.4 Leading order diagrams for the calculation of graviton emission. The curly lines represent 
a gluon, the straight lines a quark or anti-quark, and the double wavy line represents the graviton 

Before establishing our phenomenological analyses, let us present one more detail 
of the current investigation. Our aim is equally to compare matched predictions 
against NLO calculations. For the RS theory, the leading order process is the zero-jet 
sample, since spin-2 particles should decay into a pair of leptons. For ADD and 
MGM theories, the leading order process is necessarily the 1-jet sample, because 
missing transverse energy cannot not be an observable if there is no visible particle 
to allow the invisible to be tagged. Therefore, we start by presenting the investiga- 
tion carried out for RS models in fully inclusive samples for pp — >■ G + n-jets, 
where n = 0, 1, 2. Subsequently, exhibit a semi-inclusive analysis for pp —> G + 
n-jets, where n — 1, 2, 3, in which all three classes of theories presented above are 
investigated. 



5.2.1 Fully Inclusive Graviton Production for the RS Theory 

The first process to be explored is pp — >• G + X, where the graviton will eventually 
decay into leptons. Let us start the discussion with the NLO-QCD corrections to such 
process. 

5.2.1.1 NLO-QCD Corrections 

In Chap. 4 a brief discussion on NLO-QCD corrections to Z boson production within 
the SM has been presented. Currently, we follow the same strategy for the compu- 
tation of NLO-QCD corrections to the production of a spin-2 particle within the RS 
theory. We begin with the process Pi + P2 — > G + X, where Pi and P2 are the initial 
hadrons, and X is an arbitrary hadronic final state. 

The QCD hadronic cross section can be expressed as Eq. (4.20). To add NLO terms 
is necessary to consider loops and real emission of QCD radiation. Equation (4.2.2) 
represents the cross section expression including first order terms. 
The difference compared to the calculation presented in Sect. 4.2.2 is that for Z 
boson production only quark-anti-quark initial states are relevant, whereas for gravi- 
ton production it is necessary to consider also the gluon-gluon initial state. Fig. 5.4 
presents the diagrams relevant for the LO process, Fig. 5.5 shows real emission dia- 
grams contributing for the NLO calculation and finally Fig. 5.6 exemplifies diagrams 




Fig. 5.5 Real emission diagrams contributing for the NLO calculation of graviton emission process 

Fig. 5.6 Loops diagrams 
contributing for the NLO 
calculation of graviton emis 
sion process 




relevant for loop corrections. The diagrams contributing to first order corrections can 
be summarized as: 

► G + g; 

► G (one-loop); 

► G + q- 

■ G + q\ 

■ G + g; 

■ G (one-loop). 

write Eq. (4.21) explicitly for graviton emissions [16-18]: 

<t{Pu Pi) = Y,Jdxi dx 2 fi(xi) fj{x 2 ) t [afjl G + a s (A) ^jl G+x ] ■ (5.5) 
ij 




1 

9 
9 
9 
9 



We can re- 
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where ^\j\ G = Oq 5(1 — z) is the Born partonic cross section, as — as/Air and 

(7; P G+x are the contributions from loops and QCD emission. The scaling variables 
are: 



pi = X\P\ 



P2 = x 2 Pi 



S — X\X 2 S — TS, z 



(pi + piY 



Ml 



TS T 



(5.6) 



For graviton emission we have: 
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(5.7) 



where n t is the number of quark flavours considered and the "plus distribution" is 
defined as: 

(5.8) 



dzh + (z)f(z)= dzh(z)(f(z)- /(D) 



For the computation of the NLO cross section, we implement a small code based on 
Eqs. (5.5) to (5.7). A change of variable is performed to facilitate the integration: 



re y \ x 2 



~Jre 



-y- 



T — X\X 2 - 



(5.9) 



We find the jacobian and change the integration limits: 
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Table 5.1 LO and NLO results for the integrated cross section of a graviton emission process 
within the RS theory at the LHC with center-of-mass energy of 14 TeV. The first mode of the 
graviton couples to the SM fields according to Eq. (3.9). Here, it is chosen ^ = j = 3 TeV, 
= MR = m G» n f = 5 and employed the CTEQ6L1 PDF [19] 



rng GeV a L0 (pb) ctnlo (pb) K-factor 

91 4.59 x 10 3 8.73 x 10 3 1.90 

500 1.97 x 10 2 3.38 x 10 2 1.71 

800 5.77 x 10 1 9.67 x 10 1 1.67 

1000 2.99 x 10 1 4.95 x 10 1 1.65 

1200 1.67 x 10 1 2.74 x 10 1 1.63 

1500 0.78 x 10 1 1.25 x 10 1 1.61 



Table 5.2 LO and NLO results for the integrated cross section of a graviton emission process within 
the RS theory at the LHC with center-of-mass energy of 7 TeV. The first mode of the graviton couples 
to the SM fields according to Eq. (3.9). Here, it is chosen 4- = j =3TeV,f/,p = fig = mg,rif =5 
and employed CTEQ6L1PDF[19] 



m G GeV a to (pb) a NLO (pb) K-factor 



91.2 


1.61 x 10 3 


3.21 x 10 3 


1.99 


500 


2.62 x 10 1 


4.61 x 10 1 


1.76 


800 


5.24 


8.73 


1.67 


1000 


2.18 


3.52 


1.62 


1200 


0.99 


1.56 


1.57 


1500 


0.34 


0.51 


1.53 



dx\ I dx2 — / dz ~2 / dy (5.10) 

TQ Jtq/x\ J T Z J\Og^T /z 

For the leading order distribution, the integral is easy to solve. We present it therefore 
as a reference. For instance, the leading order contribution for the gluon-gluon initial 
state is exhibited in Eq. (5.7). The cross section is thus given by: 



cr(gg/qq -* G) = — r / dy fi (V^)fj (>/V-»), (5.11) 

10 ilogVTB 



where (i, j) — (g, g), (q, q), (q, q). The integral over partonic fractions presented 
just above are also named partonic luminosities. 

As a result, we obtain the K-factor for graviton emission within the RS theory. 
Tables 5.1 and 5.2 display values for the K-factor calculated for different center- 
of-mass energy at the LHC. Notice that the factorization and renormalization scales 
are set to be the mass of the graviton. 

The aim is now to obtain matching predictions for fully inclusive processes and 
compare the K-factor obtained, defined to be (Jnlo with the normalization 
factor, or N-factor, defined to be ctnlo /&match- 
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Log(Differential Jet Rate 1 -» 0) Log(Differential Jet Rate 2 -» 1 ) 

Fig. 5.7 Differential cross section as a function of the differential jet rate for 2 — »• 1 and 1 — s> 



5.2.1.2 Matching Predictions 

For the matching scheme we have chosen to use pp — > G + n-jets where n = 0, 1, 2 
(n — 2 samples are inclusive). The parameters employed for the RS theory are the 
same as for the NLO calculation. The analysis has been accomplished for: (a) the 
LHC with *Js — 14 Tev and mc — 1 Tev; (b) the LHC with yfs — 7 Tev and mc — 
500 Gev; (c) and for the Tevatron with yfs — 1.96 Tev and both niQ — 500 Gev and 
(d) niQ — 1 Tev. At generation level we have employed Q matc h = 50 Gev, and no 
cuts have been used. 

First, we make sure the matching is properly working. To this aim we plot the 
differential 2 — >• 1 and 1 — >• jet rates. Results for the LHC with «Js — 7 Tev for a 
graviton with mc — 500 Gev are displayed in Fig. 5.7. Notice the smoothness of the 
passage from when parton shower is employed to when the matrix element technique 
is used. All samples both for the Tevatron and the LHC have similar shapes for the 
differential jet rate and present as good smooth passage. 

Subsequently, in order to exemplify the predicting capability of matching technique, 
we present the most relevant distributions at the LHC in Fig. 5.8. These plots exhibit 
the differential cross section as a function of its particular observable. The integrated 
cross sections can be calculated, and are displayed in Tables 5.3 and 5.4. We normal- 
ize matched cross sections with leading order values presented in Tables 5.1 and 5.2 
for the corresponding values for the mass of the graviton, which are the underlined 
row. Results for both the K-factor (<tnlo/&lo) an d N-factor {(JNLo/o'match), as 
well as the cross sections obtained, are exhibited in Tables 5.3 and 5.4. 
Observe that values for the K- factor are very close to the normalization factors (or N- 
factors, for short). In fact, the cross section after matching analysis is expected to be 
similar to the leading order computed for pp — > G. Leading order calculation already 
includes the production of extra radiation through the PDFs (parton distribution 
functions). However, matching calculation considers explicitly the production of 
extra jets in the final state, which leads to a modification on the behaviour of relevant 
observables, such as the missing Pj, pseudo-rapidity and Ht distribution. As a result, 
we expect that matched sample, normalized to NLO calculation, will produce better 
predictions for observables which involve extra radiation than a pure parton shower 
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Fig. 5.8 £7--MLM matching results for a pp — > G + rc-jets with n = 0, 1, 2 full inclusive sample 
of the RS model with A = 3 Tev. The matched results are given by the red and blue curves for 
mi = 1 Tev at the LHC with V? = 14 Tev, and mi = 500 Gev at the LHC with *Js = 1 Tev, 
respectively. No cuts were applied in these samples. Figure extracted from [3] 



approach. In order to verify such statement, we follow with a complete analysis of 
pp —¥ G + n-jets, in which the leading order process is now mono-jet, and the 
matching is performed up to 3-jets in the final state. 



5.2.2 Semi-Inclusive Analysis for Graviton Emission Processes 

For different classes of theories, such as the ADD and MGM, it is more interesting 
to explore missing energy signature in combination with jets. We present therefore 
a similar investigation to these of the previous section. However, the leading order 
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Table 5.3 K-factor O(NL0)/&{L0) an d normalization factor (J{NLO)/ CJ match for inclusive G pro- 
duction in the RS model at the Tevatron, from pp — ► G NLO QCD calculations and parton shower 
matching with the matrix elements of pp — ► G + n jet(s) with n = 0, 1, 2, A = 3 Tev, mc = 0.5 
and 1 Tev for ^/s = 1.96 

m a Tevatron 

o-LO (pb) (tnlo (pb) K-factor Qjnatch (pb) N-factor 

500 Gev 1.039 1.416 1.36 1.059 1.34 

1 Tev 0.020 0.031 1.55 0.021 1.52 



Table 5.4 K-factor <7{NL0)/&(L0) and normalization factor o"(jv/,o)/ cr ma/c/i for inclusive G pro- 
duction in the RS model at the LHC, from pp — > G NLO QCD calculations and parton shower 
matching with the matrix elements of pp — >• G + ;i jet(s) with n = 0, 1, 2, A = 3 Tev, me = 0.5 
and 1 Tev for *Js = 7 and 14 Tev LHC, respectively 





LHC 












07,0 (pb) 


otvlo (pb) 


K-factor 


& match (pb) 


N-factor 


500 Gev 


26.238 


46.155 


1.76 


24.642 


1.87 


1 Tev 


23.803 


40.316 


1.69 


24.362 


1.65 



process should currently be the emission of one graviton in addition to one jet. An 
example of such diagrams are presented in Fig. 3.4 for the quark-anti-quark initial 
state. Similar diagrams exist for the gluon-quark and gluon-gluon initial states. 
As previously performed, we consider a general theory containing a graviton that 
couples to the SM fields following Eq. (3.9) and in which the Feynman rules are 
displayed in Fig. 3.1. We summarize the relation between the general parameter k 
(explicitly revealed in the lagrangian) and the effective scale of the three representa- 
tive classes of theories currently being analysed: 



MGM theory 
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= 2 ~ 






1 


K 


1 




A 


= 2~ 


M p i 
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ADD theory — = - = =— (5.12) 
RS theory 

For the ADD model, one extra integration over KK the mass spectrum is performed 
according to Eq. (3.40), which results the effective scale to be of the TeV order. 
The implementation of the running of the Planck mass for the MGM theory has 
been performed during the generation of the events. The technique to implement 
the running has changed from the analysis exhibited in Sect. 5.1. Previously, we 
had applied a cut on Vs. However, now we have extra jets in the final state. When 
these extra jets are soft/collinear, we could get non-physical enhancements below 
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Fig. 5.9 Example of loop dia- 
grams contributing for NLO 
calculation of the process 
pp ->■ gG 





the threshold, which would have nothing to do with the Q 2 of the scattering. A 
solution would be to apply a cut in all partonic Mandelstam variables: Vs, v— f 
and V— u. The problem is that cuts on v— ? and -J —u are not physical, and could 
spoil the infra-red divergence cancellation when NLO corrections are calculated. We 
have thus decided to apply a minimum cut of H™ n > M(^*) ~ /x* in the sum of 
transverse momentum of the jets, which is stronger than the combination of the cuts 
on Mandelston variables. In other words, a Heaviside step function-like behaviour is 
still added to the cross section. However, now for collisions with Hj < M(/i*) ~ p t , 
gravity contributions are so weak that the cross section can be approximated to zero. 
To fix the parameters related to the theories, we consider the up-to-date experimental 
constraints for each representative class [6, 7, 20-23]. For the ADD model, we 
choose A = 5 Tev at the LHC, A = 1 Tev at the Tevatron, with 8 — 2, 4, 6. For 
the MGM model, we set M(/z*) ~ /x* = 1 Tev and 2 Tev. For the RS model, we 
take niQ — \ Tev and 100 Gev, with A = 3 Tev. The latter value chosen for the 
mass of the graviton is clearly excluded by experimental searches. We have decided 
to employ it in order to check how RS theory behaves at low graviton mass limit. 
The idea is to check that no strange behaviour comes from the discontinuity of the 
propagator. 



5.2.2.1 NLO-QCD Corrections 

We present the NLO corrections to pp — >• G j + X. Such calculation has been 
performed by several authors, including [12]. Here, we summarize the main aspects. 
The NLO calculation for the cross section consists on computing loops and real QCD 
emission, as presented in the previous section. The leading order contributions are 
now the diagrams represented in Fig. 5.5. Examples of loop diagrams contributing for 
this specific calculation are presented in Fig. 5.9. Again, we can separate the calcula- 
tion in three parts: the computation of virtual corrections, real emission contributions 
and the collinear term. Virtual correction arises from the interference of the born and 
one-loop amplitudes. Real-emission contributions appear from the production of a 
QCD light jet, and the collinear term is the finite reminder of the factorization of 
collinear singularities into the PDFs. 
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Table 5.5 Integrated cross section, K-factor unlo I^LO an£ i N-factor cr^Lo/cr ma tch for semi- 
inclusive G production at the LHC with y/s = 14 Tev, from pp — s> G NLO QCD calculations and 
parton shower matching with the matrix elements of pp — ► G + n — jet(s) with n = 1, 2, 3. For 
the ADD theory, the effective scale is set as 5 Tev, and S is the number of extra dimensions. For the 
MGM theory, 1 Tev and 2 Tev stand for the values chosen for /i*. We consider A = 3 Tev for the 
RS theory, in which the mass of the first mode of the graviton is either «ig = 0.1 or mg = 1 Tev 
respectively 



LHC 





<tlo (fb) 


a NLO (fb) 


K-factor 


Pmatch (fb) 


N-factor 


ADD <5 = 2 


1.67 x 10 2 


2.76 x 10 2 


1.65 


1.35 x 10 2 


2.05 


ADD S = 4 


3.83 x 10 1 


6.22 x 10 1 


1.62 


2.67 x 10 1 


2.33 


ADD S = 6 


1.45 x 10 1 


2.36 x 10 1 


1.63 


9.47 x 10° 


2.49 


MGM 1 Tev 


2.27 x 10 3 


1.08 x 10 4 


1.86 


6.89 x 10 3 


1.56 


MGM 2 Tev 


7.19 x 10 1 


1.52 x 10 2 


2.11 


9.24 x 10' 


1.64 


RS 1 Tev 


9.65 x 10 2 


1.59 x 10 3 


1.65 


7.98 x 10 2 


1.99 


RS 100 Gev 


4.72 x 10 3 


8.51 x 10 2 


1.80 


4.70 x 10 3 


1.81 



The authors of [12] have created a code to compute NLO corrections for graviton 
emission considering the ADD theory. In [3], we have employed and extended their 
code to compute NLO corrections accounting for the three classes of theories. 
Concerning the parameters chosen for the calculation, in order to suppress SM back- 
grounds in LHC graviton searches [5] we require: 

p T > 500 GeV, \n je, \ < 4.5, and p J T et > 50 GeV. (5.13) 

Jets are defined by the kj algorithm, with the resolution parameter set to D = 0.6. 
At the Tevatron, we require 

Pt > 120 GeV, \r] Jet \ < 1.0, and p{ et > 150 GeV, (5.14) 

following the settings considered in a recent CDF study [6, 7]. Notice p J T et is the 
minimum pr required for the leading jet. Here, all jets are defined by the k j algorithm 
with D — 0.7, and are required to satisfy \r)j\ < 3.6 and p J T > 20 GeV. A second jet 
with pr > 60 GeV is vetoed. 

Again, the SM background considered is the irreducible pp —> vi l>i + jets. The 
idea is to use the background as a reference curve. A detailed experimental search 
would have to consider other types of processes for the SM background, such as the 
production of W boson decaying to a lepton and a neutrino, where the lepton is not 
tagged. 

For the LO/NLO calculations, CTEQ6L1/6M PDFs are employed together with the 
corresponding values for the strong coupling as- The renormalization and factor- 
ization scales are set to the transverse momentum of the graviton pj . We vary the 
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Table 5.6 Integrated cross section, K-factor unlo I^LO an£ i N-factor a-^Lo/cr ma tch for semi- 
inclusive G production at the Tevatron with y/s = 1.96 Tev, from pp — »• G + j NLO QCD 
calculations, and parton shower matching with the matrix elements of pp — s> G + n — jet(s) with 
n = 1, 2, 3. For the ADD theory, the effective scale is set as 1 Tev, and 8 is the number of extra 
dimensions. There is no reason for the computation with the MGM theory, since it has been proved 
not possible to be detected at the Tevatron. We consider A = 3 Tev for the RS theory, in which the 
mass of the first mode of the graviton is either mc =0.1 or mo = 1 Tev respectively 



Tevatron 





olo (fb) 


ONLO (fb) 


K-factor 


'match (fb) 


N-factor 


ADD 5 = 2 


5.22 x 10 2 


6.39 x 10 2 


1.22 


2.66 X 10 2 


2.40 


ADD S = 4 


1.82 x 10 2 


1.77 x 10 2 


0.97 


7.16 x 10' 


2.47 


ADD 5 = 6 


8.24 x 10 1 


6.91 x 10 1 


0.83 


3.59 x 10 1 


1.92 


RS 1 Tev 


1.91 x 10- 1 


1.20 x 10-' 


0.62 


6.17 x 10~ 2 


1.95 


RS 100 Gev 


2.31 x 10 2 


3.37 x 10 2 


1.46 


1.95 x 10 2 


1.73 



scales in the range between pj /2 and 2pj to show the uncertainties related to the 
scales. 

On Tables 5.5 and 5.6 we display LO and NLO values obtained for the integrated cross 
section at the LHC and at Tevatron. The K-factor is calculated and also exhibited. 
In addition, we have computed the differential cross section as a function of the 
following distributions: 

(a) Transverse missing energy (for ADD and MGM theories) or transverse energy 
of the graviton (for RS theory); 

(b) Transverse momentum of the first and second highest energetic jets; 

(c) Pseudo-rapidity of the the graviton; 

(d) Pseudo-rapidity of the first and second highest energetic jets; 

(e) The Hj distribution, which is defined by Eq. (4.19). 

Before showing and interpreting results obtained for each relevant observable, let us 
describe the implementation of the matching scheme. 

5.2.2.2 Matching Predictions 

Similarly to the NLO calculations, we analyse pp — >• G + n-jets, where n = 1, 2, 3 
with the leading order process being n— 1 . We employ MadGraph/MadEvent for 
the parton-level generation, and shower using the PYTHIA package. In the event gener- 
ation with MadGraph/MadEvent, CTEQ6L1 PDF [19] is employed, which fixes 
the values of as at the Z mass. The choices of renormalization and factorization 
scales are defined to be the transverse mass of the graviton. 

The calculation of the cross section as well as the event generation is performed 
using the following matching scales: 
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LHC Q^ n > 45 GeV, Q match > 50 GeV 
Tevatron Q™f n > 20 GeV, Q match > 30 GeV. 



(5.15) 



In Tables 5.5 and 5.6 we display the obtained matched values for the integrated cross 
section at the LHC and at Tevatron. The N-factor is calculated and exhibited in order 
to allow the comparison against the K-factor previously computed. 

5.2.2.3 Comparative Analysis and Interpretation of the Results 

Let us start the discussion with a comparative analysis of the values obtained for the 
K-factor against N-factor. Tables 5.5 and 5.6 exhibit both values for the three classes 
of theories and parameters specified in the beginning of this section. Notice we do 
not observe consistent values at the Tevatron as we observe at the LHC. Recall the 
settings used to compute such numbers: 

&7-MLM matching calculation: 

• Process: pp G + n-jets, n = 1, 2, 3. 



• PDF: CTEQ6L1 
NLO/LO calculation: 

• Process: pp — >• G + 1-jet (+ 1 jet for NLO). 

• Scales: /if — /i« — pj 

• PDF: CTEQ6L 1 /CTEQ6M 

In Chap. 4 we have seen that depending on the collider, different processes can 
have different differential luminosity contributions. At the LHC, the sub-process 
that contribute the most for graviton emission is gg — > Gg. At Tevatron this is not 
true. From Fig. 5. 10 we can see that the most contributing sub-process is qq Gg, 
which carries a large Bjorken \x\ (equal or larger than 0.5). For large values of \x | , the 
valence quarks, which are dominating, become smaller at NLO. This results in a NLO 
PDF different from the LO PDF. Notice that in our calculation, we have employed 
NLO PDF for NLO calculation, and LO PDF for LO calculation. If we also take into 
account the difference in as used to compute NLO and LO cross sections, it will be 
reasonable to obtain a smaller value for the NLO cross section. This is the cause for 
such small values obtained for the K-factor at Tevatron. 

Now, let us concentrate on the N-factors for the Tevatron samples. Notice they are 
of order 2 for a heavy graviton, due to the factorization and renormalization scales 
chosen for the NLO/LO and matching calculation. For light graviton production, 
no difference at the cross section level is expected because the scales chosen for 
&7--MLM matching and NLO/LO are similar, i.e. (pj) 2 — (nij) 2 — (pt)g + m c- 
However, for heavy graviton (pj) 2 < {m\) 2 . 
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Table 5.7 Results for the integrated cross section obtained from the LO process pp ~* Gg within 
the RS theory for mg = 1 Tev and A = 3 Tev. At the LHC it has been employed p 3 j > 500 Gev, 
and \rf et \ < 5 cuts and at Tevatron p S j > 150 Gev, and \rf et \ < 1. The PDF is set to CTEQ6L1 





&(PF = = 


P?) 


a(p F = Pr = mf) 


Ratio 


LHC 


0.957 pb 




0.803 pb 


1.19 


Tevatron 


0.203 fb 




0.087 fb 


2.33 



This can be easily observed if we perform a brief validation check in which we 
compute values for the LO cross section of the process pp — > Gg within the RS 
theory. The idea is to employ these two different scales in the calculation both at 
the LHC and Tevatron, and compare results obtained. The cuts applied within the 
analysis are the following: 

LHC p ] r > 500 GeV, |?>> et | < 5, CTEQ6L1. 

T (5.16) 
Tevatron pf > 150 GeV, |7?i et | < 1, CTEQ6L1. 

Table 5.7 presents the outcome. At the LHC we have a factor of approximately 1.2 
for the ratio cr(pj)/cr(mj). At Tevatron this ratio is much higher: 2.3. Notice that 
the ratio between both scales can be interpreted as N-factor/K-factor ratio of our real 
analysis. 

From both Tables 5.5 and 5.6 we obtain N-factor/K-factor ~1.2 for the LHC, and a 
larger ratio of approximately 2.5 at the Tevatron for the heavy graviton cases. These 
results demonstrate the reason why we obtain different values for the K-factor and 
N-factor for the current semi-inclusive analysis. They also demonstrate the consis- 
tency of our results. 

Fig. 5.10 The differential 
cross section as a function of 
the partonic center-of-mass 
energy V? for the process 
pp — > Gg considering the RS 
theory at the Tevatron. This 
plot demonstrates separately 
each sub-process contribution. 
The analysis is performed with 
m.Q = 1 Tev, A = 3 Tev and 
the cuts: $ > 150 GeV and 
|r/i et | < 1 
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Fig. 5.11 NLO/matching comparison for the ADD model at the LHC (y/s = 14 Tev) with A = 
5 Tev. The NLO results are given by the red, green and purple bands for 8 = 2, 4, 6, respectively. 
Matched results are given by the black curves in which we employ the same parameters as the for 
the NLO. The dominant Z — > vv background is also shown as a reference by the blue curve. The 
kj -MLM matched curves are normalized by the NLO results, and the normalization factors can be 
found in Table 5.5. Figure extracted from [3] 



5.2.2.4 Distributions of Relevant Observables 



Figures 5.1 1, 5.12, 5.13, 5.14, 5.15 and 5.16 exhibit the differential distributions for 
the comparison between NLO-QCD calculation and matrix element-parton shower 
matching. 

We start with a discussion of the results of the ADD theory. From Fig. 5.11, one 
can immediately see that the shapes obtained with NLO and kj-MLM matching 
agree quite well for the pseudo-rapidity of the leading jet and the graviton and for 
the transverse missing pj . For the pj distributions of the first jet, the agreement is 
satisfactory for large pj > 800 Gev region, while for pj ~ 500 Gev, NLO curves 
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Fig. 5.12 NLO/matching comparison for the ADD model at the Tevatron (^/s = 1.96 Tev) with 
A = 1 Tev for S = 2, 4, 6. The matched curve is normalized by the NLO results and the normal- 
ization factors can be found in Table5.6. Figure extracted from [3] 



quickly decrease and become unreliable. That is due to the cut on the minimum pr, 

lstiet 

which leads to inconsistency on p T lower bound at LO and NLO calculations, 
while the matching technique gives reliable results. The same argument also applies 
for Ht distributions when Hj ~ 500 Gev. However, at large Hj > 1500 Gev, the 
kj-MLM matching tends to give harder distributions, because the contributions from 
the matrix elements of G + 3 partons are also included, while this is not true for NLO. 
Figure 5. 12 presents similar results for the Tevatron. Because the energy scale of the 
Tevatron is much smaller, the agreement is in general much better.The reason is that 
the contribution from more hard jet emission included with matching does not play 
an important role for such low center-of-mass energy. 
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Fig. 5.13 Graviton mass distribution at the LHC (left) and Tevatron (right) for the ADD model 
with S = 2, 4, 6, respectively. The NLO results are shown by the red, green and purple bands, and 
the matched ones are shown by the black curves. Figure extracted from [3] 

Finally, recall that the infinite sum of KK gravitons in the ADD model has been 
considered as an extra integral over the mass states. The mass of the experimentally 
expected graviton would therefore be given by a distribution instead of one specific 
value. In Fig. 5.13 we show the results for the graviton mass distribution. If the mass 
of the graviton would be an observable, one could easily constrain the number of extra 
dimensions using these distributions. That is because the average for the graviton 
mass depends on the number of extra dimensions. The larger it is the mass, more 
extra dimensions are predicted. 

Figure 5.14 presents results for the MGM model at the LHC. Similar distributions for 
the Tevatron are out of reach for M(/j, t ) — 1, 2 Tev (and excluded for smaller values), 
therefore the corresponding results are not shown. We apply here the running of the 
Planck mass cut: H™ n > M(/i»). One can see a turning point in the distributions 
of the leading or second jet pj and transverse missing pr due the cut. The NLO 
and matched results agree quite well on shapes for distributions of pseudo-rapidity 
of the leading jet and graviton, transverse missing pj and Hj of the jets. For the 
second jet pj distribution, fc^-MLM matching gives a considerably different shape 
when compared to the NLO calculation, especially for M(/i*) = 2 Tev case. For the 
leading jet distribution, the NLO and matched results for the LHC agree well for 
large pj > M(/i*)/2, while for pj < M(/x*)/2, the NLO curves are unreliable, 

Istiet 

again due to the inconsistency on p T lower bound at LO and NLO calculations. 
Finally, Figs. 5.15 and 5.16 present results for the RS model at the LHC and Tevatron. 
This analysis aims at a fare comparison between all classes of models. In reality, RS 
graviton is expected to decay and graviton distributions would have to be re-drawn 
as a functions of its decay products. However, theoretically, distributions of the 
graviton should be similar to the ones obtained by the system formed by its decay 
products. Therefore, we can use all distributions (displayed in Figs. 5.15 and 5.16) 
in a comparative analysis against ADD and MGM theories. 

Here, the agreement between the NLO and matched results is again much better 
at the Tevatron than the LHC for the same reasons pointed for the ADD theory. 
Also, kr-MLM matching tends to give harder Hj distributions especially for large 
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Fig. 5.14 NLO/matching comparison for the MGM model at the LHC. The NLO results are given 
by the red and green bands for M(/i*) = 1,2 Tev, respectively. The matched results are given by 
the black curves with the same parameters as for the NLO ones. The dominant Z — > vv background 
is also shown as a reference by the blue curve on the p, distribution. The k j-MLM matched curves 
are normalized by the NLO results, and the normalization factors can be found in Table 5.5. Figure 
extracted from [3] 
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Fig. 5.15 NLO/matching comparison for the RS model at the LHC with A = 3 Tev. The NLO 
results are given by the red and green bands for in q = 1 Tev and 1 00 Gev respectively. The matched 
results are given by the black curves with the same parameters as for the NLO ones. Matched curves 
are normalized by the NLO results, and the normalization factors can be found in Table 5.5. Figure 
extracted from [3] 



Hj > 1500 Gev due to the fact one extra jet is considered for matching when 
compared to NLO calculation. 

To conclude, allow us to perform one final analysis. In multi-jet final states, observ- 
ables exist that can be subject to large NLO QCD corrections at TeV scales [24]. 
The Hj or the Pj of the leading jet distributions are examples of such observables. 
They are sensitive to the opening of new channels with new kinematic topologies at 
NLO, which are enhanced by parton distribution effects or by special kinematic con- 
figurations (leading to large logarithm terms). One could imagine that higher order 
corrections would be dominant and for this reason not trust the perturbative series. 
In fact, one can easily show that this is not the case for more inclusive quantities. 
This behaviour is due to the choice of the specific observable. In any case, the effects 
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Fig. 5.16 NLO/matching comparison for the RS Model at the Tevatron with A = 3 Tev, for 
ntG = 100 Tev and 1 Tev, respectively. The matched curve is normalized by the NLO results and 
the normalization factors can be found in Table 5.5. Figure extracted from [3] 



in some specific region of the phase space can be large and need to be accounted 
for. Here, we show that samples with multiple jets in the final-state with ij-MLM 
matching provide a reliable description also of this kind of effects. 
In Fig. 5. 17 we display some illustrative plots. For simplicity, we consider the RS 
model at the LHC with A = 3 Tev and rriQ — 100 Gev, with jet cuts |?7j e t| < 4.5 
and pj l > 50 Gev, and graviton cut pj > 50 GeV. Comparing the LO and NLO 
curves, it is possible to notice the appearance of very large K-factors, especially at 
large pj , leading jet pj and Hj. As a result, we can see that the matched samples 
describe extremely well the NLO shapes. 

Let us now summarize the main aspects and findings within the current study. We 
have presented a detailed comparison between predictions from NLO QCD calcula- 
tions and the tree level matrix-element matching with parton showers in the context 
of three representative theories featuring a graviton: ADD, RS and MGM. This study 
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Fig. 5.17 NLO/matching comparison for the RS Model at the LHC with nig = 100 Gev and 
A = 3 Tev with jet cuts |r/j et | < 4.5 and j?j > 50 Gev, and graviton cut > 50GeV. Here we 
show the appearance of very large K-factors explicitly observed in some observables such as the 
p^ et and Hj distributions. Figure extracted from [3] 

validates and motivates the use of a fully automatized event generator such as Mad- 
Graph/MadEvent including spin-2 particles and multi-parton jet matching for LHC 
simulations. 

In ADD and MGM theories, the graviton appears as missing energy while in RS 
model it may be observed through its decay products. The fully inclusive samples 
have been normalized to the corresponding NLO inclusive total cross sections for 
producing G+X, which has been presented in details. Subsequently, key distributions 
for the most relevant observables involving one or more jets have been compared to 
those obtained by the NLO calculation G + jet + X. We have observed a complete 
agreement in shape, as well as in normalization, between the distributions obtained 
with NLO and kr-MLM matching calculations at the LHC. At Tevatron, values 
obtained for the normalization factors are very different from K-factors. This could 
be explained by the distinct PDFs used for NLO and MLM-matching calculations 
in addition to the scale definitions employed for each analysis. Finally, we note that 
observables which are sensitive to the number of jets, such as the Hj are better 
described by the matched sample than a fixed-order calculation. 
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Chapter 6 

Gravitino Production at the LHC 



In Chap. 2 we have presented several reasons to expect new physics at the TeV scale. 
One of these is the hierarchy problem, explained in Sect. 2.4.1. Supersymmetric 
theories can naturally solve the hierarchy problem, and at the same time, shed some 
light on different SM shortcomings, providing for example candidates suitable for 
being a dark matter particle. We start the chapter with a simple description of how 
quadratic divergences arise in the SM, and the possible solution for the hierarchy 
problem provided by Supersymmetric theories. We subsequently characterize such 
theories and describe some of its main branches. Finally, phenomenological analyses 
for the production of specific super-particles are exhibited. 

We start with a simple description of how quadratic divergences arise in the SM, and 
the possible solution for the hierarchy problem provided by Supersymmetric theories. 
We subsequently characterize such theories and describe some of its main branches. 
Finally, phenomenological analyses for the production of specific super-particles are 
exhibited. 



6.1 Supersymmetry and the Hierarchy Problem 

In the SM, quadratic divergences emerge when corrections to the Higgs propagator 
are calculated. More precisely, in Sect. 2.4.1 we have demonstrated such calculation 
for a toy-model containing one fermion and one scalar particle. Recall that the cor- 
rection of the mass of a scalar particle from fermionic loops is given by Eq. (2.32). 
For a theory with N f identical fermions, we have: 



9 X f N f 

5m] = ) — 



A — 6m i log 



A 



2m , 



O 



(1/A») 



(6.1) 
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where 5m 2 is the correction for the mass of the scalar particle. In the SM, the Higgs 
receives such corrections from fermionic particles in the loop. Notice that these 
corrections are quadratically divergent, causing the hierarchy problem. 
One could think, however, that new and not-yet-observed particles exist, and must 
be taken into account when corrections to the Higgs mass are computed. 
Assume a SM-like theory, where fermions are coupled to the scalar field <j>: 

C f = (6-2) 

The symmetry is broken and we take <p — (v + h)/*/2, being h the Higgs boson. 
Suppose for instance that a number N s of new scalar fields 0, exist, where i = 1 , 
2, n. They couple to the Higgs field following the lagrangian [1-4]: 

n 

£=X \d^ 2 -m\\^ 2 -\ s \h\ 2 m 2 -2v\ s h\^i\ 2 - (6.3) 

Notice that for simplicity, the coupling \ s is the same for new scalar fields. 
When computing corrections to the Higgs from the loops of new particles, as illus- 
trated on Fig. 6.1, again divergences are encountered. Contributions from the quartic 
coupling, which are exhibited in Fig. 6. la, are: 

A 2 - J^Zml log ( A) J + o ( 1/A 2) . (6 . 4) 

Combining Eqs. (6.1) and (6.4) we already observe that if we chose Nf — 2N S and 
— \ s = A 2 , the quadratic term in A cancels. The sign difference between corrections 
computed from the fermionic and scalar fields on the loop enable such cancellation. 
Contributions from the diagram displayed in Fig. 6. lb, can be written as 

p iog (£)] +o ( i/A! )- <65) 

Combining all contributions, we observe a complete cancellation of the divergences 
if we appropriately choose the coupling constants and the masses for the new scalars 
and SM fermions to be the same: 

Fig. 6.1 Illustration of loop 
corrections from new scalars 
in the calculation of the Higgs 
propagator 



X S N S 
' 16tt 2 



(b) 



\ 2 ,v 2 N s 
16vr 2 
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(m 2 f - m 2 s ) log (Jpj + 3 m) log + O (l/A 2 ) , (6.6) 

where we have assumed Nf — 2N S ,X 2 — —X s andm^ = m s . Therefore, if we have 
a theory with two scalar particles associated with each fermion existing in the SM, no 
divergences are encountered. In other words, the problem of quadratic divergences 
of the SM is solved if the theory becomes Supersymmetric. 

In a Supersymmetric theory, new scalar fields are associated with SM fermionic 
fields, and new fermionic fields are associated with the existing (SM) bosonic fields. 
In the next section we present an introduction to the Minimal Supersymmetric Stan- 
dard Model (MSSM) and specify the particle spectrum of the theory. Here, we have 
briefly demonstrated how the addition of new particles could generate the cancella- 
tion of divergences that appear on the calculation of corrections for the mass of the 
Higgs boson. 



6.2 A Brief Introduction to the MSSM 

Supersymmetry, or SUSY for short, is a symmetry that transforms bosonic states 
into fermionic states and vice-versa. It is therefore an operator Q which performs 
the transformations: 

Q\Boson) = \Fermiori) , Q\Fermion) — \ Boson). (6.7) 

The operator Q must be an anti-commuting spinor. For generating a consistent SUSY 
theory with chiral fields (necessary due to the existence of SM fermions), Q and its 
hermitian conjugate Q} must satisfy the following algebra [5]: 

{Q, C t }aP /1 ; {Q, Q] = {Q\ Q 1 } = 0; Q] = [P", g f ] = 0. (6.8) 

is the four-momentum generator of space-time translations. 
In the SM, we identify elementary particles with irreducible representations of the 
Poincare group. The irreducible representations of the SUSY algebra are not particles, 
but supermultiplets that contain both the SM corresponding fields and its super- 
partners. Therefore, field members of the same supermultiplet should have equal 
masses and be in the same representation of the gauge group. 



6.2.1 Particle Content of the MSSM 

The Minimal Supersymmetric Standard Model, or MSSM, retains this name because 
it respects the same gauge symmetries of the SM: SU (3) x SU(2) x U(l). The 



i A f N f 
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Table 6.1 The particle content of a minimal supersymmetric theory 



Superfield 



SU(3) 



SU(2) 


U(l) 


Particle content 
SM-like 


super-partner 


2 


i 

3 


u) 


/ u,\ 

\t) 

\ a L / 


1 
1 


4 
3 


Ur 


U R 


1 


2 
3 


d-R 




2 


-1 


CO 


ID 


1 


2 




ri* 
"R 


1 





a" 
y 


a" 


3 





W i 


w' 


1 





h 


b 


2 


+1 




(1) 


2 


-1 




(I) 



£ 

G" 

VI/' 

B 
H u 



In the first block, we present the chiral supermultiplets of the MSSM. Super-partner of the SM 
fermions are named sfermions (i.e., sup, sdown, selectron, etc). The first generation of SM fermions 
are shown, while the existence of the second and third generations are implicit. In the second block, 
we present vector multiplets with their particle content: the SM gauge bosons and their corresponding 
super-partners name gluino, wino and bino. Finally, in the third block the Higgs content of MSSM 
is exhibited 



particles necessary to construct such theory are the same as these of the SM, in 
addition to their super-partners. 

A summary of the supermultiplets with symmetries explicitly displayed and their cor- 
responding particle content is exhibited on Table 6.1. Chiral supermultiplets are only 
specified for the first generation. The existence of the second and third generations 
is implicit. 

As a consequence of the SUSY transformation, SM fermionic particles should have 
scalar super-partners (both included in the chiral supermultiplet), and SM bosonic 
particles should have fermionic super-partners (included in the vector supermulti- 
plet). Furthermore, super-partners should have the same mass and quantum numbers 
as their corresponding SM particle. For example, the super-partner of the electron, the 
selectron, should have a mass of approximately 0.5 1 MeV and have the same electric 
charge as the electron. If it existed, it should have already been detected. Because 
as of today none of these super-partners have been experimentally observed, if the 
world is supersymmetric, SUSY must be broken. 

In the SM, fermionic masses are generated from yukawa terms as it is exhibited in 
Eq. (2.28). However, the term — A„ QiffruR is not allowed in SUSY theories, since the 
superpotential should only involve superfields with left-handed fermion components 
and not their conjugates in order to prevent from chiral anomalies [3, 6]. To create 



To obtain chiral fermions in SUSY theories, right-handed fermion fields are defined to 
be the complex-conjugate of the left-handed fermion fields (XjR = X'i)- Similarly, their scalar 
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such term, which is relevant for the generation of the mass of part of the SM fermions, 
another Higgs doublet is necessary to be added to the supersymmetric lagrangian. 
Therefore, the MSSM theory has two Higgs doublets: H u and H d . The most general 
super-potential can be described by [5]: 

W = fly u QH U - % d QH d - % e LH d + fxH u H d , (6.9) 

where all the gauge and family indices are suppressed. Other gauge invariant terms 
could be part of the super-potential. However, these would violate lepton and/or 
baryon number conservation. A new symmetry, R-parity, is therefore imposed in the 
MSSM: 

P R = (_l)3(B-i)+2* f (6 10) 

where B and L are the baryon and lepton numbers, s is the spin of the particle and 
Pr is the quantum conserved quantity. Terms in the lagrangian are only allowed if 
the product of the Pr of all the fields (within the specific term) is + 1 . One can check 
that this is always true for each term of Eq. (6.9). Notice all the SM particles and the 
Higgs bosons have Pr — + 1 , while all the super-partners have Pr — — 1 . 
The consequences of the introduction of R-parity are several. For instance, since 
every vertex of the theory necessarily has an even number of super-particles, they 
will always be produced in pairs. Therefore, the lightest super-particle (LSP) should 
be stable and is thus a good candidate for dark matter. As one can see, through the 
imposition of the R-parity symmetry, the MSSM provides a natural solution for the 
dark matter problem of the SM. 

From Eq. (6.9), it is possible to write the classical scalar potential for the Higgs [5]: 
V = (|/i| 2 + m 2 Hu ) (\Hy- + |//+| 2 ) + (| M | 2 + m 2 H J(\HX + \H~\ 2 ) 
+ (|H t ?l 2 + \H+\ 2 - \HX ~ \H d -\ 2 f 

+ ^\H+H°* + H^H-*\ 2 + [b (H+H- - H°H°) + ex] . (6.11) 

Notice we have two couplings g and g', and the potential should give raise to v u and 
v d , the vev of both H„ and H d . 

One can count the degrees-of-freedom of the theory. Before the symmetry breaking, 
the theory contains two complex Higgs fields, summing eight degrees-of-freedom. 
When the symmetry is broken, SM gauge fields absorb three degrees-of-freedom 
to become massive. Five degrees are left to "new fields" that arise from the Higgs 



(Footnote 1 continued) 

super-partners are defined to be their complex-conjugate (ij>* and 4>' ). Since the super-potential only 

involves tf>' , one Higgs doublet will not generate mass for all the particle spectrum related to the 
SM. 
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Table 6.2 The particle content (mass eigenstates) of the MSSM [5] explicit for the SUSY and the 
Higgs sectors 



Names 


Spin 


Mass eigenstates 


Higgs bosons 





h°, H°, A , H ± 


S quarks 





ur, ul, <1r, di 






SR, SL, Cr, Cl 






t\, ?2, b\, i>2 


Sleptons 





eR, e L , v e 






flR, JlL, 






n, h, v T 


Neutralinos 


1/2 


~xl xl xl ~xl 


Charginos 


1/2 


~xtxi 


Gluino 


1/2 


9 



Notice the gauge eigenstates of the third generation of quarks and leptons are different from the 
mass eigenstate, therefore these are represented with numerical indices 



boson: two charged Higgs H ± ; a CP-odd neutral Higgs A ; and two CP-even neutral 
Higgs h and H. 

As a result from the SUSY symmetry breaking, higgsinos and gauginos mix giving 
rise to: (i) four supersymmetric neutral mass eigenstates, the neutralinos; and (ii) two 
charged mass eigenstates, the charginos. A summary of the particle content (mass 
eigenstates) is presented in Table 6.2. 



6.3 Supergravity: The Goldstino and the Gravitino 
6.3.1 Gauge Mediated Supersymmetry Breaking 

As of today, there is no dynamical way to break SUSY which is completely adequate. 
Supersymmetry is usually broken by hand, and several mechanisms exist to enable 
spontaneous breaking. Here we briefly discuss the so-called gauge mediated SUSY 
breaking (GMSB) scenario. 

In these types of models, a small extension to MSSM is necessary. It is assumed 
that SUSY breaking occurs in a hidden sector, also known as the secluded sec- 
tor. The SM particles and their super-partners are fixed in the visible sector. The 
communication between both sectors are made through the particles contained in 
the messenger sector, usually being one singlet field in addition to fields with 
SU (3)l x SU (2)i x U (1)y quantum numbers. Soft SUSY breaking terms 2 arise 
from loop interactions involving the messenger particles. This is convenient because 
it provides an automatic and natural suppression of terms with flavour change neutral 
currents and CP- violation, which are strongly constrained by experiments. 



2 Soft breaking terms are these added to the MSSM lagrangian (in order to enable spontaneous 
symmetry breaking) that do not create new divergent contributions to the mass of the scalar particles 
such as the Higgs. 
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In the MSSM the supersymmetry breaking is introduced explicitly. Goldstone bosons 
are thus generated, and are eaten to give mass to bosons. For global SUSY, the broken 
generator is fermionic, so the Goldstone boson in this case should have spin- 1/2 and 
is usually named goldstino. 

The goldstino is the eigenstate of the fermion mass matrix in which the eigenvalue is 
zero. The derivation of the goldstino interacting lagrangian depends on the supercur- 
rent conservation and it is independent from the details of how the supersymmetry 
breaking is communicated to the visible sector [5, 7, 8]. We shall not discuss more 
theoretical aspects of the theory, it is enough to know that the appearance of the gold- 
stino is provided by the breaking of a global symmetry. If SUSY is a local symmetry, 
the gravitino will be the field associated to the invariance of such symmetry. 



6.3.2 The Super-Higgs Mechanism 

The SUSY theory invariant under local transformations is named Supergravity 
[9-12]. In this type of theories, the gravitino has spin 3/2 and is the super-partner of 
the graviton. 

In order to see explicitly how local symmetry necessarily implies gravity, let us 
consider the lagrangian of a scalar field <f> and its fermionic super-partner \ [13]: 

£ = _(^*)(S^)-ix7^x- (6-12) 

The lagrangian is invariant under global transformations: 

0^0 + e X , (6.13) 
X X-i^edfifa 

but is not invariant under local transformations e —> e(x). To maintain the invariance, 
a field has to be introduced (similarly to what happens in the SM, where gauge fields 
An are generated from the invariance of a local symmetry). An extra term Ck must 
thus be added to the lagrangian: 

c ' = j^ K S'%> (6 - 14) 

where K® = —d^tjfx 01 ~~ i/^X (^"Ja d v (f>*, and ip is a Majorana fermion with 
spin-3 /2 that transforms like 

V a ^¥a + M p ,d^ a . (6.15) 

The problem is that the sum of both lagrangians C + Ck is not yet invariant. It has an 
extra contribution which is proportional to the energy-momentum tensor. However, 
this contribution can be cancelled if we add the term C g : 
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(6.16) 



The transformation of the tensor field is 



1 - 



9nv ~ * 9p,v + 



(6.17) 



which demonstrates that any locally SUSY theory must include gravity. For the 
MSSM, we have shown that chiral supermultiplets contain fermions and their super 
partners, and vector supermultiplets contain gauge bosons and gauginos. When super- 
gravity is considered, the theory also has the gravity supermultiplet, which contains 
the graviton and the gravitino: {g^ v , ip). 

Before the supersymmetry breaking, the gravitino is massless. To become massive, 
it absorbs the goldstino and obtain spin ±1/2 components. This is the correspond- 
ing to the Higgs mechanism of the SM, and for this reason is named super-Higgs 
mechanism. The mass of the gravitino depends on how supersymmetry is broken. 
Usually, the mass of the gravitino can be estimated following from dimensional analy- 
sis [5]. Let (F) be the vev of scalar fields responsible for supersymmetry breaking. 
We thus expect the mass of the gravitino to be: 



since must vanish when SUSY is restored, {F) —*■ 0. 

Heretofore, we have avoided to introduce explicitly the SUSY lagrangian because of 
its complexity. However, since gravitino interactions will be necessary to subsequent 
phenomenological analyses, we briefly introduce the interaction lagrangian of the 
gravitino, and show that the dimensional analysis for the mass of the gravitino is 
adequate. We further demonstrate how to obtain the Feynman rules related to the 
gravitino and verify the gravitino/goldstino correspondence. This summary is based 
on Refs. [13-17]. 

Let us start writing the terms in the supergravity lagrangian where the mass term for 
the gravitino arises. The supergravity lagrangian depends on arbitrary real functions 
of scalar fields cf>* and <j>; with i, j = 1, n. It is called the Kahler function: 



where K is the Kahler potential and W is the super-potential. 3 The part of the 
lagrangian that include terms for the gravitino is given by: 



(F) 
M p i 



(6.18) 



G(<^ X ) = W,x) + ln|W(</>)| 2 , 



(6.19) 




1 



^HGij+GiGj)^^ 
(6.20) 



V2 



3 See Ref. [12] for more details. 
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where 

dG d 2 G 

Gi ~ ~£T< G U = a j. xj. ' and 9 = det{g^, v ). (6.21) 
o<Pi ocpjOcpj 

The M p i is momentarily set to 1 to simplify the format of Eq. (6.20). The second 
term specifically shows the mix of the gravitino and the goldstino r), 

ri = Grf. (6.22) 

We can redefine the fields, where the gravitino is now ip': 

^ = V - - —e GI2 d^. (6.23) 

Equation (6.20) is thus: 

C = -« G/ X*^ " ( Gy + ^GiGjj x'x', (6.24) 



V-.9 

from which we obtain the mass of the gravitino: e G / 2 . Recovering the M p i depen- 
dency, the mass of the gravitino 4 is: 

m^ = e G/2 M„, = e K/2 . (6.25) 
M p i 

One can observe from Eq. (6.24) that the gravitino ip' absorbs the goldstino rj and 
obtain mass. We can thus conclude that the dimensional analysis giving rise to 
Eq. (6.18) is adequate. Once auxiliary fields F' 

F i = e G/2 



_ j^Gj* (6.26) 

obtain a vevs, the symmetry is broken, and gravitinos obtain mass. 
Now we have seen that spontaneous supersymmetry breaking can occur. We shall 
not discuss in details the mechanism to generate it, rather, we are interested in how 
the theory presents (in the visible sector) after symmetry breaking. 



6.3.3 The Effective Theory for Light Gravitinos 

From this moment on, we use the supermultiplets in the four-components notation. 
We follow the notation presented in Ref. [14]. The four-component spinors can be 
constructed from two-component spinors. For the gravitino gaugino A and chiral 



4 From now on we drop the prime index: ij/ — ► ip, 
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fermions xl,r we have: 

=(!;)' (6 - 27) 

A (M) = f ^AO = ( A a) , (6.28) 
Xf = (j). xi D) = (0x), (6.29) 

where stands for Majorana and for Dirac. The fields satisfy either the Majorana 
or Dirac conditions depending on the given index. From now on we drop the and 
(£> ' indices to simplify the notation. 

The term in the supergravity lagrangian for the massive free gravitino, considering 
the mass term given by Eq. (6.25), is given by [14]: 

= - l -^p^^ l5lljdp ^ a - ^L* M [y< , y] V„. (6.30) 

The polarization tensor for a gravitino with momentum p, is [18]: 

A 

where the sum over the four gravitino helicities (A = ±1/2, ±3/2) is performed. 
Because we will be interested on the phenomenology of gravitino production at the 
LHC, the energy scale of gravitino production will be much larger than the gravitino 
mass. The polarization tensor can thus be simplified as [18]: 

IVO?) = - Pg, w + \p (6.32) 
3 m% 

where the first term corresponds to the A = ±3 /2 helicity states and the second to 

A = ± 1 /2. We will discuss this equation in details further ahead. 

The interaction terms of the gravitino are defined through the supercurrent by: 

C = — J^VnJ^ + h.c. (6.33) 
2M p i 



The supercurrent enables interactions of the gravitino with all the other fields. 
Although the supergravity lagrangian contains many interaction terms for the grav- 



6.3 Supergravity: The Goldstino and the Gravitino 



109 



itino, the only relevant terms for our phenomenological studies will be the following: 

-f]-fX tt FS*- (6-34) 



8M 



Pi 



Where D v 4>' — d v <\> 1 + ig A" <j) } , T,j is the generator of the gauged Lie group, <fi 
is the scalar field present in the chiral multiplet, xl/r are chiral fermions and A are 
gauginos. 

The Feynman rules for gravitino interactions arise from the lagrangian (6.34) and are 
shown in Fig. 6.2. A comprehensive derivation of the lagrangian and Feynman rules 
related to the gravitino and some discussions can be found for example in Refs. [12, 
14, 16, 19]. 

Gravitinos obtain mass from (absorbing) goldstinos after SUSY breaking. The longi- 
tudinal modes of gravitinos arise from the two helicity state components of the gold- 
stinos. As a result, gravitinos will have in total four helicity states. Equation (6.31) 
presents the polarization tensor of an external gravitino. For the case where */s 3> 
Eq. (6.32) shows that the polarization tensor can be explicitly separated in terms of 
A = ±3/2andA = ±1/2 [14, 16, 18]. From this equation, one can realize that there 
is a correspondence between the gravitino and goldstino lagrangians. The interacting 
lagrangian of the goldstino with matter can be obtained by the following substitution 
[15, 20, 21]: 




5 Other terms are irrelevant for the present analysis because we will consider processes at the LHC 
energy scale, which is much smaller than the scale at which quantum gravity corrections become 
important (M p i). 
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Fig. 6.3 Feynman rules for the goldstino interactions 
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(6.35) 



Where is the spin- 1/2 goldstino field. Considering that we are interested in processes 
with a light external gravitino, it is possible to use the relation 0ij — 0. Furthermore, 
because in supergravity high energy behaviours cancel out the total amplitude of a 
process (with one external gravitino), the derivatives on 0, % ar *d A can be replaced by 
the masses of the corresponding fields [14]. One can thus obtain the final interacting 
lagrangian of the goldstino: 



2 2 

nil — m; 
r -i t *_ 

V3/ 



im* M p i 



VXR> 



8 *j6myMpi 



fj[>f,<f\ \ a F"+h.c 



(6.36) 



and derive the Feynman rules, illustrated on Fig. 6.3. When the center-of-mass energy 
is higher than the mass of the gravitino, the gravitino is mainly composed by ±1/2 
helicity states. Therefore, at the LHC, it is sufficient to consider the Feynman rules 
of goldstinos. 

Our goal is to continue to analyse processes that could be observed at the LHC 
through the missing energy signature in combination with jets. Several such stud- 
ies employing SUSY theories have already been extensively performed. However, 
searches for the gravitino in the final state have not been largely explored due to 
limitations of simulation tools. A recent implementation of gravitinos and goldsti- 
nos in MadGraph/MadEvent [17, 22] has changed the situation. In collaboration 
with the authors we have been able to perform a phenomenological study on grav- 
itino/goldstino 6 production with multiple jets in the final state, which we present in 
what follows. 



6 From this moment on, we will refer to it as the phenomenological study of the production of 
gravitinos, which we have shown to be the same as the production of goldstinos in the low energy 
limit. 
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6.4 The Phenomenology of Gravitino Emission 

We have seen that in the context of supergravity theories, gravitinos are massive 
fields, super-partner of the graviton. In some variations of such theories, these are 
super-light particles, however not exactly massless [23, 24]. 

Similarly to what has been presented in Chap. 5, we are interested in the missing 
energy signature in combination with jets. The gravitino, being the LSP, is expected 
to be invisible. In addition, because SUSY particles are always produced in pairs 7 , 
we consider the production of two gravitinos. Therefore, we now focus mainly on the 
process p p — ► G G+«-partons, where now G represents the gravitino 8 . For simplic- 
ity yet being exhaustive, we will consider n — 3 as the maximum number of partons. 
However, super-light gravitino production can become relevant, since the cross 
section is proportional to the inverse of the gravitino mass. From Fig. 6.3, we observe 
the interaction vertices of the gravitino at the LHC energy scale. Different combi- 
nations can be considered to form processes in which the gravitino is one of the 
products of a proton-proton collision. Here, we concentrate on scenarios where the 
gluino is the NLSP. Sub-processes that contribute most significantly to our analysis 
are illustrated in Fig. 6.4. They are summarized as: 

We consider a simplified model with the Feynman rules explicitly defined in Fig. 6.3. 
The particle content is reduced, focused for the present study. Only the SM quarks and 
gluons coupled to the gravitino/goldstino and the gluino will contribute significantly. 
Cross sections related to graviton production at the LHC are suppressed by the Planck 
mass. For this reason, graviton production does not contribute for the present study. 

• Gluino- gravitino associated production: 



The first sub-process will be often referred hereby as "gluino-gravitino associated 
production", and the second as "double gluino resonance". In principle, the zero-jet 
process should also exist, in which two gravitinos are directly produced. However, 
here we specifically focus on processes with at least one gluino, which decays pro- 
ducing the gravitino. 



The decay width of a gluino (g) to a gravitino (G) and a gluon (g) is given 
by [14, 17]: 



pp — >• Gg^f G G g + extra partons; 



(6.37) 



• Double-gluino resonance: 



pp — ► g g — s- G g G g + extra partons. 



(6.38) 




in~ 

9 



(6.39) 



48nM pl m 2 £ 



7 We consider R-parity conserving theories. 

8 Within the present study we consider processes in which at least one gravitino decays from a 
gluino, as illustrated in Fig. 6.4. 
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1-parton 2-partons 3-partons 




Fig. 6.4 Example of processes contributing to pp — ► GG+ jets. The three-lines initial state 
represents a proton. The double-straight-line on the final state represents a gravitino. The curly 
line is a gluon, and the curly-line crossed by a straight line represents the gluino. The filled circle 
symbolizes that all the processes that can possibly generate such final state are included in the 
illustration of the diagram. Notice that the diagrams are separated in columns characterized by the 
number of final state partons. The diagonals represent matched diagrams. The second diagram of 
the last column is not matched in order for the highest multiplicity to be the same for both processes 



Notice that the decay width depends both of the gravitino and the gluino mass. For a 
fixed value of the gravitino mass, the decay width rapidly increases with the increase 
of the the gluino mass. 

Let us explore first the cross sections of each sub-process presented above. However, 
because we consider Br(g — >■ Gg) = 1, it is enough to compute the 2^2 sub- 
processes that give rise to pp — > Gg and pp — > gg. 

We can now start by estimating the cross section for pp — > gG. One has to consider 
all the three possible initial states: pair of gluons, quark-anti-quark and quark-gluon 
(anti-quark-gluon). Each vertex including gravitino is inversely proportional to the 
gravitino mass. The result is thus inversely proportional to the square of the gravitino 
mass: 

*(pp -> gG) oc -L (6-40) 

m~ 
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Fig. 6.5 Cross section of pp — s> gG at the LHC separated for each contributing sub-process. 
To obtain the curves, we have employed MADGRAPH/MADEVENT with CTEQ6L1 PDF. The mass 
of the gravitino is fixed at 10~ 13 GeV, the mass of the squarks at 3TeV and fi = {ir = [ip = 
(mg +zn§)/2 ~ m-g/2 

which has several consequences for our analysis. Equation (6.40) also depends on the 
gluino mass. However, each sub-process has a different dependency on m~„, and for 
this reason it is not explicitly written. Figure 6.5 exhibits the cross section dependency 
on the gluino mass for a fixed value for the gravitino mass. An exact calculation of the 
amplitudes and cross section for each sub-process can be found for instance in [25]. 
Similarly, the cross section for producing a pair of gluinos can be calculated. As 
expected, it does not depend on the mass of the gravitino, since no interacting vertex 
with the gravitino exists. Figure 6.6, illustrates the dependency of the cross section 
on the gluino mass. 

The cross sections of both sub-processes have a particularity. If we fix the gluino mass, 
the cross section for the double-gluino resonance is determined. However, the cross 
section for gluino-gravitino associated production process varies depending also of 
the mass of the gravitino. For analysing processes with two gravitinos and jets in 
the final state, one has to take into account contributions from both sub-processes. 
The mass of the gravitino determines therefore which sub-process contribution is 
dominant. 

In principle, it would be necessary to divide the simulation of pp — > GG + n-jets 
in three stages: simulation of dominant gluino-gravitino production, simulation of 
dominant double-gluino resonance; and the simulation of equally contributing sub- 
processes. However, the matching technique automatically accounts for both sub- 
processes at the same time, independent of the region of the parameter space. This 
is one substantial advantage of using such technique. There is therefore no need for 
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Fig. 6.6 Total cross section of pp — > gg separately in gg -» gg, qq — > gg contributing sub- 
processes. To obtain the curve we have employed MadGraph/MadEvent with CTEQ6L1 PDF, 
fixed the mass of squarks at 3 TeV, the mass of the gravitino at 10~ 13 GeV and /i = /j/e = [ip = nig 



simulating each sub-process separately, and we can thus treat the whole parameter 
space at once. 

Our aim is therefore to perform a complete analysis on gravitino production in 
combination with jets, in which we make use of the matching technique presented on 
previous chapters. We know that apart from normalization issues, the employment of 
the matching technique can correctly predict the shape of relevant observables and 
can cover all regions of the parameter space. 



6.4.1 Parameter Definitions 

In the present analysis, the gluino mass is fixed at 800 GeV, value chosen based on 
recent experimental searches [26]. The gravitino mass is selected according to the 
three cases below: 

CaseA. a(pp — > gg) < <r(pp -*■ Gg) (6.41) 
CaseB. o(pp -*■ gg) » cr(pp — >• Gg) (6.42) 
CaseC. <j(pp — > gg) > cr{pp Gg) (6.43) 

We observe from Eq. (6.42) that the value for the gravitino mass in "case B" (m B ~ ) 

G 

is selected as to fix a(pp — > gG) & &(pp — > gg)- We define one heavier and one 
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lighter gravitino, as stated by Eqs.(6.41) and (6.43), referred here as "case A" and 
"case C". Their masses and cross sections are related to case B as: 

A 1 B - 1 

m% — ~ m a ~, cta(pp ->■ gG) » a(pp ->■ gg), 
° n « « z (6.44) 

m 6 =nm f' a c(pp ^ gG) ^ n 1 a(pp ^ gg), 

For our investigation we take « = 3, i.e., = 1 x 10~ 13 GeV, = 3 x 10~ 13 

& G G 

GeV and m c . = 9 x 10" 13 GeV. 

G 

The SUSY model with gravitinos has been implemented in the FeynRules pack- 
age [27], and we make use of such implementation [28]. It can be used by Mad- 
Graph/MadEvent version 5 thanks to the UFO and ALOHA packages [29, 30] 
which translate the Feynman rules obtained from FeynRules into a universal output 
that can be used by event generators, such as MadGraph/MadEvent. In particular, 
we have worked for the development of the Aloha package, which is detailed in the 
Appendix C. 

The complete analysis is performed for the LHC with *Js = 14 TeV and we employ 
the CTEQ6L1 PDF. No cut is used for the generation of events, except for a minimum 
jetp T : 

p T > 20 GeV. (6.45) 

For the matching technique, we employ the kj-MLM scheme, use 
MadGraph/MadEvent version 5 [31] for the matrix element simulation and 
Pythia for parton shower and hadronization. We choose Q match = 100 GeV, for jet 
clustering we use the anti-^r -algorithm with R — 0.5 [32], and impose a cut in the 
minimum pr of the jets: 

pr>50GeV. (6.46) 



6.4.2 Checks and Validations 

Before displaying the results for a specific analysis, it is crucial to validate the imple- 
mentation. Since we have two distinct sub-processes contributing separately for the 
final cross section, such checks are indispensable. 

The first check is made as a separate study. First, we consider both sub-processes 
with one extra jet, i.e. pp — > GG + n-jets, with n = 1, 2 jets for gluino-gravitino 
associated production and n — 2, 3 jets for double-gluino resonance. We employ 
the parameters presented above and make use of the matching technique. We plot the 
distributions for different relevant observables. To check the validity of our imple- 
mentation, besides the total matched simulation, we can separately simulate the 
contribution of each sub-process: a G ~ and o~ g ~g- We can check therefore, that the sum 
of both contributions agree with the overall matched distribution. 




Fig. 6.7 Differential cross section as a function of the Hj (left) and p T (right). Here we exhibit 
separately the contributions of both pp — » gg and pp -> sub-processes that are included in the 
evaluation of pp — »• GG j+ 1 extra jet. The sum of the contributions agree with matched sample. 
We employ a pj l > 20 GeV cut, and the results are for the LHC with = 14 TeV 



These plots are interesting to demonstrate the contribution of each sub-process for 
the three sets of parameters chosen. For example, Fig. 6.7 presents the differential 
cross section as a function of the Hj (left), and of the fir (right). We can observe 
that for case C, the double-gluino resonance dominates, as requested in Eq. (6.43), 
and the sum of the contributions is very similar to this specific sub-process. For case 
A the opposite happens: the sum of the contributions is mainly given by the gluino- 
gravitino production, as Eq. (6.41) show. Finally, for case B, we can check that both 
contributions are equally relevant. 

Different from NLO calculations, the matching technique only considers real- 
emission contributions. For high-energy collisions processes with initial/final state 
radiation become most likely and the matching technique can treat correctly these 
real contributions. However, virtual corrections are not considered, which are needed 
to cancel the infrared divergences that could arise in the calculation of real QCD- 
emission corrections. Therefore, with the employment of the matching technique, we 
can obtain the correct shape for the distributions of the differential cross section, but 
we cannot predict the K-factors needed to correctly rescale the final cross section: 

da(GG +jets) = da(gg + jets) + K-gda(gG +jets). (6.47) 

We thus assume Kgg ~ ~ 1 in the following studies. 



6.4.3 Results 



We now consider pp — > GG + 1,2, 3 -jets. The relevant diagrams are illustrated 
in Fig. 6.4 separately in columns for 1-jet, 2-jets and 3-jets contributions. These 
diagrams exemplify the main sub-processes existent. Notice that diagrams with no 
initial-state-radiation (ISR) are illustrated in the first line, with one ISR in the second 
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Fig. 6.8 Differential cross section as a function of the relevant observables. The signal is pp — >■ 
GG + rc-jets with n = 1, 2, 3. The mass of the gluino is 800 GeV, and the mass of the gravitino is 
1 x 10~ 13 GeV in case A, 3 x 10" 13 GeV in case B and 9 x 10~ 13 GeV in case C. The irreducible 
background, pp —> vv + re-jets with re = 0, 1,2, 3, is shown for reference. The simulation is 
performed using the kr-MLM matching technique at the LHC with = 14TeV and CTEQ6L1 
PDF 



line and two ISR in the third line. Therefore, diagrams that should be matched are 
represented in the diagonals. 

Because there are two distinct sub-processes relevant for the analysis, PYTHIA cannot 
automatically correctly set the highest jet multiplicity. The reason is that the leading 
order contribution for each sub-process has different multiplicity. It is necessary 
therefore to define the number of additional jets as an input of the analysis. We have 
set it to be one, which removes the explicit contributions of the second diagram of 
the last column of Fig. 6.4 that has two radiated jets. 

The irreducible background is given by the matched sample of pp — >• vv + n-jets, 
where n can be 0, 1,2 and 3. To simulate the background we have used a cut in the 
missing Ej, 

E T > 200 GeV, (6.48) 
in order to improve the statistics in the high energy region. 

Figure 6.8 collects the results obtained, in which we display the differential cross 
section as a function of: 



(i) the pt of the leading jet; 
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Fig. 6.9 In the left, the fraction of events (N events / N to tal) for different multiplicities in each sample 
is displayed. In the right, the final number of events expected as a function of the jet multiplicities is 
shown. The signal is pp — »• CG + rc-jets with n = 1, 2, 3. The mass of the gluino is 800 GeV, and 
the mass of the gravitino is 1 x 10~ 13 GeV in case A, 3 x 10 -13 GeV in case B and 9 x 10~ 13 GeV in 
case C. The irreducible background, pp ~* ^i> + n-jets with n = 0, 1, 2, 3, is shown for reference. 
The simulation is performed using the kj-MLM matching technique at the LHC with *Js = 14 
TeV, L = 100 ftr 4 and CTEQ6L1 PDF 



(ii) the pr of the second highest jet; 

(iii) the p T \ 

(iv) the H T . 

From Fig. 6.8 (i) and (ii), one can immediately see that the shapes for cases A, B and 
C are similar, and the curves differ mainly in overall magnitude. However, shapes 
for the pr and Hj distributions are slightly different in the high-energy region. 
Case A tends to give harder distributions than cases B and C. This is due to the most 
contributing sub-process of each sample. In case A, the dominant contribution comes 
from the gluino-gravitino associated production, and in case C from the double- 
gluino resonance. The $t is a combination of the pr of both gravitinos. In case A, 
one gravitino is directly produced and the other comes from the decay of the gluino. 
In case C, both gravitinos come from the decay of the gluino. A more detailed analysis 
related to the kinematics is performed further ahead, when we present the correlation 
between relevant observables. 

In Fig. 6.8 we also display the irreducible background curve as reference. The curves 
for cases A and B lie above the background distribution and could thus indicate 
signs of new physics at the LHC. However, for case C more complete analyses are 
necessary. 

One indication of new physics and of the dominant sub-process contribution can 
be given by the jet rates. Case A, gluino-gravitino associated production, has one 
hard-jet in addition to extra jets coming from QCD radiation. Case C, double-gluino 
resonance, has two hard-jets and contributions from radiated extra jets. Therefore, it 
is expected that the number of events with one jet is larger for case A than for case 
C, being case B intermediate. This is what we observe in the left graph of Fig. 6.9. 
Moreover, distinct behaviours are noticed for each sample. The background has the 
shape of decreasing stairs, where most of the events have 1-jet, while signal exposes 
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Fig. 6.10 Scatter plots for cases A, B and C. These illustrate the correlation between the pj of the 
leading jet and the pj . For case A, in which gravitino-gluino production is dominant, most of the 
points are distributed following a straight diagonal line. For case C however, where double-gluino 
resonance is dominant, points are spread into an ellipsoid region around the tfj = p r stJct diagonal 



a peak. The position of the peak varies for the signal samples. Case A has a similar 
number of events with 1 and 2-jets. The peak is at events with 2-jets for case B, and 
in case C the maximum is reached for events with 3-jets. If signals of new physics 
are observed at the LHC, one could therefore use the shape of jet multiplicities to 
differentiate between theories, and/or to constrain parameters. However, the right plot 
of Fig. 6.9 exhibits the expected number of events as a function of the jet multiplicities 
for the present analysis (in which we use L — 100 fb _1 ), showing the difficulty in 
observing any signal when accounting for the background. 

In order to curb the background, we investigate the correlation between observables. 
One interesting result is exhibited in Fig. 6.10. These are scatter plots with p T and 
the p'j distribution in both axes, in which each point represents one event. The three 
different cases considered are displayed. 

For case A, where gluino-gravitino production is dominant, points are concentrated 
around the main p 1 ^ 1 ^ 1 — p'j diagonal. This is expected since for this process, the 
hard-jet is the gluon decayed from the gluino. Extra jets in events with more than one 
jet are most likely to have been generated from QCD radiation, which results in jets 
with a reasonably low pj . In this case, the gravitino directly produced has a fairly 
high pj, larger than that of the gravitino decayed from the gluino. The tendency is 
therefore to obtain the hard-jet approximately back-to-back 9 with the missing energy, 
which results in events concentrated in the diagonal region as we observe from the 
left plot of Fig. 6.10. 

For case B, we have both sub-processes equally contributing. We can thus observe 
an intertwine behaviour between cases A and C. Figure 6. 10 shows the evolution 
of the contributions from gluino-gravitino dominant production to double-gluino 
resonance, and how it affects the pj of the particles involved. 
For case C, the main contribution comes from the double-gluino resonance, in which 
we have two hard-jets. Now, the probability is larger for both gluinos to be back- 
to-back, and the pj of the leading and the second jets are likely to be of the same 



Jets are defined to be back-to-back when they have equal momenta but opposite directions. 
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order of magnitude, but opposite directions. A similar argument can be used for 
both gravitinos produced in combination with each hard-jet. However, we cannot 
distinguish between the gravitinos, as the missing energy is experimentally obtained 
from visible particles. Since the missing pj is a combination of the pj of both 
gravitinos, in this case, there is a larger probability that they will (partially) cancel 
each other. The tendency is therefore to have events more spread around the diagonal, 
with the high p T and low pj region populated. That is what we observe in the 
third plot of Fig. 6. 10. 

To reduce the background we use scatter plots to search for the most populated region 
with signal events and scarcely any background events. Figure 6. 11 shows the back- 
ground superimposed to the signal — case B (equally contributing sub-processes). 
We can thus impose the following cut 

p T > 500 GeV, /4 stje ' > 500 GeV, (6.49) 

which is represented by the grey lines in Fig. 6.11. The differential cross section as a 
function of the relevant observables, accounting for the cuts given by Eq. (6.49), are 
collected in Fig. 6.12. 

From Fig. 6. 12 we observe that consequences for the observables related to the pj of 
the particles do not change considerably at the high pj region with the application 
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Fig. 6.11 Scatter plot of the signal pp -> GG +1,2, 3-jets (case B), in red, and the background 
pp — > vv + 0, 1, 2, 3-jets, in black. Each point represents one event. Here, we can easily see that 
the background is concentrated in the small Pj Je and jpj region. To reduce the background we 
have applied the cuts specified by Eq. (6.49), i.e. we only accept events in the region represented by 
the upper black box 




Fig. 6.12 Differential cross section as a function of the relevant observables after cuts: pj > 
500 GeV and p^ 6 ' > 500 GeV. The signal is pp — ► GG + n-jets wifhn = 1, 2, 3. The irreducible 
background, pp -> vv + «-jets with ;1 = 0, 1,2, 3, is shown for reference. The mass of the gluino 
is 800 GeV, and the mass of the gravitino is 1 x 10~ 13 GeV in case A, 3 x 10~ 13 GeV in case B and 
V x 10" 13 GeV in case C. The simulation is performed for the LHC with ~Js = 14 GeV 



of the cuts. One could possibly observe signals of new physics for cases A or B, 
however the curve for case C still runs together with the background. 
Nevertheless, Fig. 6.11 shows that applying the cuts specified in Eq. (6.49) should 
remove most of the background. We can thus plot the fraction of events observed after 
such cuts as a function of the jet rate. From the left plot of Fig. 6. 1 3 we observe that 
most of the background is removed, whereas a great part of the signal remains. The 
shape of the distributions change from before the imposition of cuts. The background 
distribution has now a peak for events with 2-jets, and the distribution for case A 
obtains the decreasing stair shape. Let us stress that the left plot of Fig. 6. 13 only 
confirms that a great part of the background is removed. However, it is the plot with 
number of events as a function of the jet multiplicity, in the right, that can show us 
the distributions fairly compared. We observe that case A is possibly be detectable 
at the LHC. For cases B and C, an excess at higher multiplicities (more than 4 jets) 
could be indication of new physics. 

To summarize, the main difficulty of the present analysis lies in the fact that two 
different sub-processes are substantially relevant. Normally, they would have to be 
considered separately and a careful analysis would be necessary to gather all the 
contributions into one final result. This process would be susceptible to unintentional 
errors due to several reasons, such as an under/overestimation of one of the sub- 
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Fig. 6.13 In the left, we show the ratio of events N e vents(,AC) I ^total(BC) f° r different multiplicities 
in each sample, where A 'events(AC) i s me number of events after the imposition of the cuts specified 
in Eq.6.49 and N, ola i(BC) is the total number of events before the application of any cut. In the 
right, the final number of events expected as a function of the jet multiplicity is displayed. The 
signal is pp — > GG + ;i-jets with n = 1, 2, 3. The irreducible background, pp vv + n-jets 
with n = 0, 1, 2, 3, is shown for reference. The mass of the gluino is 800GeV, and the mass of 
the gravitino is 1 x 10~ 13 GeV in case A, 3 x 10~ 13 GeV in case B and 9 x 10~ 13 GeV in case C. 
The simulation is performed for the LHC with = 14 GeV and L = 100 fb _1 



processes for example. A significant upshot of our technique is therefore the ability 
of treating both sub-processes at once, at any region of the parameter space. 
If signals of new physics emerge at the LHC as an excess of events with high jet 
multiplicities, more refined analyses will be necessary to identify the theory and 
constrain parameters. Essentially, a complete control of jet parameters are required, 
since jet rate distributions can dramatically change if one parameter is modified. 
Nevertheless, once jet parameters are dominated, and the mass of the gluino is fixed, 
one could use jet multiplicities to constrain the mass of the gravitino. We observe that 
different peaks in Fig. 6. 13 could allow us to set limits in the mass of the gravitino, 
and identify what sub-process is dominant. 

The mass of the gluino must be fixed in order for the analysis with the jet rates to be 
useful. At the LHC gluinos, being the NLSP, could be identified through different 
signatures. Within the strategy of the present study, gluinos should decay into one jet 
and one gravitino, which is the invisible particle. Therefore, gluinos could be iden- 
tified through analyses that use the distributions of the differential cross section as a 
function of the pj of the leading jet. Because of the kinematics of the processes being 
considered, the peak observed in the signal distributions should lie approximately at 
nig /2. This is what we observe in Fig. 6.8. For case C, where double-gluino resonance 

is dominant, the peak is clearly at p^ et = 400 GeV. For case A the determination of 
the peak is more difficult, because the dominant contribution is the gluino-gravitino 
associated production, but still possible. As we have stated, once the mass of the 
gluino is fixed, one can use the shape of jet rate distributions as an indication of the 
mass of the gravitino. 
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Chapter 7 
Epilogue 



7.1 Summary and Conclusions 

The SM describes very accurately the experimental data we have so far, however, 
chances are high that signs of new physics will appear soon at the LHC. Several 
BSM theories have been proposed as an alternative to the SM at high energies. If 
the theory extension could incorporate gravity in the quantum description of Nature, 
one of the greatest puzzles of particle physics would be unravelled. Theories with 
such attributes are therefore specially interesting. One can cite for example higher 
dimensional theories, such as large extra dimensions or warped extra dimensions. 
Another option are 4-dimensional theories with a separate sector that contains many 
particles that do not interact with the SM particles through the usual weak, strong 
and electromagnetic forces. One example is the massless graviton model presented 
in Chap. 3. There are also theory extensions that allow for extra symmetries, such as 
supersymmetry and supergravity, both introduced in Chap. 6. 
Given the large number of existing scenarios, it will certainly not be trivial to deter- 
mine the nature of new physics once it emerges at high energy hadron colliders. 
Simulations, and search strategies will thus play a key role in the determination of 
such theories, and moreover, in the treatment of new physics. Therefore, the neces- 
sity of general simulation tools that quickly predict experimental observables have 
become essential. Furthermore, the demand for new search strategies has substan- 
tially increased with the advance of computer technology and simulation tools aimed 
at the treatment of particle physics focused on theory extensions to the SM. 
The present work has been intended to accomplish a complete phenomenological 
analysis at hadron colliders of BSM theories that incorporates gravitational interac- 
tions. To this aim, the development of new simulation tools and search strategies has 
been necessary, and performed. We have carried out an analysis that covers all the 
steps to go from a theoretical framework to the experimental observables expected 
at hadron colliders. 

The format of this thesis has also followed the "theories — > experiments path". We 
have started with an brief introduction followed by an overview of the SM and some 
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of the open questions that we hope to be (partially) answered by the LHC. In the 
third chapter, we have introduced theory extensions to the SM that contains graviton 
particles. A rather theoretical description has been presented, in which we show the 
difference between theories that include massive gravitons to these with massless 
gravitons. In particular, it has been highlighted that the propagator of such particles 
has a unique behaviour when the limit from a massive to a massless graviton is 
considered. 

The connection to a more experimental point-of-view has started with the calculation 
of differential cross sections for graviton production at hadron colliders. We have 
considered different types of theories, containing both massive and massless gravi- 
tons. We have observed that the unique behaviour of the propagator of such theories, 
known as the van Dam-Veltman discontinuity, does not appear at cross section level. 
This is an interesting consequence, since the opposite would be expected. We cannot 
therefore use the discontinuity of the graviton propagator as hint for the determi- 
nation of the mass of the graviton particle if it is observed at the LHC or at future 
colliders, and more phenomenological analyses are necessary. 
The proximity to more experimental searches has increased in Chap. 4, when we have 
presented the LHC and defined observables. Very energetic particles are the conse- 
quences of high energy collisions, such as these happening at the LHC. When the 
center-of-mass energy increases, final-state full of jets become more likely. There- 
fore, it is a natural choice to study signatures that predict numerous particles as a 
result from the collision. 

For this reason, in Chap. 5, we have focused on graviton production in combination 
with jets. We have started with a simple phenomenological investigation of processes 
with one graviton and one jet at the LHC. We have used the Monte-Carlo event gen- 
erator MadGraph/MadEvent for the simulation of events, which included recent 
implemented helicity routines for correctly treating massive spin-2 particles. We 
have manually implemented the helicity routines needed by massless spin-2 parti- 
cles, and the MGM model within the event generator. With such implementations, 
MadGraph/MadEvent has become complete in terms of spin-2 particles, i.e. could 
automatically generate events for any type of spin-2 particles. Results for this par- 
ticular analysis have been obtained as a function of relevant observables at the LHC, 
and compared to the expected SM background. 

However, at the LHC larger number of jets in the final state are expected. Therefore, 
despite the significant results obtained with the mono-jet analysis, more realistic 
investigations needed to be performed. Because the emission of QCD radiation needs 
to be taken into account when multiple jets are the final product of a collision, different 
types of tree-level simulation descriptions should be used simultaneously: the matrix 
element approach and the parton-shower description. To this aim, modifications in 
the routines used by MadGraph/MadEvent were necessary. We have performed 
such modifications, and used the merge/matching technique to investigate graviton 
production in inclusive multi-jet samples at hadron colliders. 
A very general analysis has been carried out. Since the development of the routines 
used by MadGraph/MadEvent has been made aiming on the automatic generation 
of events involving spin-2 particles, we could investigate the differences/similarities 
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between three main classes of theories that contain the graviton particle. In the ADD 
and MGM theories, the graviton appears as missing energy, while in RS model it 
may be observed through its decay products. 

We have checked and validated the robustness of matching results with a comparison 
against next-to-leading order simulation of graviton + jet production. The fully inclu- 
sive samples have been normalized to the corresponding NLO inclusive total cross 
sections for producing G + X. Key distributions for the most relevant observables 
involving one or more jets have then been compared to those obtained by the NLO 
calculation G + jet + X. The overall agreement in shape as well as in normalization 
between the NLO observables and these obtained by making use of the matching 
technique is excellent. In particular, notice that a unique normalization factor (the 
one normalizing the overall inclusive sample) suffices to describe not only the shape 
but also the normalization of jet distributions. It demonstrates that the use of matched 
samples in experimental analyses is straightforward and accurate. We therefore make 
use of such technique once more in Chap. 6. Finally, we have seen that observables 
which are sensitive to the number of jets, such as the Hj, are better described by the 
matched sample than a fixed-order calculation. 

Since we are interested in the missing energy signature generated from theories that 
include gravitational interactions, supersymmetry within the context of simplified 
models of supergravity should be taken into account. Such theories comprise grav- 
itino particles, which are spin-3 /2 particles, super-partner of the graviton. In some 
variations of supersymmetric theories, gravitinos are super-light, invisible particles, 
whose productions can become relevant at hadron colliders. In Chap. 6 we thus con- 
sider the production of two gravitinos in multi-jet final states. 
To obtain the predictions for relevant experimental observables, we have made use 
of the matching technique, which takes into account simultaneously hard and soft 
jets. Because two different sub-processes separately contribute for our analysis, the 
use of the matching tool has been a major advantage. It has not only facilitated and 
prevented us from errors, but has also saved substantial computational time. That is 
due to the fact that we do not need to simulate each sub-process separately, nor worry 
about the region of the parameter space that a particular sub-process is dominant. 
Results obtained for gravitino production are encouraging. We have observed that, 
if such new physics signs are detected at the LHC, one could use observables related 
to the pr of the jets and/or pj to estimate the value for the mass of the gluino. Once 
this is fixed, jet multiplicities allow us to predict the mass of the gravitino, and the 
leading contributing sub-process. 

The methods and ideas employed in the course of the present research are at the cut- 
ting edge in high-energy phenomenology. We have helped to develop and advance 
tree-level simulation tools and search strategies to easily connect theory with experi- 
mental observables. With them, we have performed a complete analysis on processes 
involving missing energy signature for theories that incorporate the gravitational 
force. We now look forward to their use in the current and forthcoming searches at 
hadron collider. 
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7.2 Outlook 

The advancement of simulation methods and computer technology have defined the 
beginning a new era for the phenomenology of particle physics. Phenomenological 
investigations, are now essential for the search of BSM theories, and hints of their 
existence may appear very soon at the LHC. With the large number of possible 
SM extensions, it is necessary to know in advance the best strategy to identify new 
particles and connect them with the most favourable theory which these particles 
could belong to. 

For this reason, simulation tools and the search strategies developed and employed 
here can be very fruitful for future researches. First, the implementation of spin-2 
helicity routines can be used by MC event generators to simulate processes containing 
such particles. Several BSM theories can be considered, where similar analyses 
to those presented in Chaps. 5 and 6 can be performed. Moreover, phenomeno- 
logical investigations could be modified and extended to more model-independent 
approaches. 

The automation of helicity routines necessary for calculating the matrix element 
associated with any new model presents an indispensable step completed towards 
the effort to connect theory with experimental data. Any theory can now be eas- 
ily implemented into a MC event generator, accounting that the implementation of 
helicity routines are no longer restricted to the same Lorentz structures as these of 
the SM. Obtaining predictions for experimental observables (considering tree-level 
calculations) has become feasible for anyone, from a person with a more theoretical 
background to someone with a more experimental point of view. There are several 
theories that could be implemented and studied with the employment of such tool. 
For instance, there are extra dimensional theories with SM-like particles propagat- 
ing in the bulk, different SM extensions based on specific symmetries, as well as 
simplified models focused on the treatment of new particles. 

One example of the use of these tools would be a signature-based research aiming the 
discovery of a new particle. The idea would be to consider, in a model-independent 
manner, all the channels that could result in a particular signature. Similarly to the 
investigation that has been performed for the graviton production in Chap. 5, where 
we consider three representative classes of theories that could produce a graviton, 
we could concentrate in final states that include any type of invisible particles and 
build a strategy to differentiate theories using final observables. One goal would be 
the identification of the observable(s) that best discriminate classes of theories. 
A consequence of the non-observation of new physics so far, is the realization that 
some processes need more accurate predictions than these that can be provided by LO 
calculations. Since a few years ago, a huge effort is being made for the development 
of both the automatic NLO computations, and the MC event generators that take into 
account such corrections. Recent ideas, in addition to the advancement of simulation 
techniques and machines, enabled the creation of such NLO simulation tools, which 
are however still in (the final stage of) development [ 1 -3 ] . Once NLO event generators 
are available, an "accuracy revolution" will begin, and we will have the possibility of 
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performing several phenomenological investigations. In particular, we could extend 
the analysis made here for mono-jet graviton production at NLO, and obtain the 
prediction for the relevant observables with larger precision and, most importantly, 
the correct normalization for the cross section. 

Although NLO predictions provide accurate results for the cross section, it is still nec- 
essary to consider QCD emission and an hadronization description to correctly pre- 
dict observables at the detector level. The LO matching/merge description described 
in Chap. 4, and employed in Chaps. 5 and 6, is interesting because it considers the 
best features of both approaches. A novel description combining NLO event gener- 
ation and parton-shower description have been developed [4-9] and great progress 
has been made since their actual conception. Another idea would therefore be the 
use of such tool to perform similar investigations to these presented in the course of 
this work. With the use of such tools, one could obtain the distribution of relevant 
observables with the correct shapes and cross section normalization. 
Finally, more severe investigations will be necessary to constrain parameters of 
BSM theories with the LHC data once they are available. Simulation tools and 
search/identification strategies will help to accelerate this process and enable more 
theories to be looked into simultaneously. There is certainly still a long way ahead, 
when it comes to interpreting and understanding results that will soon be provided 
by the LHC, but that can also be the reason why the next years may be the most fun 
and exciting years for modern particle physics! 
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Appendix A 

Massless Graviton Implementation 



Here, we present in detail the implementation of massless graviton routines within 
MadGraph/MadEvent. Initially, we alter helicity routines (HELAS) for spin-2 
particles in order to handle massless particles. We subsequent create a model file 
for the particular theory described on Chap. 3. All information characterizing imple- 
mentation can be found in Ref. [1] and is reproduced in what follows. 



A.l HELAS Routines Modifications 



The HELAS routines for massive spin-2 particles have been implemented by K. 
Hagiwara et al. [2] in 2008. A summary of the routines can be found on Table A. 1 . In 
order to have massless spin-2 particles in MadGraph/MadEvent, we have modified 
the tensor wave function subroutine TXXXXX, and four subroutines for the off-shell 
tensor were added to the former implementation. 



A.1.1 TXXXXX 

This subroutine computes the graviton tensor wave function, namely e ltu (p, A) and 
e^ip, A)* in terms of the graviton four-momentum p and its helicity A. 
The helicity states of the graviton can be expressed by 

e^(p, ±2) = e^p, ±)e v {p, ±) 

e^{p, ±1) = i= [e?(p, ±)e>, 0) + e?(p, 0)e v (p, ±)] (A.l) 
e^ip, 0) = -J= [e^p, +)e v (p, -) + e^p, -)e v {p, +) + 2e»(p, 0)e>, 0)] , 
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Table A.l Summary of Spin-2 HELAS routines from [2] 



Vertex 


Inputs 


Output 


Subroutine 


SST 


SST 


Amplitude 


SSTXXX 




ST 


S 


HSTXXX 




ss 


T 


USSXXX 


FFT 


FFT 


Amplitude 


IOTXXX 




FT 


F 


FTIXXX, FTOXXX 




FF 


T 


UIOXXX 


VVT 


VVT 


Amplitude 


WTXXX 




VT 


V 


JVTXXX 




vv 


T 


UWXXX 


FFVT 


FFVT 


Amplitude 


IOVXX 




FVT 


F 


FVTIXX, FVTOXX 




FFT 


V 


JIOTXX 




FFV 


T 


UIOVXX 


VVVT 


VVVT 


Amplitude 


VWTXX 




VVT 


V 


JWTXX 




VVV 


T 


UWVXX 


VVVVT 


GGGGT 


Amplitude 


GGGGTX 




GGGT 


G 


JGGGTX 




GGGG 


T 


UGGGGX 



where e^ip, A) is the vector boson wave function, already pre-defined on the HELAS 
convention. For completeness it is important to state here that the spin-2 wave func- 
tions are traceless, transverse, orthogonal and symmetric, obeying therefore the fol- 
lowing relations: 

e%(p, A) = 0, e^ip, X)e liv {p, A')* = 5\x 

e^(p, A) = e^(p, A) Pfi e^(p, A) = p v e^(p, A) = 

The sum of helicity states for the graviton is introduced on the routines to be given 
by 

^ tuvip, A)e Q /3<J>. A)* - V^apip), (A3) 

A=±2 

where V/j, va p is the massless graviton propagator defined on Eq. (3.73). While the 
massive graviton has five helicity states of polarization, the massless spin-2 particles 
will only have two physical states. The sum of polarization for the massless graviton 
is therefore modified on the routines to set A = ± 1 , to zero allowing only A = ±2 
to contribute. 
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A.1.2 UIOXXX, UVVXXX, UIOVXX, uwvxx 



These subroutines compute off-shell tensor currents U by the FFT, WT, FFVT, and 
VWT vertices respectively. The main modification here was the inclusion of the zero 
mass graviton propagator to be the one showed in Eq. (3.73) for the massless case. 
The UIOXXX routine computes the bi-spinor spin-2 current from the fiowing-in and 
flowing-out fermions calculated by the FFT vertex. 

The flowing-in FI and flowing-out FO fermions are defined to be arrays which 
contain their four-momenta: 



FI 



/Pl\ 
P) 



. ™=( P o plplpl) 



(A.4) 



The output is a complex array defined by 



9t 



FO [-Vfiv(Pl+ Po - 2m j) + 7 M (p/ + p )y + j„(pi + p )fj] FI 



where q is the four-momentum related to the massless propagator, g t — ^ and A 
is the effective scale of the model. For the MGM theory, A = M p i(p t ) ~ p t . 
The UUVXXX is the subroutine that computes an off-shell spin-2 current from the 
two flowing-out vector bosons calculated by the WT vertex. The graviton output is 
the complex array: 

T a!3 = 9t—^— [(»»v + Pi ' Po)Ciivpo + D^pvipi, po) 

+ r 1 E^ P a(pi,Po)]v i p V^ (A.5) 

where where q is the four-momentum related to the massless propagator, g t = -j, A 
is the effective scale of the model, and the parameters C^ap, D^ap and E^ va p are 
defined on Table 3. 1 . 

The UIOVXX subroutine computes an off-shell spin-2 current calculated by the FFVT 
vertex. The graviton output complex array is given by: 

T a ^ = gt—j-to^tipc ~ <W) FO [gfPL + 9fPR)\ FI, (A.6) 

where q is the four-momentum related to the massless propagator, g, — Jr, A is the 
effective scale of the model, the parameter Cp_ va fi is defined on Table 3. 1 and <A and 
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(2) 

gy are the SM FFV couplings. For example, if we have a quark-anti-quark-gluon- 

(1) (2) 

graviton vertex, g *\ = g \ = g$. 

Finally, the UVWXX subroutine computes an off-shell spin-2 current calculated by 
the VWT vertex. The graviton tensor is given by 

T = — gv [Cfwpaipa ~ Pb)\ + C^axiPb ~ Pc)p 

+ C^Xp(Pc - Pa)a + F^x] V? V§ V c \ (A.7) 

where q is the four-momentum related to the massless propagator, g, — j-, A is 
the effective scale of the model, the parameters C^ va p and Fu,vpo\ are defined on 
Table 3.1 and gv is the SM VW coupling. 



A.2 MGM Model File 

The model file was created for MadGraph/MadEvent version 4. It contains all the 
information regarding particles, couplings and interactions and follows the structure 
presented on [2] for theories containing a massive spin-2 particle. A summary of the 
model content is presented below. 

The particle content of the mode is similar to the SM, in addition to a massless tensor: 

# This is a special data file which contains particles of the 

# model . The format for entering new particles is : 

# Name anti_Name Spin Linetype Mass Width Color Label PDG 

y y T D ZERO ZERO S G 39 

The graviton is represented by the letter y, and its PDG number is 39. 
The new coupling associated to the graviton particles are: 

# The ordering of particles is very important. 

# For FFV interactions, must use particles, not anti_particles . 
# 

# 3-particle vertices are entered as follows: 

# particlel particle2 particle3 coupling_name coupling_order 
# 

g g y GTV QTD 
a a y GTV QTD 
z z y GTV QTD 
w- w+ y GTV QTD 
q q y GTF QTD 
1- 1- y GTF QTD 
h h y GTS QTD 

# 

# 4-particle vertices are entered as follows: 
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# parti part2 part3 part4 2xcoupling_names 2xcoupling_orders 
# 

g g g y G GTV QCD QTD 
q q g y GG GTFV QCD QTD 
qd qd a y GAD GTFV QED QTD 
qu qu a y GAU GTFV QED QTD 
1- 1- ay GAL GTFV QED QTD 
qd qd z y GZD GTFV QED QTD 
qu qu z y GZU GTFV QED QTD 
1- 1- z y GZL GTFV QED QTD 
vl vl z y GZN GTFV QED QTD 
qd qu w- y GWF GTFV QED QTD 
qu qd w+ y GWF GTFV QED QTD 
vl 1- w+ y GWF GTFV QED QTD 
1- vl w- y GWF GTFV QED QTD 

Here, all types of quarks are represented by g, up-quark family is qu and down-quark 
family is qd. For the leptons, electron, muon and tau are defined as 1-, and the three 
types of neutrinos are represented by vl. The interaction order is the proportionality 
of the vertex with the g$ coupling for QCD, g e coupling for QED, and g, coupling 
for QTD, which is the new interaction defined for spin-2 particles. The couplings 
depend on the particles of the vertex. 

The new couplings related to the spin-2 particle are GTV, GTF and GTS. They are 
defined to be: 

1 1 1 

GTV= — , GTS = , GTF = (A.8) 

A 2A 4A 

where A is the effective scale of the theory. 
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New Routines for ADD Theories in MadGraph 



B.l 5-Point Vertices 

In Chap. 3 we have obtained the Feynman rules and vertices for a theory with a 
general spin-2 particle coupled to the SM theory. These are exhibited on Fig. 3.1. 
We can check the existence of 5-point vertices, such as the ggggG vertex or the 
qqggG. Previously, MadGraph/MadEvent version 4 [ 1 ] did not have the automatic 
generation of 5-point diagrams, which was necessary for the correct evaluation of 
processes with more than 1 jet in the final state. Routines to automatic generate 
5-point vertices have been implemented in collaboration with Q. Li [2]. Here we 
reproduce details of such implementation. 1 

Recall that the ggggG vertex is described on Fig. 3.1. It can be re-organized and 
written as: 

£ ggggT = - GTGC 2 f abe f de T"»* x [\r 1 ^9 a ' P *9 h ' <T *9'?9t* ~ 9? 9?], 

(B.l) 

with the coupling constants GT = — S and GC = g s . In HELAS such vertex can be 
encoded as GGGGTX computing the portion of the amplitude of the ggggT ver- 
tex from four gluon (G) polarization vectors and a tensor (T) boson wavefunction 
corresponding to the color structure f abe f cde ^ as it is illustrated on Fig. B. la. 
In order to implement this vertex in MadGraph, which supports at most 4-point ver- 
tices, we reinterpret the scattering in the t, u and s channels of the non-propagating 
octet tensor boson t& as an auxiliary particle. With ?a we can decompose the 5- 
point vertex into two 3-point vertices: ggtA and t/JhT , which can then be generated 
by MadGraph automatically. The color structure of the vertex f ahe f cde i s main- 
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Fig. B.l Illustration of the 
decomposition of the 5-point 
vertex into two 3 -point ver- 
tices to be automatically 
treated in MADGRAPH ver- 
sion 4 



Table B.l Summary of the HELAS sub-routines introduced for the implementation of the 5-point 
vertex in MADGRAPH extracted from [2] 

Sub-routine Summary 

WTCXX computes the amplitude of the ggf^t vertex 

JVTCXX computes the off-shell vector current J made from the interactions of a Vector 

gluon and an auxiliary Tensor boson by the gg?A vertex 

UWCXX computes the off-shell tensor current U for the auxiliary tensor , made from 

two flowing-out Vector gluons by the ggtA vertex 

TTTXXX computes the amplitude of the t A t A T vertex 

UTTAXX computes the off-shell non-propagating tensor current U, made from the 

flowing-out graviton Tensor and the auxiliary Tensor by the t\ t& T vertex 

UTTBXX computes the off-shell graviton tensor current U, made from the two flowing-out 

auxiliary Tensors by the tAtAT vertex 




tained. Figure B.l illustrates the decomposition. Observe that we assign a flow to tA, 
i.e. we treat tA different from its antiparticle in MadGraph. The idea is to forbid the 
interaction of tA and its antiparticle to avoid appearance of additional diagrams for 
99 ->• 99- 

In order to introduce the auxiliary particle tA, several HELAS sub-routines have 
been modified and some have been created. Table B.l presents a summary of all 
modifications. We follow with a full description of each routine. 



B.l.l GGGGTX 

The subroutine GGGGTX computes the portion of the amplitude of the ggggT vertex 
from four Gluon polarization vectors and a Tensor boson wavefunction corresponding 
to the color structure f abe f ccle ^ which should be called as 



CALL GGGGTX(VA, VB, VC, VD, TC, GC, GT, VERTEX) (B.2) 
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The inputs VA ( 6 ) , VB ( 6 ) and VC ( 6 ) are complex six-dimensional arrays which 
contain the Vector boson wavefunctions, and their momenta as 

p% = (3teVA(5), SfteVA(6), 3mVA(6), 3mVA(5)), 
pi = (9teVB(5), SfteVB(6), 3mVB(6), SmVB(5)), 
p£ - (9teVC(5), SfteVC(6), 3mVC(6), 3mVC(5)). 

The input TC ( 18 ) is a complex 18-dimensional array which contains the wavefunc- 
tion of the Tensor boson: 

r r y+ 1 ."+ 1 = TC (4/i + i/+l), (B.3) 

and its four-momentum as 

//' - (SReTC(17), sRgTC(18), 9fmTC(17), 9fmTC(17)), 

The output VERTEX is a complex number in units of GeV: 

VERTEX = -GT GC 2 G^, pX a5 V£VgV£V& (B.4) 

with G^na-xs given on Table3.1, and we use the notation 

V% = VAGH-l), vg = vboh-i), 
v£ = vcoh-i), v£ = vdoh-1). 

In order to insert automatically 5-point vertices in MadGraph, we introduce non- 
propagating colourless tensor boson t& as an auxiliary particle. With the portion 
of the 5-point vertex with color structure f abe f cde can be reduced automatically. 
After introducing particle Ia , several new HELAS subroutines were added. They are 
summarized on Table B.l, and can be called on the following way: 



B.1.2 VVTCXX 

This subroutine computes the amplitude of the gg?A vertex, 

CALL WTCXX(VA, VB, TC, GC, VMASS, VERTEX) 

VMASS represents the vector gluon mass which is zero (although it does not play 
any role here, we keep it as an input argument in accordance with the subroutine 
WTXXX, for convenience of MadGraph). What we compute here is 

VERTEX = GC V^VuaT^V^VS. (B.6) 



140 



Appendix B : New Routines for ADD Theories in MadGraph 



B.1.3 JVTCXX 

This subroutine computes an off-shell vector current J made from the interactions 
of a Vector gluon and an auxiliary Tensor boson by the ggf a vertex, and should be 
called as 

CALL JVTCXX(VC, TC, GC, VMASS, WIDTH, JVTC) 

The input VC ( 6 ) is the wavefunction and momentum of the gluon. The output 
JVTC ( 6 ) gives the off-shell vector current multiplied by the gluon propagator, 
which is expressed as a complex six-dimensional array: 

GC 

JVTC( + 1) = ^ atl V a T^, (B.7) 

<r 

and 

JVTC(5) = V(5) + TC(17), (B.8) 
JVTC(6) = V(6) + TC(18). (B.9) 

Here the momenta q are 

qf 1 = (3'leJVTC(5),3 f leJVTC(6),3OTjVTC(6),3mJVTC(5)). 



B.1.4 UVVCXX 

This subroutine computes an off-shell tensor current U for the auxiliary tensor tA, 
made from two flowing-out Vector gluons by the ggtA vertex, and should be called 

as 

CALL UWCXX(VA, VB, GC, VMASS, XM, XW, UWC) 

The inputs XM and XW are two dummy arguments for which we keep them in 
accordance with the subroutine UWXXX. The output UWC ( 18 ) is a complex 18- 
dimensional array: 

T a0 = GCV%vjj (B.10) 
for the first 16 components of UWC, and 

UWC(17) = VA(5) + VB(5), (B.ll) 



UWC(18) = VA(6) + VB(6). 



(B.12) 
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This subroutine computes the amplitude of the tA £a t vertex, called as 

CALL TTTXXX(TC, TIC, T2C, GT, VERTEX) 

The inputs TIC (18) and T2C(18) are complex 18-dimensional arrays which 
contain the wavefunction and momenta for the auxiliary tensors. What we compute 
here is 

VERTEX = GT 'T^Tf G „ A5 , (B.13) 
with 7*1,2 defined from TIC ( 18 ) andT2C(18) as in Eq. (B.3). 

B.1.6 UTTAXX 

This subroutine computes an off-shell non-propagating tensor current U, made from 
the flowing-out graviton Tensor and the auxiliary Tensor by the Za /AT vertex, and 
should be called as 

CALL UTTAXX (TC, TIC, GT, UTTA) 
The output UTTA ( 18 ) is a complex 18-dimensional array: 

Tf = GT n^T^T^G^xs (B.14) 
for the first 16 components of UTTA, and 

UTTA(17) = T(17) + T1(17), (B.15) 

UTTA(18) = T(18) + T1(18). (B.16) 

B.1.7 UTTBXX 

This subroutine computes an off-shell graviton tensor current U, made from the two 
flowing-out auxiliary Tensors by the tA tA^ vertex, and should be called as 
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CALL UTTBXX(T1C, T2C, GT, TMASS, TWIDTH, UTTB) 

The inputs TMASS and TWIDTH are the graviton mass and width, mj and Yj. The 
output UTTA ( 18 ) is a complex 18-dimensional array: 

T ap = gt - TC X TfG^ pc7XS (B.ll) 

q A — nij + imj\ j 

for the first 16 components of UTTB, and 

UTTB(17) = Tl(17) + T2(17), (B.18) 

UTTB(18) = T1(17) + T2(18). (B.19) 

Here q is the momentum of the off-shell tensor boson given in Eqs. (B.18) and (B.19) 
as 

= (9teUTTB(17), 3teUTTB(18), 3mUTTB(18), 3mUTTB(17)). 

And 

^a,i//S is the polarization summation tensor (3.77) [2]. 



B.2 Summation Over the Mass Spectrum 

On Chap. 2 it has been shown that for the ADD theory we expect to have an infinite 
number of spin-2 particles couples to the SM. Because their masses lie very close to 
each other experimental searches would hope to detect signs of gravity through the 
sum of such states. 

The MC simulation would therefore have to perform an extra integral (accounting for 
all the KK modes of the graviton) for each point of the phase-space. The amplitude 
integration over the mass density function is specified on Eq. (3.40). 
For virtual graviton exchange, the tower of KK graviton modes leads to the summa- 
tion of their propagators, which is divergent since ADD is only an effective theory. 
It is possible to introduce a cutoff for the highest KK modes and replace the sum- 
mation a factor on the effective scale, as it was first proposed by Han et al. [3]. This 
was the method chosen for the implementation in MadGraph, which can be easily 
performed by setting the graviton mass to zero and the extra dimensional scale to 
be the same as the effective scale. Note however, that the propagator of the massive 
graviton is still employed here, and setting the mass of the graviton to zero has only 
been accomplished locally aiming on executing such trick. 

Nevertheless, it is more complicated for graviton emission, because it is actually 
necessary to perform the integration described by Eq. (3.40). We have implemented 
the integration through the introduction of new vertices and two auxiliary dummy 
particles x\ and x%. These are massless colourless particles, in which the first is 
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Table B.2 Summary of new HELAS sub-routines introduced for the implementation of the KK 
mass integration in MADGRAPH/MADEVENT within the framework of the ADD model. Table 
extracted from [2] 

Sub-routine Summary 

PXXXXX stores the helicity and momentum of the auxiliary particles x\ and X2 

TPSXXX computes the amplitude of X1X2T vertex from a set of parameters specified by 

the theory 

UPSXXX computes an off-shell graviton tensor current U, made from the flowing-out 

auxiliary Scalar and Pesudo-particles, by the X1X2T vertex 



defined as a scalar, and the second a pseudo-tensor. MadGraph is subsequently 
modified to include a new particle type P, which is only involved in the interaction 
with graviton: xiX2 T - 

For generate a process (exclusively for the ADD theory) in which a spin-2 particle is 
emitted, one should request for example pp —> x\X2 + X, instead of pp — >• T + X. 
The phase spaces are related by 

dm\ 

d®(xix 2 + X) = %rd<i>(T + X). (B.20) 

I6n z 

We can thus perform the mass integration with the help of the new HELAS subrou- 
tines for the X1X2T vertex summarized in Table B.2. 



B.2.1 PXXXXX 

This subroutine stores the helicity and momentum of the auxiliary spin-2 Pseudo- 
particle for further usage, and should be called as 

CALL PXXXXX(P, XM, NHEL, NST, PC) 

The input P ( : 3 ) is a real four-dimensional array which contains the four- 
momentum p^ of the pseudo particle, NHEL (= ±2, ±1, 0) specifies its helicity 
A, NST specifies whether the boson is in the final state (NST = 1) or in the initial 
state (NST = -1). XM is its mass but does not play any role here. 
The output PC ( 1 8 ) is a complex 1 8-dimensional array, among which only the 
following matter, namely 

PC( 1) = NHEL, (B.21) 

and 

(PC(17), PC(18)) = NST (P(0) + z'P(3), P(l) + z'P(2)). (B.22) 
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B.2.2 TPSXXX 



This subroutine computes the amplitude of the X1X2T vertex, and should be called as 

CALL TPSXXX(TCO, PC, SC, GP, TMASSO, TWIDTH, VERTEX) 

SC ( 3 ) is a complex three-dimensional array which contains the momentum for x\ . 
GP(1:2) here are not coupling constant, but contain the ADD model inputs: 

GP(1) = A + iS, 

(B.23) 

GP(2) = M low + iM up , 

M] ow and M llp are the lower and upper limits for graviton mass integration. 
TMASSO (mjo for below) is only temporarily used and got from the MG input 
file param_card . dat, while the real graviton mass raj is defined as the invariant 
mass of xl and x2, i.e., 

m T = \ qi +q 2 \ t (B.24) 
with qi and 172 are the momenta of x\ and x%. 

Also note the input TC is indeed calculated from the other side of Feynman diagram, 
and contains the polarization summation tensor Va Vp(7 , which can be used to project 
out the real graviton wave function, due to the relations of the 5 helicity states of the 
tensor boson [4]: 

7V,a/?0>)= X ^,Ap,X)e a 3(p,X)*, (B.25) 

A=±2, ±1.0 



and 



^ v {p, A)e,«y(p, A')* = 5\\>, (B.26) 



with p = q\ + qi, and the helicity A is got from the pseudo-particle's wavefunction, 
orEq. (B.21). 

Now we can get the true graviton tensor wave function by calling the subroutine 

CALL TXXXXX(P, mj, A, +1, TC) (B.27) 
Finally, the output VERTEX is 

VERTEX = — p(m T )(p 2 - m 2 TQ + im TQ l WW^Tg 1 ' (B.28) 
raj 

x 9(m T - M, ow )6(M up - m T ), 
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where 6 represents the Heaviside step function. In equation (B.28), we include the 
mass density factor, the inverse graviton propagator, and the compensation factor for 
the decay phase space. 



B.2.3 UPSXXX 

This subroutine computes an off-shell graviton tensor current U, made from the 
flowing-out auxiliary Scalar and Pesudo-particles, by the X1X2T vertex, and should 
be called as 

CALL UPSXXX(PC, SC, GP, TMASSO, TWIDTH, UPS) 

The output UPS is a complex 18-dimensional array, which is indeed got simply by 
calling inside UPSXXX 

CALL TXXXXX(P, m T , A, +1, UPS). (B.29) 
and including the step functions in Eq. (B.28). 
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Appendix C 

ALOHA Automatic Libraries of Helicity 
Amplitudes for Feynman Diagram 
Computations 



Aloha is a package for the automation of the HELAS (HELicity Amplitude Subrou- 
tines) library [1] corresponding to the Feynman rules of any in quantum field theory 
Lagrangian. 

From the input given by the Universal FeynRules Output (Ufo) [2], Aloha produces 
the complete set of routines, wave-functions and amplitudes, that are needed for the 
computation of Feynman diagrams. The code is written in Python, but the output 
language can be chosen by the user between Python, C++ and Fortran. Here, 2 
we present how ALOHA works, and show a sample application implemented in the 
MadGraph 5 framework for theories containing spin-2. 



C.l Introduction 



We have seen that due to the numerous quantity of possible BSM theories that could 
show up at the LHC, extensive efforts are being made to automate the process of 
going from the Lagrangian of a new theory, to simulating events that can be then 
compared with data. Aloha is one important piece of that process with the ability to 
automatically create helicity amplitude routines necessary for calculating the matrix 
element associated with any new model. Such an ability will be useful to many groups 
writing custom generators, however, at present the greatest impact is on the recently 
released MadGraph 5 suite [4]. 

Helicity amplitudes methods [5-1 1 ] are apowerful way to evaluate the squared matrix 
element of any process. Several Monte-Carlo (MC) event generators, matrix-element 
based, use such routines already at the initial stage of the calculation. 



2 This appendix is a reproduction of excerpts and tables of [3], Reprinted from Computer Physics 
Communications, Vol. 183/10, Priscila de Aquino, William Link, Fabio Maltoni, Olivier Mattelaer 
and Tim Stelzer, ALOHA: Automatic libraries of helicity amplitudes for Feynman diagram com- 
putations, Pages No. 2254 - 2263, Licence numbers 3126460314522 and 3126460539382 (2012), 
with permission from Elsevier. 
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The HELAS [1] library has been a particularly successful implementation of an 
helicity amplitude method for tree-level processes. Amplitudes are generated by ini- 
tializing a set of external wavefunctions using their helicity and momenta. These 
wavefunctions are combined based on the particle interactions in the Lagrangian 
to determine the wavefunctions of the internal lines (propagators). Once all of the 
wavefunctions are determined, they are combined to calculate the complex number 
corresponding to the amplitude for the diagram. These amplitudes can then be added 
and squared to give the required result. The set of routines needed to calculate ampli- 
tudes for the SM at tree-level was released in the original HELAS and successfully 
employed for many phenomenological studies thereafter. 



C.2 Automation of the Method 

Aloha basically uses the same strategy as these of HELAS. These routines consists 
of three categories: 

External particle routines: Wavefunction routines corresponding to each external 
leg of a Feynman diagram. They are computed for a given momentum and helicity. 
The type of the particle must be given as an input. The choices are: vector, scalar, 
tensor or incoming/outgoing fermion. 

Wavefunction/Off-shell routines: The propagator wave function, which is gener- 
ated when one leg of a vertex has no associated external wave-function. These 
routines return both the internal wavefunction and the momentum of the propa- 
gator. 

Amplitude/On-shell routines: The amplitude routines, which is computed after all 
the legs of a diagram are associated to a wave-function. These can be used to 
compute the matrix-element of the diagram by the associated amplitude routines. 

To illustrate how Aloha works from the beginning, let us start with the one example 
of vertex produced by the Ufo. The Ufo output for the e + e~j vertex is: 

FFV = Lorentzfname = 1 FFV ' , 

spins =[2,2,3], 
structure = ' Gamma ( 3 , 2 , 1 ) ' ) 

The spins attribute contains the list of spins of each particle (present in the vertex) 
in the 2 *S+1 scheme. 3 The structure property collects the analytical expression 
for the Lorentz structure in the vertex. The arguments of the Lorentz object refer to 
the particle associated with each index of the object. The list of available objects and 
the Python convention for different indices are explained in Table C.l. 
The Ufo convention requires that: 



In the case of a fermion, the first particle is considered as an incoming fermion and the second 
as an outgoing fermion. For interactions containing more than two fermions, the convention is that 
the vertex should be defined for a given fermion flow. The outgoing fermion follows directly its 
associated incoming fermion. 
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Table C.l The greek indices 
stands for lorentz indices 
while the lower case latin 
indices stands for spin 
indices. The indices N stands 
for the label of the particle 



Elementary Lorentz structures 



Charge conjugation matrix: C a b 
Epsilon matrix: t >lvf "' 
Gamma matrices: (a^ab 
Gamma5 matrix: (cr')ab 
(Spinorial) Kronecker delta: 5 a & 
Minkowski metric: 
Momenta of particle N: p^, 
Right chiral projector: ( -4?- ) 

V J ab 

Left chiral projector ( —S— ) 

V / ab 

Sigma matrices: (a tlv ) a b 
Scalar wavefunction: (f>^ 
Spinor wavefunction: tp^ 
Vector wavefunction: 
Spin2 wavefunction: Tjt' 



C(a,b) 

Eps i lon(fj,,v,p,a) 
Gamma(/i, a, b) 
Gamma5 (a, b) 
Identity (a, b) 
Metric(p,^) 
P(//,N) 

Proj P (a, b) 
ProjM (a, b) 

Sigma(/i,^,a, b) 

Scalar(N) 

Spinor(a,N) 

Vector(p,N) 

Spin2(/j,v,N) 



i. positive indices are linked to the particle numbers; 

ii. all Lorentz indices should be listed before spin indices; 

iii. all repeated indices should be summed with the appropriate matrix; 

iv. and all momenta should be incoming. 

From a given Ufo model structure, ALOHA creates the analytical expression linked 
to the associated HELAS routine by contracting the expression with a set of wave- 
functions. 

Now, we exemplify how ALOHA works by relating a particular vertex produced by 
ALOHA and HELAS. The structure shown in the FFV python object introduced above 
is Gamma (3,2,1). 



C.2.1 Amplitude/On-shell Routine 

The associated amplitude routine, which was previously named IOVby the HELAS con- 
vention, is now: 

Spinor (2) * Gamma ( 3 , 2 , 1 ) * Spinor (1) * Vector (3) 
which correspond to the analytical expression: 

FFV_00/>i,V>2,e3 A1 ,sO = g*U27 fl v l e fl . (C.l) 

The analytical module of ALOHA expands Eq. (C.l) component by component. The 
result is written to a C++, Python or Fortran file, language chosen by the user. 
For instance, the Fortran output would be: 
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C This File is Automatically generated by ALOHA 

C The process calculated in this file is : 

C Gamma (3,2,1) 

C 

SUBROUTINE FFV1_0 ( F 1 , F2 , V3 , COUP , VERTEX ) 

IMPLICIT NONE 

DOUBLE COMPLEX Fl(*) 

DOUBLE COMPLEX F2 ( * ) 

DOUBLE COMPLEX V3 ( * ) 

DOUBLE COMPLEX COUP 

DOUBLE COMPLEX VERTEX 

VERTEX = COUP* ( (F2 (1) * ( (Fl (3) * ( (0, -1) *V3 (l) + (0, 1) *V3 (4) ) ) 
S +(F1(4)*((0, 1)*V3(2)+V3(3)))))+((F2(2)*((F1(3)*((0, 1) 
$ *V3 (2) -V3 (3) ) )+(Fl (4) * ( (0, -1) *V3 (l) + (0, -1) *V3 (4) ) ) ) ) 
$ +( (F2 (3) * ( (Fl (1) * ( (0, -1) *V3 (1) + (0, -1) *V3 (4) ) )+(Fl(2) 
$*((0, -1) *V3 (2) -V3 (3) ) ) ) ) + (F2 (4) * ( (Fl (1) * ( (0, -1)*V3(2) 
$ +V3(3))) + (F1(2)*((0, -l)*V3(l) + (0, 1)*V3 (4) )))))) ) 

END 

The convention for the naming of Aloha's routines is denoted by NAME_X, where 
NAME is the name provided by Ufo and X is either for an amplitude routine or the 
number corresponding to the off-shell particle for a wavefunction routine. 



C.2.2 Wavef unction/Off -shell Routines 

To produce the necessary off-shell routines, Aloha loops over the list of particles and 
replaces the on-shell wavefunction with a propagator. For particle 1, the incoming 
spinor is replaced by a propagator, computing 4 : 

Spinor(2) * Gamma ( 3 , 2 , 1 ) * SpinorPropagator ( 1 ) * Vector (3) 
which is analogous to: 

- p -\- ffii 

FFV_1(V>2, e 3M , 5, mi, Ti) = 5*Vf27" ■ ' , 2 : — W (C2) 

p L — m\ + im\\ i 

FFV_1 is the wavefunction of the first spinor, and would be FOV in the original 
HELAS language. The momentum for the propagator is computed by using conser- 
vation of momentum at the vertex. 
Similarly, for the outgoing spinor, we have: 

FFV_2 (?/>!, e 3 „, g, m 2 , V 2 ) = i . { 7 ^ 2 ^ (C .3) 

p A — rrij + im,2\ 2 



4 The propagator object is defined internally to ALOHA. 
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The syntax used by Aloha reads: 

SpinorPropagator ( 2 ) * Gamma ( 3 , 2 , 1 ) * Spinor(l) * Vector (3) 

which can be related to the FIV routine in the original HELAS language. 
Finally, Aloha also creates the routine associated to the off-shell vector. This rou- 
tine returns a vector wavefunction from a pair of incoming and outgoing fermion 
wavefunctions, and it is written by ALOHA as: 

Spinor(2) * Gamma ( 3 , 2 , 1 ) Spinor(l) * VectorPropagator ( 3 ) 

analogous to the analytic expression FFV_3, which would be JIO in the original 
HELAS language: 

FFV_3(V>i, ip 2 , g, »»3. r 3 ) = ip 2 i In ipi ■ -i — ; j r ^ ■ ( c - 4 ) 

p- — m\ + imji 3 

Following the original HELAS convention, it is used the unitary gauge for massive 
particles and the light cone gauge for massless particles. 



C.2.3 External Particle Routines 

Aloha is not designed to create external particles routines, e.g., IXXXXX and 
OXXXXX for incoming and outgoing fermions, VXXXXX for vectors, SXXXXX for 
scalars and TXXXXX for spin-2 particles. These from HELAS are used. There is a 
relatively small number of such routines needed, and for this reason is reasonable to 
use a static library. This is incorporated into the ALOHA package, and it is available 
in Python, C++ and Fortran languages. 

Two additional types of special routines are needed in order to cover the full range 
of routines required for any type of interaction (and model): 

• Conjugate Routines: Required by the presence of majorana fermions and of 
fermion flow violation [12]. In the case of the previous example, the Lorentz 
expression would have to be replaced by: 

C(-2,2) * Spinor(2) * Gamma (3, 2,1) * 

Spinor(l) * Vector(3) * C(3,-3) 

To distinguish conjugate routines from the standard ones, ALOHA adds to their 
name the letter C followed by the number of the fermion pair that is conjugated. 
In this case the name of the conjugate routine would be FFVC1_X with X being 
either , 1 , 2 or 3 depending on the number corresponding to the amplitude 
or off-shell particle. Note that_l (_2) returns the wavefunctions associated to the 
incoming (outcoming) fermion. 
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• Multiple Couplings: In the Ufo scheme each Lorentz structure is associated to a 
single coupling. In practice, this means that some interactions are linked to more 
than one Lorentz structure. However, it is convenient for the readability and the 
compactness of the code to have one single routine per interaction. Therefore, 
Aloha defines a series of wrapper functions which call different routines associ- 
ated by their couplings. The name of these wrapper routines are the concatenation 
of the name of the associated routines combined by an underscore. If all of the 
associated routines start with the same prefix (with letter S, F, V, T) then this pre- 
fix is only kept the first time. For example, the wrapper function of FFV1_X and 
FFV2_X will be called FFVl_2_X. 

In some cases, this method provides redundancy. Consider a four-Higgs-interaction. 
Due to the symmetry of the Lorentz structure, the four wavefunction routines will 
each have the same content. This type of symmetry is detected automatically by 
Aloha and the computation is done only once. However, in order to simplify the 
work of the matrix element generator, all four routines are explicitly defined and 
have their own names. The first one is computed and the remaining three are written 
as an alias of the first one. 



C.3 ALOHA Output 

The major advantage of ALOHA comes from using it in combination with a MC 
event generator. To this aim, ALOHA offers two options. The first is to use ALOHA as 
a standalone program. The user can call Aloha from the shell: 

./bin/aloha UFOjath [ -f Fortran | Python | CPP] [-o output_dir] 

where UFO_path is the path to the UFO model directory, and the two possible 
arguments are 

• the language in which the routine should be written (default is Fortran). 

• the path for the output directory (which is by default UFO_path/LANGUAGE, 
where LANGUAGE can be either Fortran, Python or C++). 

The second possibility is to internally link the Aloha code to a MC event generator. 
This option is only possible for a generator written in Python, and is currently used 
by MadGraph version 5. For each process MadGraph requests a set of helicity 
routines, which is produced by Aloha. This first implementation of Aloha inside 
a matrix element generator played an essential role in fully validating the package. 



C.4 Checks and Validation 

To validate Aloha, we use two different versions of MadGraph to compute values 
for the square of the matrix element, given a specific phase-space point of a known 
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process. We employ MadGraph version 4.4.44, which uses the original HELAS rou- 
tines, and compare the results with those from MadGraph version 5.1.1.0, which 
uses the new Aloha routines. 

The comparison is made by specifying a point of phase space and a center-of-mass 
energy, after fixing properly the set of external parameters. The values for each matrix 
element squared are computed, and the difference is given by AD: 

AD = 2 x \ MiAh0HA ^ 2 ~ |X(HELAS)| 2 
X |7W(ALOHA)| 2 + |A1(HELAS)| 2 ' 

Due to effects of numerical precision, one can not expect exact agreement between 
the two codes. Instead, the test is considered to pass if values of the squared matrix 
element computed by Aloha and HELAS agree to at least five orders of magnitude, 
i.e., AD < 1 x 1(T 5 . 

Such tests are carried out for different processes based on the Standard Model (SM), 
the Minimal Supersymmetric Model (MSSM) and the Randall-Sundrum type I (RS- 
I) model, the latter including spin-2 particles. Additional tests are made by checking 
gauge and Lorentz invariance of each process. We observe a great agreement for all 
processes, as exemplified in Table C.2 for the SM and the MSSM. 



C.4.1 Spin-2 Validation Checks 

We have seen that several models with extra dimensions have recently been proposed. 
Examples of such theories include the Large Extra Dimensional model (or ADD) 
[13-15] and Randall-Sundrum (RS) models [16, 17]. These models contain spin-2 
particles allowed to propagate through higher dimensional space-time that couple 
to the SM particles in a very particular way. To complement the implementation 
of the ALOHA helicity routines, it is interesting to test and validate spin-2 routines 
comparing it to the HELAS implementation [18]. 

Using FeynRules, we implement a reduced version of the RS model [19]: it contains 
only the first mode of the spin-2 particle explicitly defined (which is denominated y), 
which couples to SM particles. Validation tests, similar to the ones presented above 
for the SM and MSSM, are performed. Approximately 500 2 ^ 2 and 2000 2^3 
processes are tested, and some examples can be found in Table C.3. We obtain a full 
agreement between routines obtained with the use of HELAS and ALOHA. 



C.4.2 Goldstino Validation checks 

Similarly to what has been performed for the SM, MSSM and RS theories, a sim- 
plified SUSY theory that includes goldstino interactions has been implemented in 
MadGraph/MadEvent via the FeynRules package [20]. Since the goldstino is 
the helicity ±1/2 component of the gravitino, it can be treated by ALOHA. A com- 
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Table C.2 Example of the results obtained by the validation tests for SM and MSSM processes 



SM processes 



MSSM processes 



Process 



v. 



77 

7T -> i/, W~ 
b e~ -> b e~ 

M ^ J/ 



7 e + 



„+ 



v e W~ — > ~i e~ 
v e Z — > v e Z 
19 -*■ tt 
gt-+tZ 
tt -> W + W~ 
b W~ -s- b W~ 

W+ w~ z z 

Z Z -» e+ e~ 
e~ W~ -> e- W" : 
f c Z — ► 7 e~ W + 
i f v e ^ "f uv e 
tu e ^e + t W~ 
v e W+ ->■ bb e+ 
7 -> e + e~ 
» Z — ► e + e~ h 
d d -*■ H v e v e 
HZ -» pf ; 
r ; — >■ w « Z 

; ->• e~ f i> e 
£>f -> 7&f 
/ Z ^ ddt 
b Z^ b ZZ 
W+ W~ ^ dd H 



AD 

1.512 x 1(T 7 
4.541 x 10~ 8 
4.545 x 1(T 8 
7.310 x 10~ 8 
3.911 x 1(T 8 
5.458 x 10~ 8 
7.164 x 1CT 8 
4.545 x 10~ 8 
< 1.0 x 10- 10 
7.528 x 10~ 8 
3.173 x 10~ 8 
1.833 x 10~ 8 
6.291 x 10~ 9 
1.666 x 10~ 7 
1.162 x 10~ 7 
9.767 x 10~ 8 
2.692 x 10~ 8 
1.179 x 10~ 7 
8.095 x 10~ 8 
8.430 x 10~ 8 
6.045 x 10~ 8 
3.816 x 10~ 8 
2.566 x 10~ 8 
2.457xl0- 8 
3.962 x 10~ 8 
4.744 x 10~ 8 
7.409 x 10~ 8 
2.224 x 10~ 8 
8.196 x 10~ 8 
9.519 x 10~ 8 



Process 



W 4 



4 

H 2 

e+ W- 
v e v e - 
e+v e - 



-> 4 W+ 
-+ H 3 W+ 
-+e+ H- 
H+ H~ 



- R 



e L X l_ 

sv e Xj — > u e W~ 
H~ sv e — ► v e xj~ 
e+ H~ -> e+ W~ 

H2 X\ -> e~ su e 

4 

urW -> guitW 
n\ Z — » ur u Z 

7 — > e~ su e v e 
e~ xf -» v e W + xj~ 
sv e — > e + v e v e 



Hi v e 



„+ 



n\ sv e — > v e v e v e 
gu L -? e~ e\ u 
go go ->■ Ur ur Z 
go ur u Z Z 
uv e — > g sv e Hi 



e + u L Z 
Mr -> e~ gu 
e~ — >• n\uu 
v e v e — > e\ v e x^ 



AD 

5.544 x 
1.108 x 
5.285 x 
6.214 x 
5.987 x 
1.792 x 
6.679 x 
1.984 x 

3.028 x 
3.714 x 
1.039 x 
7.992 x 
2.863 x 
1.527 x 
6.045 x 

5.029 x 
8.805 x 
4.001 x 
3.355 x 
6.548 x 
1.116x 
2.176 x 
3.677 x 
1.089 x 
1.526 x 
1.233 x 
1.542 x 
2.204 x 
2.122 x 
4.661 x 



10- 8 
10- 8 

10- 9 
10- 9 

10- 8 

10- 9 
10- 9 
10- 9 

10- 8 
10- 8 
10- 8 

10- 9 

10- 8 

10- 7 
10- 7 
10- 9 

10- 10 

10- 9 
10- 7 

10- 8 
10 -io 

10- 9 
10- 10 

10- 7 

10- 10 

10- 8 

10- 8 

10- 7 

10- 8 

10" 8 



parison of the values of the matrix element squared obtained by Aloha against these 
obtained with the use of HELAS [21, 22] has been performed. A full agreement is 
again obtained. 

The next step is to add the treatment of spin-3/2 particles by ALOHA, which is 
not yet completed. However, since in Chap. 6 we have shown that goldstino parti- 
cles are the relevant components of gravitino production at the LHC, for the phe- 
nomenological studies of gravitino emission (at the LHC) it is enough that Mad- 
Graph/MadEvent handle goldstino particles. 
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Table C.3 Example of the 
results obtained by the 
validation tests for RS 
processes 



RS processes 



Process 


AD 


Q U t U V 


< 1.0 x 10 -10 


b y — > b g 


< 1.0 x 10 -10 


t y ^ Qt 


< 1.0 x 10~ 10 


y y 


< 1.0 x 10~ 10 


a t v t v 

y i — ?- i v 


< 1.0 x 10 -10 


a a — > o v 

as j j 


< 1.0 x 10~ 10 


ty -*■ ggt 


9.766 x 10" 10 


gg-> ggy 


3.979 x 10- 10 


bg-*bg_y 


8.009 x 10-'° 


bb bby 


1.864 x 10~ 9 


hi ' -*■ biy 


7.411 x 10" 10 


tt-+ bby 


1.928 x 10" 10 
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